Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,969 A→C 5.4% intergenic (+367/+247) 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,9690AC 5.4% 717.0 / 21.2 244intergenic (+367/+247)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (75/151);  new base C (6/9);  total (82/162)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

ATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCG  >  NC_000913/3269820‑3270117
                                                                                                                                                     |                                                                                                                                                    
aTGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                                                                                                                                      >  2:4109001/1‑150 (MQ=255)
aTGAAAAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTAGGGTATAGCGTCGTGGACAGTCATTCATAATTCAGCCCCACCAAAAGCAAAAACCCGCCGa                                                                                                                                                      >  2:1053923/1‑150 (MQ=18)
    aTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCg                                                                                                                                                  >  1:2145468/1‑150 (MQ=255)
     tAAAGCGCACTCAACAATCTTAAACCCGCCGGATCGGTTTATTACGTACTTCAGGTAAAACTGACCGATACGCCGATTTTTTTTGGGTATAGTGTCGTTGACAGTCTTTCTTCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                                 <  1:967559/150‑1 (MQ=18)
     tAAAGCGCACTCAACAATCTGACACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                                 <  2:3890247/150‑1 (MQ=255)
     tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCCATCATCTTTATGCCCCACCAAAAGAAAAAACCCGCCGAAGCgg                                                                                                                                                 >  1:4219225/1‑150 (MQ=18)
     tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                                 <  1:4085191/150‑1 (MQ=255)
       aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGACCCTCCAAAAGTAAAAACCCGCCGAAGCGGGt                                                                                                                                               >  1:3480646/1‑150 (MQ=255)
        agggcACTCCACAATCTGAAACCCGCCGTACCGGTTTATTACGTACATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                                                                                                              <  2:1681952/147‑1 (MQ=14)
         gcgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGttt                                                                                                                                             <  1:2755967/150‑1 (MQ=255)
           gcACTCAACAATCAGAAACCCGCCGTATCGGTTTATTACGTACAGCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                                                                                                           <  2:3480646/150‑1 (MQ=14)
            cACTCCACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTa                                                                                                                                          <  2:2710001/150‑1 (MQ=18)
              cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                                        <  1:1654099/150‑1 (MQ=255)
              cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                                        <  1:2988805/150‑1 (MQ=255)
               tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                       <  2:2981999/150‑1 (MQ=255)
               tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                       <  2:1497243/150‑1 (MQ=255)
               tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                       >  1:997156/1‑150 (MQ=255)
               tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                       >  2:152922/1‑150 (MQ=255)
               tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                                       >  1:352822/1‑150 (MQ=255)
                caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTActtt                                                                                                                                      <  2:3381301/150‑4 (MQ=17)
                  acaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                                                                                                                                            >  1:2427398/1‑142 (MQ=255)
                   caaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                                   >  1:1032651/1‑150 (MQ=255)
                       cTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                                                                                                               <  1:3978276/150‑1 (MQ=255)
                       cTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                                        <  2:3218373/141‑1 (MQ=255)
                        tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                              <  1:3619142/150‑1 (MQ=255)
                        tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                              <  1:1585470/150‑1 (MQ=255)
                         gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                                             <  1:4138593/150‑1 (MQ=255)
                              ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                                                                                                        <  1:145128/150‑1 (MQ=255)
                              ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTGTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                                                                                                        >  1:4018775/1‑150 (MQ=255)
                               cgccgGAACGGTTTATTACGTACATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGACGTGGACAGTCATTCAACTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGCAGCGGGGTTTTACGTAAAACAGGTGAAAATg                                                                                                                       >  2:868325/1‑150 (MQ=14)
                                 ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                                     <  1:574485/150‑1 (MQ=255)
                                  cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                                    <  1:3383325/150‑1 (MQ=255)
                                     aaCGGTTTATTACGTACATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGAt                                                                                                                 <  1:597823/150‑1 (MQ=9)
                                     aaCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGAt                                                                                                                 <  2:3979782/150‑1 (MQ=255)
                                     aaCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                                 <  1:2205462/150‑1 (MQ=255)
                                       cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATCGTTCAGCCCCACCAAAAGTAAAAACCCCCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACAGACCGATaa                                                                                                               >  2:3937226/1‑150 (MQ=14)
                                         gTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                                             >  1:279301/1‑150 (MQ=255)
                                           ttattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAACTCAATAGTGTGTATGCCCCTAAAAAAGTAAAAAACCGCCTAATAGGGTTTTTACGGAAAACAGGTGAAACAGACCGATAAGGCg                                                                                                           >  2:162726/1‑150 (MQ=14)
                                                  gTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCAATCAACTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGGGTTTTAAGTAAAAAAGGTGAAACTGACCGATAATCCGcattctt                                                                                                    >  2:1987364/1‑150 (MQ=14)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGggttctgtc                                                                                                  <  2:2225368/150‑4 (MQ=255)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                  <  2:2569104/150‑1 (MQ=40)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                  <  2:2197577/150‑1 (MQ=40)
                                                     cATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                                 <  1:1044665/150‑1 (MQ=39)
                                                      aTCAGGTAAAACTGACCGATAAGCCTCTTTCTTTTTGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                <  1:804167/150‑1 (MQ=255)
                                                      aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                <  1:3781863/150‑1 (MQ=38)
                                                      aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                <  2:3956723/150‑1 (MQ=38)
                                                      aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                <  2:2284591/150‑1 (MQ=38)
                                                      aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                <  1:2216215/150‑1 (MQ=38)
                                                       tCAGGTAAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAACATTGTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTTAGTAAAACAGGTGGAACTGACCGATAAGGCGCTGTCGTTTggg                                                                                               >  1:1293547/1‑150 (MQ=11)
                                                           gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                                                           <  1:1133557/150‑1 (MQ=37)
                                                           gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                                                           <  1:2923398/150‑1 (MQ=37)
                                                           gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATATTAATTCCCCTCCAAAATTAAAAAAACGCCGAAGCCGGTTTTTTAGTAAAACAGGTGTAAAGTACAGAAAAGCCGCTTTCGTGTGGGtaat                                                                                           >  1:3224922/1‑148 (MQ=2)
                                                                aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTc                                                                                      <  1:1937394/150‑1 (MQ=37)
                                                                  tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                                    <  1:3126896/150‑1 (MQ=37)
                                                                     ccGATCAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTCTGCCCCTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGa                                                                                 <  2:2442328/150‑1 (MQ=255)
                                                                     ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGa                                                                                 <  1:1679507/150‑1 (MQ=37)
                                                                     ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGa                                                                                 <  2:3409018/150‑1 (MQ=37)
                                                                          aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                                            <  2:1856895/150‑1 (MQ=37)
                                                                             ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCACCCCCACCAAAAAAAAAAACAAGCCGAAGACCCTTTTAAAGAAAAACAGGTGAAACTGACCGAAAAGCCGCATTCTTCTGGTTAAAGTCTCGTGGAcgcgcatg                                                                         >  2:3816900/1‑143 (MQ=2)
                                                                                ctttctttTGGTTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattcat                                                                      <  1:2948417/150‑1 (MQ=37)
                                                                                 tttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATc                                                                     <  1:180775/150‑1 (MQ=37)
                                                                                      tttGGGTATCGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCt                                                                <  2:4207871/150‑1 (MQ=37)
                                                                                         gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcc                                                             <  2:2173141/150‑1 (MQ=37)
                                                                                         gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcc                                                             <  2:1447270/150‑1 (MQ=37)
                                                                                         gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcc                                                             <  2:3184568/150‑1 (MQ=37)
                                                                                          ggTTTAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGccc                                                            <  2:354684/150‑1 (MQ=255)
                                                                                          ggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGccc                                                            <  2:1363679/150‑1 (MQ=37)
                                                                                            tataGTGTCGTTGACAGTCATTCATATTTCTGCCCCTCCAAAATTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCt                                                          <  2:2047780/150‑1 (MQ=255)
                                                                                            tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCt                                                          <  1:1207381/150‑1 (MQ=37)
                                                                                                gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaa                                                      <  1:4188711/150‑1 (MQ=37)
                                                                                                  gtCGTGGACAGTCTTTCATCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGCATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                    <  2:3963300/150‑1 (MQ=255)
                                                                                                  gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                    <  1:2842434/150‑1 (MQ=37)
                                                                                                  gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                    <  1:1664585/150‑1 (MQ=37)
                                                                                                   tCGTGGACAGTCATTCATCTTTCTGCCCCTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                   <  1:1599234/150‑1 (MQ=37)
                                                                                                      tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGGTTTTACGTAAAACAGGTGAAACTGACCGATACGCCGCTGTCTTTTGGGGATAGTGTCGTGGACAGTCCTTCATCTTTCTGCCCCACCAAAAGTaa                                                >  1:3851213/1‑150 (MQ=11)
                                                                                                             tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                         <  2:1151169/150‑1 (MQ=37)
                                                                                                                  catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGGAGAGTGTCGAGGCCAGGCCATCATCGGAACGCGCCGCCAATACCAAGAACccgccgcc                                    >  2:625420/1‑148 (MQ=11)
                                                                                                                     cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                 >  1:3047674/1‑150 (MQ=37)
                                                                                                                     cTTTCTGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                 <  2:1252359/150‑1 (MQ=37)
                                                                                                                      tttCTGCCCCTCCAACAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                <  1:4232313/150‑1 (MQ=21)
                                                                                                                          tGCCACTCCAAAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                            <  2:1807860/150‑1 (MQ=255)
                                                                                                                           gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                               >  2:90416/1‑66 (MQ=255)
                                                                                                                           gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                           <  1:4005751/150‑1 (MQ=37)
                                                                                                                            ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGTGTCGTGGACAGTCATTCATCTTGCTGCCCCGCCAACAGCACAAACCCGCCGAAGAggggggtt                          >  2:2393937/1‑145 (MQ=11)
                                                                                                                            ccacTCCAAAATGATAAACCCGCCGAATCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGCAAAAACCCGCTTCGGCGGGttttt                          <  2:2863538/147‑1 (MQ=11)
                                                                                                                              ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                        >  1:2707728/1‑150 (MQ=37)
                                                                                                                              ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                        >  1:1069927/1‑150 (MQ=37)
                                                                                                                                tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCAGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTGTTTACGt                      >  1:3427680/1‑150 (MQ=255)
                                                                                                                                 ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGGGTCGTGGACAGTCAATCATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCCAAGCGGGTTTTTTCgga                     >  1:3213546/1‑148 (MQ=11)
                                                                                                                                   aaaaTTATAAACCCGCCGTAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCCAAACCCCGCCGAAGCGGGTTTTTACGTaaa                   <  1:3186856/150‑1 (MQ=11)
                                                                                                                                    aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc                  <  1:307735/150‑1 (MQ=38)
                                                                                                                                       gTAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTCTAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCCAAAGCAAAAACCCCCCGAAGCGGGTTTTTACGTAAACCAg               >  2:782955/1‑150 (MQ=11)
                                                                                                                                        tAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg              >  2:2477825/1‑150 (MQ=11)
                                                                                                                                           aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa           >  1:3315648/1‑150 (MQ=21)
                                                                                                                                            aaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATAGTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGGAAATCAGGTGaa          >  1:2263643/1‑150 (MQ=255)
                                                                                                                                              cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc        >  1:192993/1‑150 (MQ=37)
                                                                                                                                               ccgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACt       >  1:1391148/1‑150 (MQ=37)
                                                                                                                                               ccgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCGTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAGCAGGTGAAACt       >  2:2934605/1‑150 (MQ=21)
                                                                                                                                                 gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCCAAGCGGGTTTTTACGTAAAACAGGTGAAACTgc     >  1:1889599/1‑149 (MQ=21)
                                                                                                                                                  ccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAAACAGGTGAAACTGAc    >  1:2556440/1‑150 (MQ=21)
                                                                                                                                                     |                                                                                                                                                    
ATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCG  >  NC_000913/3269820‑3270117

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: