Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,970 A→G 6.3% intergenic (+368/+246) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,9700AG 6.3% 572.4 / 21.5 189intergenic (+368/+246)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (52/125);  new base G (5/7);  total (57/132)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 192 total aligned reads displayed.

TGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAT  >  NC_000913/3269821‑3270119
                                                                                                                                                     |                                                                                                                                                     
    tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                                   <  1:1029638/150‑1 (MQ=255)
    tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                                   <  2:2347628/150‑1 (MQ=255)
    tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                                   <  2:1398762/150‑1 (MQ=255)
    tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                                   <  1:2590592/150‑1 (MQ=255)
    tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAACTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                                   >  1:4061294/1‑150 (MQ=255)
    tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCGGTCTTTGGGGTATAGTGTCGTGGACCGTCATTCATCGTTCTGCCCCTCCAACAGGAAAAACCCGCCGAAGCgg                                                                                                                                                   >  2:2411311/1‑150 (MQ=255)
       aGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtt                                                                                                                                                <  1:1416597/150‑1 (MQ=255)
        gcgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGttt                                                                                                                                               <  1:728776/150‑1 (MQ=255)
        gcgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGttt                                                                                                                                               <  2:1220905/150‑1 (MQ=255)
        gcgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                                                                                                                               >  2:1588847/1‑150 (MQ=255)
        gcgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGGGTCCGGTCATACAGCTTTCTGCCCCGACATAAGTAAATACCCGCCGAATCGGGttt                                                                                                                                               >  2:1417421/1‑150 (MQ=18)
          gcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                                             <  1:200872/150‑1 (MQ=255)
          gcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTAAATCAGTTAAAACTGACCTATAAGCCGCTTTATTTTGGGTATAGTGTCTTGGACAGTCTTTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGtttttt                                                                                                                                             <  1:4250514/150‑1 (MQ=17)
           cACTCAACACTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTa                                                                                                                                            <  2:1405112/150‑1 (MQ=18)
           cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                                            <  1:592146/150‑1 (MQ=255)
            aCTCAACAATCTGAAACCCGCCGGACCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAACCGCTTTCTTTTGGTTATAGTGTCGTGGACTGTCATTCATCTTTCTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                                           <  1:4221664/150‑1 (MQ=17)
            aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTAc                                                                                                                                           <  1:248958/150‑1 (MQ=255)
              tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                         >  1:2168015/1‑150 (MQ=255)
              tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                         >  1:2667613/1‑150 (MQ=255)
              tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTGTTTACgg                                                                                                                                         >  1:2851158/1‑149 (MQ=255)
              tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAAGTCAATCATCTTTCAGCCCCTACCAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                         >  2:1253179/1‑150 (MQ=18)
                aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                                                       >  1:3965946/1‑150 (MQ=255)
                  caaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaat                                                                                                                                     <  2:2310593/150‑2 (MQ=255)
                   aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAc                                                                                                                                    >  2:720215/1‑150 (MQ=255)
                        gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                                               <  2:4290787/150‑1 (MQ=255)
                        gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                                               <  1:2401991/150‑1 (MQ=255)
                         taaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGTCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTGATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                                                                                                                              <  2:3831928/149‑1 (MQ=255)
                           aCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                                            <  1:3754993/150‑1 (MQ=255)
                           aCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                                            <  2:3064828/150‑1 (MQ=255)
                             ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                                                                                                          >  2:606671/1‑150 (MQ=255)
                              cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTg                                                                                                                         >  1:214132/1‑150 (MQ=255)
                                ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                                       <  1:236060/150‑1 (MQ=255)
                                 cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                                      <  1:3142614/150‑1 (MQ=255)
                                      cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                                 <  2:636799/150‑1 (MQ=255)
                                                cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                                       <  2:2596720/150‑1 (MQ=255)
                                                cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                                       <  2:868128/150‑1 (MQ=255)
                                                cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                                       <  2:1527817/150‑1 (MQ=255)
                                                cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                                       <  2:396643/150‑1 (MQ=255)
                                                 gTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTGATTGATATTTGTGGCCCTCCAAAAGGTAAAAACCGCCGGAGGGGGTTTTTTGGTTAATCAGGTGTTAATTTACGATAAGCCgggttgtt                                                                                                      >  1:1950400/1‑143 (MQ=2)
                                                 gTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTGTTTACGTAAAACAGGTGAAACTGAACGATAAGCCGcgttctt                                                                                                      >  1:3854449/1‑150 (MQ=255)
                                                 gTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGggttctg                                                                                                      >  1:2049522/1‑149 (MQ=255)
                                                   aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                    <  2:286916/150‑1 (MQ=40)
                                                   aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                    <  1:2605968/150‑1 (MQ=40)
                                                    cATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                                   <  1:2244780/150‑1 (MQ=39)
                                                      tCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGGAAAACAGGTGGAACTGACCGATAAGCCGCTGTCGGTTggg                                                                                                 >  1:1296964/1‑150 (MQ=11)
                                                        aGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGta                                                                                               <  2:1238034/150‑1 (MQ=37)
                                                          gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                                                             <  2:1924730/150‑1 (MQ=37)
                                                                    ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGa                                                                                   <  2:2118480/150‑1 (MQ=37)
                                                                    ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGa                                                                                   <  1:1444125/150‑1 (MQ=37)
                                                                      gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                                                 <  1:715788/150‑1 (MQ=37)
                                                                         aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                                              <  1:2015655/150‑1 (MQ=37)
                                                                           gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcat                                                                            <  1:2433650/150‑1 (MQ=37)
                                                                            ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcatt                                                                           <  1:1617224/150‑1 (MQ=37)
                                                                            ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcatt                                                                           <  1:3888075/150‑1 (MQ=37)
                                                                            ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGTAAAAACCCGCCGAAGCCGGTTTTTACGTAAAACAGGTGAAACAGACCGATAAGCCCCTTTATTTTGGGTATAGTGTCGTGGACACtcaat                                                                           >  2:3670407/1‑148 (MQ=11)
                                                                            ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTCAGCCCCTCCAAAAGTAAAAACCCCCCGAAGCCGGTTTTTACGTAAAAAAGGTGAAACAGACCGATAATCCGCTTTCTTTTGGGTATAGAGTCGTGGACACtcctt                                                                           >  2:1508692/1‑150 (MQ=11)
                                                                                    ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTc                                                                   <  2:2165568/150‑1 (MQ=37)
                                                                                        gggTATAGTGTCGTGGACAGTCTTTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcc                                                               <  1:3505953/150‑1 (MQ=37)
                                                                                          gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcccc                                                             <  1:3749173/150‑1 (MQ=37)
                                                                                          gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcccc                                                             >  2:4234811/1‑150 (MQ=37)
                                                                                           tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCt                                                            >  2:4328113/1‑150 (MQ=37)
                                                                                            ttaGTGTCGTGGACAGTCATTCAACTTTCTGCCCCTCCAAAAGTACAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                                                           <  1:1846607/149‑1 (MQ=255)
                                                                                             taGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAAATTTCTGCCCCTcc                                                          >  2:914177/1‑150 (MQ=255)
                                                                                                tgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaa                                                       <  2:3438646/150‑1 (MQ=37)
                                                                                                tgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaa                                                       >  2:516609/1‑150 (MQ=17)
                                                                                                ggtCGTGGTCATTCTTTCATCTTTCTGCCCCTCCAAAAGTTTAATGCCGGCGAAGCGGGTTTTTACGTAATACTGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaa                                                       <  2:3167938/149‑1 (MQ=11)
                                                                                                 gtCGTTGACAGTCATTCATCTTTCTGCCCCTACAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGTTGAAACTTACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGTGGACAGTCATTCATCTTTCTGACCCTCCaaaa                                                      <  2:2000868/150‑1 (MQ=11)
                                                                                                 gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                      <  2:4096402/150‑1 (MQ=37)
                                                                                                 gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                      <  1:1031296/150‑1 (MQ=37)
                                                                                                  tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                     >  1:277085/1‑150 (MQ=37)
                                                                                                      ggACAGTCAGTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaa                                                 <  2:2220475/150‑1 (MQ=37)
                                                                                                        aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaaa                                               <  1:3008653/150‑1 (MQ=37)
                                                                                                         cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAc                                              <  2:2787252/150‑1 (MQ=37)
                                                                                                         cAGTCATTCATCTTTCTGCCCCTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAc                                              <  2:2787254/150‑1 (MQ=37)
                                                                                                           ctcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAccc                                            <  1:3246098/149‑1 (MQ=37)
                                                                                                            tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                           <  1:3893732/150‑1 (MQ=37)
                                                                                                              attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc                                         <  2:2532334/150‑1 (MQ=37)
                                                                                                                tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATTGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                       <  1:3322029/150‑1 (MQ=37)
                                                                                                                 catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGACACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                      >  1:1187333/1‑150 (MQ=255)
                                                                                                                 catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCCTTCATCTGTCGGCCCCTCCAAAAGCAAAAACccgccgca                                      >  1:92926/1‑148 (MQ=255)
                                                                                                                  atatttCTGCCCCTCAAAAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                     <  2:3999106/147‑1 (MQ=255)
                                                                                                                     tttCTGCCCCTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                  <  2:3327757/150‑1 (MQ=37)
                                                                                                                     tttCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCgg                                  <  2:2696876/150‑1 (MQ=255)
                                                                                                                     tttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACATGTGAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                  <  1:720215/150‑1 (MQ=255)
                                                                                                                          gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                             <  2:1803841/150‑1 (MQ=37)
                                                                                                                          gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                             >  2:2028207/1‑150 (MQ=37)
                                                                                                                           ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggggtgtt                            >  2:2395629/1‑145 (MQ=255)
                                                                                                                             ccTCAAAAATCAAAAACCCGCCGTAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                          <  2:2727999/150‑1 (MQ=11)
                                                                                                                              cTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCCAAAGCAAAAACCCGCCGAAGCGGGGTTTTACg                         >  2:1355318/1‑150 (MQ=255)
                                                                                                                              cTCCAAAAGTAAAAACCCGCCGAAGCGGGGTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATCGTGTCGTGGCCAGTCATTCATCGTGCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGGTGTTACg                         >  2:3391048/1‑150 (MQ=11)
                                                                                                                                  aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                     <  1:3836251/150‑1 (MQ=37)
                                                                                                                                     aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa                  <  1:2086136/150‑1 (MQ=39)
                                                                                                                                      gTAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                 >  2:530509/1‑150 (MQ=11)
                                                                                                                                      gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                 <  2:4088576/150‑1 (MQ=25)
                                                                                                                                       tAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                >  2:965501/1‑150 (MQ=11)
                                                                                                                                       tAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGGCCTTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCCCCGCAGCGGGTTTTTTCGGAAAACAgg                >  1:433480/1‑150 (MQ=11)
                                                                                                                                        aaaaaCCCGCCGAAGCGGGTTTTTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt               <  1:1570119/150‑1 (MQ=37)
                                                                                                                                         aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTGTCTTTTGGGTCTAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTGTTTACGTACACCAGGTg              >  1:1670486/1‑150 (MQ=11)
                                                                                                                                          aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGa             <  1:2031618/150‑1 (MQ=37)
                                                                                                                                          aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGa             <  1:180063/150‑1 (MQ=37)
                                                                                                                                                     |                                                                                                                                                     
TGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAT  >  NC_000913/3269821‑3270119

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: