Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,188,777 A→G D1143G (GAC→GGC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,188,7770AG99.0% 356.3 / ‑4.6 98D1143G (GAC→GGC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (47/50);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAA  >  NC_000913/4188637‑4188924
                                                                                                                                            |                                                                                                                                                   
aGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGa                                                                                                                                            >  1:1658741/1‑150 (MQ=255)
aGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGa                                                                                                                                            >  2:1048201/1‑150 (MQ=255)
   gcAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCa                                                                                                                                        <  2:869487/151‑1 (MQ=255)
        aCTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCAcgtcg                                                                                                                                   <  1:1002688/151‑1 (MQ=255)
          ttCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACgtcgtc                                                                                                                                 >  2:1543360/1‑151 (MQ=255)
           tCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTcc                                                                                                                                >  1:1180470/1‑151 (MQ=255)
             cTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTcc                                                                                                                                >  1:290591/1‑149 (MQ=255)
             cTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCg                                                                                                                               >  1:960614/1‑150 (MQ=255)
               gCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGaa                                                                                                                             >  1:1477834/1‑150 (MQ=255)
               gCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGaa                                                                                                                             <  2:88538/150‑1 (MQ=255)
               gCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGaa                                                                                                                             <  2:982685/150‑1 (MQ=255)
               gCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGa                                                                                                                              <  2:154612/149‑1 (MQ=255)
                ccGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGaa                                                                                                                             >  2:1240655/1‑149 (MQ=255)
                ccGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGaa                                                                                                                             >  2:1612179/1‑149 (MQ=255)
                 cGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGaa                                                                                                                             >  1:1069374/1‑148 (MQ=255)
                  gggTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAagag                                                                                                                         <  2:412556/151‑1 (MQ=255)
                   ggTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGc                                                                                                                        <  2:432488/151‑1 (MQ=255)
                    gTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGcc                                                                                                                       >  1:899384/1‑151 (MQ=255)
                    gTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGcc                                                                                                                       >  2:1023664/1‑151 (MQ=255)
                      aaaGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCgg                                                                                                                     >  1:190094/1‑151 (MQ=255)
                      aaaGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCg                                                                                                                      <  1:1000785/150‑1 (MQ=255)
                         gCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCa                                                                                                                   >  2:1344917/1‑150 (MQ=255)
                          cGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAAt                                                                                                                 >  1:840609/1‑151 (MQ=255)
                           gATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAAt                                                                                                                 >  2:451110/1‑150 (MQ=255)
                           gATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATc                                                                                                                <  2:1176867/151‑1 (MQ=255)
                           gATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATc                                                                                                                >  2:748398/1‑151 (MQ=255)
                              tGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCt                                                                                                              <  2:597795/150‑1 (MQ=255)
                               gTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTgg                                                                                                            <  2:1412310/151‑1 (MQ=255)
                               gTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTgg                                                                                                            >  1:1390352/1‑151 (MQ=255)
                               gTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTg                                                                                                             <  1:732239/150‑1 (MQ=255)
                               gTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTg                                                                                                             <  2:575585/150‑1 (MQ=255)
                                 tCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGc                                                                                                           >  1:1297782/1‑150 (MQ=255)
                                  cAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCt                                                                                                          <  2:550898/150‑1 (MQ=255)
                                  cAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCt                                                                                                          >  1:931371/1‑150 (MQ=255)
                                      tGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGaaa                                                                                                      >  2:1041799/1‑150 (MQ=255)
                                      tGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAAt                                                                                                     >  1:692715/1‑151 (MQ=255)
                                       ggAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATc                                                                                                    >  1:1597828/1‑151 (MQ=255)
                                              ggCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGt                                                                                              >  2:1422219/1‑150 (MQ=255)
                                                 gTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATc                                                                                           <  1:1582165/150‑1 (MQ=255)
                                                  tACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGt                                                                                         >  1:1062300/1‑151 (MQ=255)
                                                   aCAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGtt                                                                                        <  1:1015059/151‑1 (MQ=255)
                                                   aCAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGt                                                                                         <  1:138224/150‑1 (MQ=255)
                                                      gATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTcc                                                                                     <  1:1092641/151‑1 (MQ=255)
                                                       aTCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGTGGTATCGTTTCCt                                                                                    >  2:710439/1‑151 (MQ=255)
                                                       aTCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCt                                                                                    <  2:996817/151‑1 (MQ=255)
                                                       aTCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCt                                                                                    <  1:522337/151‑1 (MQ=255)
                                                           gCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCgg                                                                                >  1:1474346/1‑151 (MQ=255)
                                                           gCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCgg                                                                                <  1:356085/151‑1 (MQ=255)
                                                                   gACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGaa                                                                         <  2:1415282/150‑1 (MQ=255)
                                                                   gACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGaa                                                                         <  2:1393428/150‑1 (MQ=255)
                                                                    acacCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAAc                                                                       <  2:28538/151‑1 (MQ=255)
                                                                     cacCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAAcc                                                                      >  1:881720/1‑151 (MQ=255)
                                                                       cccTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCa                                                                     <  2:479231/150‑1 (MQ=255)
                                                                       cccTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCa                                                                     >  1:1059170/1‑150 (MQ=255)
                                                                       cccTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCa                                                                     >  2:1440479/1‑150 (MQ=255)
                                                                        ccTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaaa                                                                   <  1:1401270/151‑1 (MQ=255)
                                                                        ccTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaa                                                                    >  1:1355781/1‑150 (MQ=255)
                                                                        ccTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaa                                                                    >  2:1433288/1‑150 (MQ=255)
                                                                         cTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaa                                                                    >  2:364783/1‑149 (MQ=255)
                                                                           ggCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGt                                                                <  1:1046026/151‑1 (MQ=255)
                                                                             cgcgTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTa                                                               >  1:840981/1‑150 (MQ=255)
                                                                              gcgTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTaaa                                                             <  2:386479/151‑1 (MQ=255)
                                                                              gcgTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTaaa                                                             <  2:1072716/151‑1 (MQ=255)
                                                                                gTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcg                                                           >  1:522122/1‑151 (MQ=255)
                                                                                  aTTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcgtc                                                         <  1:131879/151‑1 (MQ=255)
                                                                                  aTTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcgtc                                                         <  2:1658849/151‑1 (MQ=255)
                                                                                     ccGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACgtcgtc                                                      <  1:87925/151‑1 (MQ=255)
                                                                                       gCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTg                                                    >  2:809303/1‑151 (MQ=255)
                                                                                       gCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTg                                                    <  2:1699071/151‑1 (MQ=255)
                                                                                       gCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTg                                                    >  2:501867/1‑151 (MQ=255)
                                                                                          ggAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGtt                                                 >  2:808819/1‑151 (MQ=255)
                                                                                           gAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTa                                                <  2:833205/151‑1 (MQ=255)
                                                                                            aaTCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTAt                                               <  1:1698793/151‑1 (MQ=255)
                                                                                             aTCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATc                                              >  2:99459/1‑151 (MQ=255)
                                                                                             aTCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATc                                              <  1:1612081/151‑1 (MQ=255)
                                                                                               ccGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCAc                                            >  2:863133/1‑151 (MQ=255)
                                                                                               ccGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCAc                                            >  1:1238521/1‑151 (MQ=255)
                                                                                                 ggcggTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCAccc                                          <  2:1197183/151‑1 (MQ=255)
                                                                                                      tACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGt                                      >  2:1107118/1‑150 (MQ=255)
                                                                                                       aCCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTa                                     <  1:1344824/150‑1 (MQ=255)
                                                                                                          aaGGACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACg                                 >  1:475780/1‑151 (MQ=255)
                                                                                                             gACATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTa                              <  2:1069459/151‑1 (MQ=255)
                                                                                                              aCATCACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGt                               >  1:893788/1‑149 (MQ=255)
                                                                                                                  cACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGa                          <  2:679474/150‑1 (MQ=255)
                                                                                                                  cACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGAt                         >  2:859883/1‑151 (MQ=255)
                                                                                                                  cACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGAt                         <  2:542859/151‑1 (MQ=255)
                                                                                                                  cACCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGAt                         <  2:899452/151‑1 (MQ=255)
                                                                                                                   aCCGGTGGTCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATc                        <  2:1390444/151‑1 (MQ=255)
                                                                                                                   aCCGGTGGTCTGCCGCGCGTTGCGGGCCTGGTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGAGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGGAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATc                        >  2:1423654/1‑151 (MQ=255)
                                                                                                                      ggtggtCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGt                     <  2:1092335/151‑1 (MQ=255)
                                                                                                                         ggtCTGCCGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTAc                   <  2:1062381/150‑1 (MQ=255)
                                                                                                                               ccGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAag               >  2:1421915/1‑148 (MQ=255)
                                                                                                                               ccGCGCGTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAgaga            <  2:360267/151‑1 (MQ=255)
                                                                                                                                   gcgTTGCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAgatgat        >  1:1667762/1‑151 (MQ=255)
                                                                                                                                        gCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCg    <  2:1696851/150‑1 (MQ=255)
                                                                                                                                        gCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCg    <  1:1686232/150‑1 (MQ=255)
                                                                                                                                        gCGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCg    <  1:1240352/150‑1 (MQ=255)
                                                                                                                                         cGGGCCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGaa  <  1:1500960/151‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
AGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAA  >  NC_000913/4188637‑4188924

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: