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.

TTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTC  >  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)
      gAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTaaa                                                                                                                                                  >  2:71736/1‑150 (MQ=255)
         ggAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACaa                                                                                                                                               >  2:819031/1‑150 (MQ=255)
           aaTAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtg                                                                                                                                             <  2:897512/150‑1 (MQ=255)
           aaTAAAGAACGTAGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGTGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTAACAAAAATAATCTAATAAAAAGCACATTGTATTCATTAAACAAtt                                                                                                                                             >  2:1585971/1‑149 (MQ=255)
            aTAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtgt                                                                                                                                            >  1:567735/1‑150 (MQ=255)
            aTAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtgt                                                                                                                                            >  1:1646162/1‑150 (MQ=255)
             tAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAtgtg                                                                                                                                           <  2:620359/150‑1 (MQ=255)
                 gAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttt                                                                                                                                      >  1:695323/1‑150 (MQ=255)
                 caaCGTTGAATCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCTTGCAGATTACTGAGGATTTTCAATCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTTTAATGAATACAATGTGCttt                                                                                                                                      <  1:270443/149‑1 (MQ=255)
                    cGTTGAATCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATTAAAATCACATTGTTTAATTAATACAATGTGCttttta                                                                                                                                   <  2:452078/150‑2 (MQ=255)
                    cGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCtttttt                                                                                                                                   <  1:1597630/150‑1 (MQ=255)
                     gTTGAAGCTGCCATAACACACTCCTGAAATTGTAAAAAATGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttttttt                                                                                                                                  <  2:1414130/150‑1 (MQ=255)
                     gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttttttt                                                                                                                                  >  1:448347/1‑150 (MQ=255)
                     gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttttttt                                                                                                                                  <  1:1130128/150‑1 (MQ=255)
                     gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttttttt                                                                                                                                  <  1:388416/150‑1 (MQ=255)
                        gAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTaaa                                                                                                                               <  2:656623/150‑1 (MQ=255)
                         aaGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAAt                                                                                                                              <  2:2018224/150‑1 (MQ=255)
                          aGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAAtt                                                                                                                             >  1:1678589/1‑150 (MQ=255)
                           gCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTc                                                                                                                            <  1:1515277/150‑1 (MQ=255)
                               ccATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATaa                                                                                                                        >  1:178815/1‑150 (MQ=255)
                                  tAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGaa                                                                                                                     >  2:232656/1‑150 (MQ=255)
                                    acacacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAAtt                                                                                                                   >  1:792344/1‑150 (MQ=255)
                                       cacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTg                                                                                                                <  2:1346914/150‑1 (MQ=255)
                                       aacTCCTGAAAGTGTAAAAAAGTGCGGTCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTg                                                                                                                <  1:984489/149‑1 (MQ=255)
                                           ccTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTAtt                                                                                                            >  1:872538/1‑150 (MQ=255)
                                              taaaGTTTAAAAAAGGGCGGGAAGCTTACTGATGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAat                                                                                                         <  1:202271/149‑1 (MQ=255)
                                              taaaGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTAGAAAATATGATCATTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAat                                                                                                         <  2:1617781/149‑1 (MQ=255)
                                              gAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTTAACACTGTGCTTTATTTAAATTCATAAGAATTTTGTATTTat                                                                                                         >  1:360023/1‑150 (MQ=255)
                                               aaaTTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAatt                                                                                                        <  1:1262765/150‑1 (MQ=255)
                                                  gtgtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAatttat                                                                                                     >  1:1065914/1‑150 (MQ=255)
                                                  gtgtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAatttat                                                                                                     >  1:510430/1‑150 (MQ=255)
                                                       aaaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaa                                                                                                >  1:1724338/1‑150 (MQ=255)
                                                        aaaaGGGCTGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaaa                                                                                               <  2:235612/150‑1 (MQ=255)
                                                        aaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaaa                                                                                               <  1:713788/150‑1 (MQ=255)
                                                          aaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAAc                                                                                             <  1:1817387/150‑1 (MQ=255)
                                                          aaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAAc                                                                                             <  2:447468/150‑1 (MQ=255)
                                                                         aCTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTc                                                                              <  1:1967239/150‑1 (MQ=255)
                                                                          cTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATAAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTcc                                                                             >  2:1823857/1‑150 (MQ=255)
                                                                           tGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCAGATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCt                                                                            <  2:480619/150‑1 (MQ=255)
                                                                            gAGGATTTTGATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACTATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCtt                                                                           <  2:794485/150‑1 (MQ=255)
                                                                              ggATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAt                                                                         <  1:1257788/150‑1 (MQ=255)
                                                                               taTTTTAATTAGAAAATATTAACATTATTATTGATTATCAAAATTAATCTAATAAAAATAACATTGTATTAATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTTTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAtt                                                                        <  1:1810344/149‑1 (MQ=255)
                                                                               gATTTTCATTAGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAAAAAAAATCACATTGTATTCATTAAACAATGTGCTTTTTTTATATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAtt                                                                        <  1:1818516/150‑1 (MQ=255)
                                                                                 ttttCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTTTTAATTTATTTTAAAACTAAATTAATCATTTCATTAtttt                                                                      >  2:317161/1‑150 (MQ=255)
                                                                                  tttCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAttttt                                                                     >  2:912365/1‑150 (MQ=255)
                                                                                    tCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTcc                                                                   <  2:723277/150‑1 (MQ=255)
                                                                                    tCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTcc                                                                   <  2:158979/150‑1 (MQ=255)
                                                                                      aTTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                                                                 >  1:693103/1‑150 (MQ=255)
                                                                                       ttCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                                                                 >  1:2015358/1‑150 (MQ=255)
                                                                                          gAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaat                                                             >  1:788971/1‑150 (MQ=255)
                                                                                                  gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGATTTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaataataattc                                                     >  1:400573/1‑148 (MQ=255)
                                                                                                  gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATga                                                     >  1:427506/1‑149 (MQ=255)
                                                                                                  gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGc                                                     >  2:1295223/1‑150 (MQ=255)
                                                                                                  gATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGc                                                     >  2:1861844/1‑150 (MQ=255)
                                                                                                  gATCAGGATTATTGATAATCAAAATTAATCTAATAAAAAGCACAGGGGATTCATTAAACAATGTGCTTTTTTTAAATTCATATGAATTTTGTATTAATTTATTTTAAATCTGAATTAAGCATTTCCTTATTTTTCCGGATATAATAATGc                                                     >  2:295617/1‑150 (MQ=255)
                                                                                                           tattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCaa                                            <  1:911004/150‑1 (MQ=255)
                                                                                                           tattGATTATCAAAAATAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCaa                                            >  2:198372/1‑150 (MQ=255)
                                                                                                            attGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAAt                                           >  2:1337515/1‑150 (MQ=255)
                                                                                                              tGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATaa                                         >  2:1440723/1‑150 (MQ=255)
                                                                                                                aTTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAt                                       <  2:1264862/150‑1 (MQ=255)
                                                                                                                   aTCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGAt                                    >  1:1983920/1‑150 (MQ=255)
                                                                                                                   aTCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTATCCGGAAATAATAATGCGTCGCGCAATAAATGAt                                    <  2:1502199/150‑1 (MQ=255)
                                                                                                                    taaaaaTTAATATAATAAAATGAAAATTTTATTAATTAAAAAATTTGATTTTTTAAAATTCATAACAATTTTGTATTAATTAATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATc                                   <  1:166164/148‑1 (MQ=255)
                                                                                                                    tCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATc                                   <  2:883281/150‑1 (MQ=255)
                                                                                                                    tCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATc                                   >  2:139139/1‑150 (MQ=255)
                                                                                                                    tCAAAATTAATCTAATAAAAAGCAAATTGTTTTCATTAAACAATTTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCTGAAATAATAATGCGTCGCGCAATAAATGATc                                   <  2:446637/150‑1 (MQ=255)
                                                                                                                     cAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCa                                  >  1:1196740/1‑150 (MQ=255)
                                                                                                                     cAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAAATTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCa                                  >  2:508843/1‑150 (MQ=255)
                                                                                                                      aaaaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCAt                                 >  1:523594/1‑150 (MQ=255)
                                                                                                                        aaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTCCTGAAATAATAAAGAGTCGATCAATAAATGATCAttt                               >  2:360286/1‑150 (MQ=255)
                                                                                                                               cTAATTAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                        <  1:1021476/150‑1 (MQ=255)
                                                                                                                               cTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                        >  2:1964015/1‑150 (MQ=255)
                                                                                                                               cTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                        >  2:347566/1‑150 (MQ=255)
                                                                                                                               cTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                        >  2:4119/1‑150 (MQ=255)
                                                                                                                               cTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATAAATCATTTCCTTATTTTTCCGGAAAAAATAATGCTTCGCGCAATAAATGATCATTTCCacgac                        >  1:1794713/1‑150 (MQ=255)
                                                                                                                                taataaAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgacg                       >  2:870315/1‑150 (MQ=255)
                                                                                                                                taataaAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCAGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgacc                       >  2:1115799/1‑149 (MQ=255)
                                                                                                                                 taaaaaaaGCACATTGTTTAATGAATACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCTCAATAAATGATCATTTCCacgacga                      >  2:1726954/4‑150 (MQ=38)
                                                                                                                                 aataaAAAGCAAAATGTATTCATTAAACAATGAGCTTTTAATAAATTAATAAGAATATTGTATTAATTTATATTAAAACTTAATTAATCATTTCCTTATTTTTCCGGATAAAATAATGCGTCGCGCAATAAATGATCATTACCacgacga                      >  1:1670670/1‑150 (MQ=255)
                                                                                                                                  ataaAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgacgac                     <  1:402945/150‑1 (MQ=255)
                                                                                                                                    aaaaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATTCGTCGCGCAATAAATGATCATTTCCACGACGAcgc                   >  2:1013489/1‑150 (MQ=255)
                                                                                                                                    aaaaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGAcgc                   >  2:744625/1‑150 (MQ=255)
                                                                                                                                     aaaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGAcgcg                  >  2:1658709/1‑150 (MQ=255)
                                                                                                                                      aaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCTGAAATAATAATGCGTCGCGCAATAAATCATCATATCCACGACTAcgctt                 >  1:1235687/1‑148 (MQ=255)
                                                                                                                                      aaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGt                 >  2:1799288/1‑150 (MQ=255)
                                                                                                                                      aaaGCACATTGTTTAATGAATACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGt                 >  1:1242682/1‑150 (MQ=255)
                                                                                                                                       aaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTa                <  2:510430/150‑1 (MQ=255)
                                                                                                                                         gCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACTACGCGTaaa              >  2:190331/1‑150 (MQ=255)
                                                                                                                                          cacaTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAg             >  1:1189075/1‑150 (MQ=255)
                                                                                                                                                ttATTCATTATACAATGTGATTTTTTTAAATTCTTAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATGATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATg       <  2:1126598/149‑1 (MQ=255)
                                                                                                                                                gTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATg       >  2:1820424/1‑150 (MQ=255)
                                                                                                                                                 tATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGc      <  1:561117/150‑1 (MQ=255)
                                                                                                                                                    tCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGt   <  2:1422767/150‑1 (MQ=255)
                                                                                                                                                     cATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTc  <  2:178815/150‑1 (MQ=255)
                                                                                                                                                     cATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAAAAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGAggaa  >  2:822346/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTC  >  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: