Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,757,942 C→A 5.3% intergenic (‑115/+75) yiaY ← / ← selB L‑threonine dehydrogenase/selenocysteinyl‑tRNA‑specific translation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,757,9420CA 5.3% 408.3 / 9.0 141intergenic (‑115/+75)yiaY/selBL‑threonine dehydrogenase/selenocysteinyl‑tRNA‑specific translation factor
Reads supporting (aligned to +/- strand):  ref base C (74/59);  new base A (4/4);  total (78/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 145 total aligned reads displayed.

TTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTC  >  NC_000913/3757793‑3758091
                                                                                                                                                     |                                                                                                                                                       
ttCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTc                                                                                                                                                         >  1:1420079/1‑150 (MQ=255)
 taacaGAAGGTATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGTGCAGCTTACTGATGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCa                                                                                                                                                        <  1:1729120/148‑1 (MQ=255)
  cacaGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCAt                                                                                                                                                       >  2:1531503/1‑150 (MQ=255)
  cacaGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCAt                                                                                                                                                       >  2:813819/1‑150 (MQ=255)
     aGAAGGAATAAAGAACGTTGAAGCTGCCATACCACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCGTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATATAATAAAAATCACATTGTATTCATTaa                                                                                                                                                    >  2:1277876/1‑150 (MQ=255)
      gAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTaa‑a                                                                                                                                                  >  2:71736/1‑150 (MQ=255)
         ggAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACaa                                                                                                                                               >  2:819031/1‑150 (MQ=255)
           aaTAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAAtg                                                                                                                                             <  2:897512/150‑1 (MQ=255)
           aaTAAAGAACGTAGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGTGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTAACAAAAATAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAAtt                                                                                                                                             >  2:1585971/1‑149 (MQ=255)
            aTAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAAtgt                                                                                                                                            >  1:567735/1‑150 (MQ=255)
            aTAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAAtgt                                                                                                                                            >  1:1646162/1‑150 (MQ=255)
             tAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAAtgtg                                                                                                                                           <  2:620359/150‑1 (MQ=255)
                 gAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑ttt                                                                                                                                      >  1:695323/1‑150 (MQ=255)
                 caaCGTTGAATCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCTTGCAGATTACTGAGGATTTTCAATCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGT‑TTAATGAATACAATGTGC‑ttt                                                                                                                                      <  1:270443/149‑1 (MQ=255)
                    cGTTGAATCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATTAAAATCACATTGT‑TTAATTAATACAATGTGC‑ttttta                                                                                                                                   <  2:452078/150‑2 (MQ=255)
                    cGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑tttttt                                                                                                                                   <  1:1597630/150‑1 (MQ=255)
                     gTTGAAGCTGCCATAACACACTCCTGAAATTGTAAAAAATGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑ttttttt                                                                                                                                  <  2:1414130/150‑1 (MQ=255)
                     gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑ttttttt                                                                                                                                  >  1:448347/1‑150 (MQ=255)
                     gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑ttttttt                                                                                                                                  <  1:1130128/150‑1 (MQ=255)
                     gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑ttttttt                                                                                                                                  <  1:388416/150‑1 (MQ=255)
                        gAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTaaa                                                                                                                               <  2:656623/150‑1 (MQ=255)
                         aaGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAAt                                                                                                                              <  2:2018224/150‑1 (MQ=255)
                          aGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAAtt                                                                                                                             >  1:1678589/1‑150 (MQ=255)
                           gCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTc                                                                                                                            <  1:1515277/150‑1 (MQ=255)
                               ccATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATaa                                                                                                                        >  1:178815/1‑150 (MQ=255)
                                  tAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGaa                                                                                                                     >  2:232656/1‑150 (MQ=255)
                                    acacacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAAtt                                                                                                                   >  1:792344/1‑150 (MQ=255)
                                       cacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTg                                                                                                                <  2:1346914/150‑1 (MQ=255)
                                       aacTCCTGAAAGTGTAAAAAAGTGCGGTCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTg                                                                                                                <  1:984489/149‑1 (MQ=255)
                                           ccTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTAtt                                                                                                            >  1:872538/1‑150 (MQ=255)
                                              taaaGTTTAAAAAAGGGCGGGAAGCTTACTGATGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAat                                                                                                         <  1:202271/149‑1 (MQ=255)
                                              taaaGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTAGAAAATATGATCATTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAat                                                                                                         <  2:1617781/149‑1 (MQ=255)
                                              gAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTTA‑ACACTGTGC‑TTTATTTAAATTCATAAGAATTTTGTATTTat                                                                                                         >  1:360023/1‑150 (MQ=255)
                                               aaaTTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAatt                                                                                                        <  1:1262765/150‑1 (MQ=255)
                                                  gtgtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAatttat                                                                                                     >  1:1065914/1‑150 (MQ=255)
                                                  gtgtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAatttat                                                                                                     >  1:510430/1‑150 (MQ=255)
                                                       aaaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaa                                                                                                >  1:1724338/1‑150 (MQ=255)
                                                        aaaaGGGCTGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaaa                                                                                               <  2:235612/150‑1 (MQ=255)
                                                        aaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaaa                                                                                               <  1:713788/150‑1 (MQ=255)
                                                          aaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAAc                                                                                             <  1:1817387/150‑1 (MQ=255)
                                                          aaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAAc                                                                                             <  2:447468/150‑1 (MQ=255)
                                                                         aCTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTc                                                                              <  1:1967239/150‑1 (MQ=255)
                                                                          cTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATAAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTcc                                                                             >  2:1823857/1‑150 (MQ=255)
                                                                           tGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCAGATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCt                                                                            <  2:480619/150‑1 (MQ=255)
                                                                            gAGGATTTTGATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACTATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCtt                                                                           <  2:794485/150‑1 (MQ=255)
                                                                              ggATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAt                                                                         <  1:1257788/150‑1 (MQ=255)
                                                                               taTTTTAATTAGAAAATATTAACATTATTATTGATTATCAAAATTAATCTAATAAAAATAACATTGTATTAATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTTTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAtt                                                                        <  1:1810344/149‑1 (MQ=255)
                                                                               gATTTTCATTAGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAAAAAAAATCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTATATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAtt                                                                        <  1:1818516/150‑1 (MQ=255)
                                                                                 ttttCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTTTTAATTTATTTTAAAACTAAATTAATCATTTCATTAtttt                                                                      >  2:317161/1‑150 (MQ=255)
                                                                                  tttCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAttttt                                                                     >  2:912365/1‑150 (MQ=255)
                                                                                    tCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTcc                                                                   <  2:723277/150‑1 (MQ=255)
                                                                                    tCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTcc                                                                   <  2:158979/150‑1 (MQ=255)
                                                                                      aTTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                                                                 >  1:693103/1‑150 (MQ=255)
                                                                                       ttCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                                                                 >  1:2015358/1‑150 (MQ=255)
                                                                                          gAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaat                                                             >  1:788971/1‑150 (MQ=255)
                                                                                                  gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGATTTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaataataattc                                                     >  1:400573/1‑148 (MQ=255)
                                                                                                  gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATga                                                     >  1:427506/1‑149 (MQ=255)
                                                                                                  gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGc                                                     >  2:1295223/1‑150 (MQ=255)
                                                                                                  gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGc                                                     >  2:1861844/1‑150 (MQ=255)
                                                                                                  gATCAGGATTATTGATAATCAAAATTAATCTAATAAAAAGCACAGGGGATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATATGAATTTTGTATTAATTTATTTTAAATCTGAATTAAGCATTTCCTTATTTTTCCGGATATAATAATGc                                                     >  2:295617/1‑150 (MQ=255)
                                                                                                           tattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCaa                                            <  1:911004/150‑1 (MQ=255)
                                                                                                           tattGATTATCAAAAATAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCaa                                            >  2:198372/1‑150 (MQ=255)
                                                                                                            attGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAAt                                           >  2:1337515/1‑150 (MQ=255)
                                                                                                              tGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATaa                                         >  2:1440723/1‑150 (MQ=255)
                                                                                                                aTTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAt                                       <  2:1264862/150‑1 (MQ=255)
                                                                                                                   aTCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGAt                                    >  1:1983920/1‑150 (MQ=255)
                                                                                                                   aTCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTATCCGGAAATAATAATGCGTCGCGCAATAAATGAt                                    <  2:1502199/150‑1 (MQ=255)
                                                                                                                    taaaaaTTAATATAATAAAATGAAAATTTTATTAATTAA‑AAAATTTGA‑TTTTTTAAAATTCATAACAATTTTGTATTAATTAATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATc                                   <  1:166164/148‑1 (MQ=255)
                                                                                                                    tCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATc                                   <  2:883281/150‑1 (MQ=255)
                                                                                                                    tCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATc                                   >  2:139139/1‑150 (MQ=255)
                                                                                                                    tCAAAATTAATCTAATAAAAAGCAAATTGTTTTCATTAA‑ACAATTTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCTGAAATAATAATGCGTCGCGCAATAAATGATc                                   <  2:446637/150‑1 (MQ=255)
                                                                                                                     cAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCa                                  >  1:1196740/1‑150 (MQ=255)
                                                                                                                     cAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAAATTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCa                                  >  2:508843/1‑150 (MQ=255)
                                                                                                                      aaaaTTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCAt                                 >  1:523594/1‑150 (MQ=255)
                                                                                                                        aaTTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTCCTGAAATAATAAAGAGTCGATCAATAAATGATCAttt                               >  2:360286/1‑150 (MQ=255)
                                                                                                                               cTAATTAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                        <  1:1021476/150‑1 (MQ=255)
                                                                                                                               cTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                        >  2:1964015/1‑150 (MQ=255)
                                                                                                                               cTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                        >  2:347566/1‑150 (MQ=255)
                                                                                                                               cTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                        >  2:4119/1‑150 (MQ=255)
                                                                                                                               cTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATAAATCATTTCCTTATTTTTCCGGAAAAAATAATGCTTCGCGCAATAAATGATCATTTCCacgac                        >  1:1794713/1‑150 (MQ=255)
                                                                                                                                taataaAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgacg                       >  2:870315/1‑150 (MQ=255)
                                                                                                                                taataaAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCAGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgacc                       >  2:1115799/1‑149 (MQ=255)
                                                                                                                                 taaaaaaaGCACATTGT‑TTAATGAATACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCTCAATAAATGATCATTTCCacgacga                      >  2:1726954/4‑150 (MQ=38)
                                                                                                                                 aataaAAAGCAAAATGTATTCATTAA‑ACAATGAGC‑TTTTAATAAATTAATAAGAATATTGTATTAATTTATATTAAAACTTAATTAATCATTTCCTTATTTTTCCGGATAAAATAATGCGTCGCGCAATAAATGATCATTACCacgacga                      >  1:1670670/1‑150 (MQ=255)
                                                                                                                                  ataaAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgacgac                     <  1:402945/150‑1 (MQ=255)
                                                                                                                                    aaaaaGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATTCGTCGCGCAATAAATGATCATTTCCACGACGAcgc                   >  2:1013489/1‑150 (MQ=255)
                                                                                                                                    aaaaaGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGAcgc                   >  2:744625/1‑150 (MQ=255)
                                                                                                                                     aaaaGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGAcgcg                  >  2:1658709/1‑150 (MQ=255)
                                                                                                                                      aaaGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCTGAAATAATAATGCGTCGCGCAATAAATCATCATATCCACGACTAcgctt                 >  1:1235687/1‑148 (MQ=255)
                                                                                                                                      aaaGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGt                 >  2:1799288/1‑150 (MQ=255)
                                                                                                                                      aaaGCACATTGT‑TTAATGAATACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGt                 >  1:1242682/1‑150 (MQ=255)
                                                                                                                                       aaGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTa                <  2:510430/150‑1 (MQ=255)
                                                                                                                                         gCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACTACGCGTaaa              >  2:190331/1‑150 (MQ=255)
                                                                                                                                          cacaTTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAg             >  1:1189075/1‑150 (MQ=255)
                                                                                                                                                ttATTCATTAT‑ACAATGTGA‑TTTTTTTAAATTCTTAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATGATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATg       <  2:1126598/149‑1 (MQ=255)
                                                                                                                                                gTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATg       >  2:1820424/1‑150 (MQ=255)
                                                                                                                                                 tATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGc      <  1:561117/150‑1 (MQ=255)
                                                                                                                                                    tCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGt   <  2:1422767/150‑1 (MQ=255)
                                                                                                                                                     cATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTc  <  2:178815/150‑1 (MQ=255)
                                                                                                                                                     cATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAAAAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGAggaa  >  2:822346/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTC  >  NC_000913/3757793‑3758091

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑