Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,968 G→T 6.0% intergenic (+366/+248) 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,9680GT 6.0% 623.6 / 21.5 192intergenic (+366/+248)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (54/126);  new base T (5/7);  total (59/133)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 192 total aligned reads displayed.

ATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAA  >  NC_000913/3269820‑3270110
                                                                                                                                                    |                                                                                                                                              
 tGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGaa                                                                                                                                              >  1:2165568/1‑150 (MQ=255)
     tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                          <  2:1398767/150‑1 (MQ=255)
     tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                          <  1:1029638/150‑1 (MQ=255)
     tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                          <  2:2347628/150‑1 (MQ=255)
     tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCAGCCGAAGCgg                                                                                                                                          >  1:3999106/1‑150 (MQ=255)
     tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCGGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                          <  1:1548941/150‑1 (MQ=255)
     tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCGGTCTTTGGGGTATAGTGTCGTGGACCGTCATTCATCGTTCTGCCCCTCCAACAGGAAAAACCCGCCGAAGCgg                                                                                                                                          >  2:2411311/1‑150 (MQ=255)
        aGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtt                                                                                                                                       <  1:2594785/150‑1 (MQ=255)
        aGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtt                                                                                                                                       <  1:1416597/150‑1 (MQ=255)
        aGCGCAATCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtt                                                                                                                                       <  1:3152395/150‑1 (MQ=255)
         gcgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGttt                                                                                                                                      <  2:1220905/150‑1 (MQ=255)
          cgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGtttt                                                                                                                                     >  1:3261021/1‑150 (MQ=255)
           gcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                                    <  1:200872/150‑1 (MQ=255)
            cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTa                                                                                                                                   <  2:1405140/150‑1 (MQ=255)
             aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTAc                                                                                                                                  <  1:248942/150‑1 (MQ=255)
             aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAACCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                                  <  1:4221711/150‑1 (MQ=255)
               tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                <  1:659637/150‑1 (MQ=255)
               tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                >  1:206323/1‑150 (MQ=255)
               tCAACAATCTGAAACCCACCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                >  2:1253160/1‑150 (MQ=255)
                   caaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                            <  1:2276572/150‑1 (MQ=255)
                   caaGCTTAAACCCGCCGGAGCGGTTTTTTACGTACATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaac                                                                                                                            <  2:2168015/150‑2 (MQ=11)
                    aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAc                                                                                                                           >  2:720215/1‑150 (MQ=255)
                      tCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                                                                                         <  1:3844684/150‑1 (MQ=255)
                         gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                                      <  2:4290787/150‑1 (MQ=255)
                         gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                                      <  1:1890223/150‑1 (MQ=255)
                          taaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGTCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTGATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                                                                                                                     <  2:3831928/149‑1 (MQ=255)
                            aCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                                   <  2:3064828/150‑1 (MQ=255)
                              ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                                                                                                 >  2:606671/1‑150 (MQ=255)
                               cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAGTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTg                                                                                                                >  1:214047/1‑150 (MQ=255)
                                  cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                             <  1:2192136/150‑1 (MQ=255)
                                       cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                        <  2:636799/150‑1 (MQ=255)
                                         gTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                                      <  1:2903907/150‑1 (MQ=255)
                                             attaCGTACATCAGGTAAAACTGACCGATAAGACGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTATGCCCCTCCCCAAGAAAAAACCAGCCGCACCGGGTGTTTCAGTAAAACAGGTGAAACTCACCGCTAAGCCggg                                                                                                  >  2:2233679/1‑148 (MQ=14)
                                                 cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                              <  2:396643/150‑1 (MQ=255)
                                                 cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                              <  2:2596720/150‑1 (MQ=255)
                                                 cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                              <  2:868128/150‑1 (MQ=255)
                                                  gTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTGTTTACGTAAAACAGGTGAAACTGAACGATAAGCCGcgttctt                                                                                             >  1:3854449/1‑150 (MQ=255)
                                                  gTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACACGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctg                                                                                             >  1:2049362/1‑149 (MQ=255)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                           <  2:286916/150‑1 (MQ=40)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                           <  1:2605968/150‑1 (MQ=40)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                           <  1:1971126/150‑1 (MQ=40)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                           <  2:1373531/150‑1 (MQ=255)
                                                     cATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAATCGTGTTTTTACGTAAAACCGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                          <  2:1522275/150‑1 (MQ=255)
                                                     cATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGTAAAAAACCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACAGACCCATAAGCCGCGTTCTTTTg                                                                                          >  1:811470/1‑150 (MQ=14)
                                                     cATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAACTTTCTGCCCCTCCAAAACTAAAAACCCGCCGAACCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                          <  1:526207/150‑1 (MQ=255)
                                                       tCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGGAAAACAGGTGGAACTGACCGATAAGCCGCTGTCGGTTggg                                                                                        >  1:1296964/1‑150 (MQ=11)
                                                         aGGTAATACTGTCCGATAAGCCGCTTTCTTTTGGGTTTAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGta                                                                                      <  2:2094026/150‑1 (MQ=255)
                                                               aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt                                                                                <  1:1588847/150‑1 (MQ=37)
                                                                     ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGa                                                                          <  1:1444125/150‑1 (MQ=37)
                                                                       gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                                        <  1:722198/150‑1 (MQ=37)
                                                                          caggagCTTTCTGTTTGGTATAGAGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTCAAACAGGTGAAACTGTCCGATTAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                                     <  1:4164472/145‑1 (MQ=11)
                                                                          aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                                     <  1:2015683/150‑1 (MQ=37)
                                                                             ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcatt                                                                  <  1:3059617/150‑1 (MQ=37)
                                                                             ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGTAAAAACCCGCCGAAGCCGGTTTTTACGTAAAACAGGTGAAACAGACCGATAAGCCCCTTTATTTTGGGTATAGTGTCGTGGACACtcaat                                                                  >  2:3670407/1‑148 (MQ=11)
                                                                                        tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                       <  1:4331418/150‑1 (MQ=37)
                                                                                           gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcccc                                                    >  2:4234811/1‑150 (MQ=37)
                                                                                           gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcccc                                                    <  1:3749356/150‑1 (MQ=37)
                                                                                            tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCt                                                   >  2:4328113/1‑150 (MQ=37)
                                                                                             ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                                                  <  2:1478101/150‑1 (MQ=37)
                                                                                             ataGTGTCGTGGACAGTCATTCATCTGTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                                                  <  1:1846327/150‑1 (MQ=37)
                                                                                              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                                             <  1:1031296/150‑1 (MQ=37)
                                                                                                  gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                             <  2:4096402/150‑1 (MQ=37)
                                                                                                  gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                             <  1:2547047/150‑1 (MQ=37)
                                                                                                      tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaa                                         <  1:2389424/150‑1 (MQ=37)
                                                                                                      tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaa                                         <  1:372498/150‑1 (MQ=37)
                                                                                                          cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAc                                     <  2:2787252/150‑1 (MQ=37)
                                                                                                          cAGTCATTCATCTTTCTGCCCCTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAc                                     <  2:2787254/150‑1 (MQ=37)
                                                                                                             tcattcatCTTTCTGCCCCTCAAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                  <  1:3893728/150‑1 (MQ=255)
                                                                                                              cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc                                 <  2:1936327/150‑1 (MQ=37)
                                                                                                              cattcatCTTTCTGCCCCTCAAAAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc                                 <  1:2733771/150‑1 (MQ=255)
                                                                                                               attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc                                <  2:2532334/150‑1 (MQ=37)
                                                                                                               attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAAAGTGTCGTGGACAGTCAATAATCATTCAGCCCCTCCAAAAGCAAAAAcccccc                                >  1:3894987/1‑150 (MQ=11)
                                                                                                                  catCTTTCTGCCCCTCCAGAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgccgcc                             >  2:1340983/1‑148 (MQ=19)
                                                                                                                  catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTTCGTAAAACTGGTGAAACTGACCGATATGCCGCTTTCTTTGGGGTATAGTGTCGTGGACAGTCATTCTTCGTTCTGCCCCTCCAAAATCAAAAACCCGCCGaa                             >  1:1687632/1‑150 (MQ=11)
                                                                                                                  catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGACACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                             >  1:1187333/1‑150 (MQ=255)
                                                                                                                   atatttCTGCCCCTCAAAAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                            <  2:3999106/147‑1 (MQ=255)
                                                                                                                   atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTTCCCCTCCAAAAGCAAAAACCCGCCGAAg                            >  2:1572639/1‑150 (MQ=37)
                                                                                                                   atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACACTCAATCATCTTTCTGCCCCTACAAAAGCAAAAACCCGCCGAAg                            >  2:1158106/1‑150 (MQ=255)
                                                                                                                      tttCTGCCCCTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                         <  2:3327757/150‑1 (MQ=37)
                                                                                                                      tttCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCgg                         <  2:2696876/150‑1 (MQ=255)
                                                                                                                      tttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                         <  1:719767/150‑1 (MQ=37)
                                                                                                                            ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggggtgtt                   >  2:2395629/1‑145 (MQ=255)
                                                                                                                             cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                  <  1:2849654/150‑1 (MQ=37)
                                                                                                                             cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGGCATTCATCGTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTa                  >  2:1795303/1‑150 (MQ=11)
                                                                                                                              ccTCAAAAATCAAAAACCCGCCGTAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                 <  2:2727999/150‑1 (MQ=11)
                                                                                                                               cTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                >  2:1838103/1‑150 (MQ=37)
                                                                                                                               cTCCAAAAGTAAAAACCCGCCGAAGCGGGGTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATCGTGTCGTGGCCAGTCATTCATCGTGCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGGTGTTACg                >  2:3391048/1‑150 (MQ=11)
                                                                                                                                   acaaGTAAAAACCCGCCGCATCGGGTTTTTACGTACATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa            <  2:2133568/148‑1 (MQ=255)
                                                                                                                                   aaaaGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa            <  1:3836250/150‑1 (MQ=37)
                                                                                                                                      aGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa         <  1:3355116/150‑1 (MQ=255)
                                                                                                                                       gTAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg        >  2:530509/1‑150 (MQ=11)
                                                                                                                                       gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg        <  2:4088576/150‑1 (MQ=25)
                                                                                                                                        tAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg       >  2:965501/1‑150 (MQ=11)
                                                                                                                                        tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg       >  1:1618051/1‑150 (MQ=17)
                                                                                                                                          aaaaCCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCATTCTGACCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAACCAGGTg     >  1:2430227/1‑150 (MQ=11)
                                                                                                                                                    |                                                                                                                                              
ATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAA  >  NC_000913/3269820‑3270110

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: