Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,454 T→A Q33L (CAG→CTG)  rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,4540TA96.1% 1443.4 / 30.1 416Q33L (CAG→CTG) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base T (5/11);  new base A (202/197);  total (208/208)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

Only 100 of 417 total aligned reads displayed.

GCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGT  >  NC_000913/2867306‑2867604
                                                                                                                                                    |                                                                                                                                                      
    cgcgCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTcta                                                                                                                                                  <  3:1603666/151‑1 (MQ=255)
      cgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTctac                                                                                                                                                 >  7:1388997/1‑150 (MQ=255)
        aaaaaTAAAGTTCTTCTTCGGCTGTTAACAGTGGTGAATAACCAATCTCACCATGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTaa                                                                                                                                              <  1:551488/150‑1 (MQ=255)
         aaaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAg                                                                                                                                             <  5:874666/151‑1 (MQ=255)
          aaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAgg                                                                                                                                            >  5:132156/1‑151 (MQ=255)
            aTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGcc                                                                                                                                          <  2:1306233/151‑1 (MQ=255)
            aTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGcc                                                                                                                                          >  3:123041/1‑151 (MQ=255)
              aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCt                                                                                                                                         <  6:630359/150‑1 (MQ=255)
               aaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCCGTTCTACTAAGGCCttt                                                                                                                                       >  6:877954/1‑151 (MQ=255)
                aCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCtttt                                                                                                                                      <  4:158377/151‑1 (MQ=255)
                aCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCtttt                                                                                                                                      >  7:308489/1‑151 (MQ=255)
                aCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCtttt                                                                                                                                      <  6:295182/151‑1 (MQ=255)
                  ttcttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGCTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCg                                                                                                                                    >  6:185895/1‑151 (MQ=255)
                  ttcttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCg                                                                                                                                    >  5:222330/1‑151 (MQ=255)
                   tcttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCg                                                                                                                                    <  5:676097/150‑1 (MQ=255)
                   tcttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCg                                                                                                                                    <  2:1656493/150‑1 (MQ=255)
                   tcttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGt                                                                                                                                   >  2:450879/1‑151 (MQ=255)
                    cttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTc                                                                                                                                  >  2:312612/1‑151 (MQ=255)
                     ttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTc                                                                                                                                  >  2:724889/1‑150 (MQ=255)
                       cttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCaa                                                                                                                                >  3:201923/1‑150 (MQ=255)
                         tcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCaaaa                                                                                                                              <  5:579059/150‑1 (MQ=255)
                            gCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCt                                                                                                                          <  7:605181/151‑1 (MQ=255)
                            gCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCt                                                                                                                          >  6:1118895/1‑151 (MQ=255)
                             ccGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCt                                                                                                                          >  8:870198/1‑150 (MQ=255)
                              cGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTc                                                                                                                         >  8:682728/1‑150 (MQ=255)
                              cGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGGTCTAGTAAGGCCTTTACGTCAAACGCCTc                                                                                                                         >  4:82305/1‑150 (MQ=255)
                               gTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCaa                                                                                                                       <  4:366611/151‑1 (MQ=255)
                                ttAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAAc                                                                                                                      >  7:527285/1‑151 (MQ=255)
                                 tAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt                                                                                                                     >  6:1604061/1‑151 (MQ=255)
                                    cAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCg                                                                                                                  >  4:452649/1‑151 (MQ=255)
                                     aGTGCTGAATAACCAAGCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAATGCGTTATCACTGGGGTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGt                                                                                                                 <  3:66844/151‑1 (MQ=255)
                                       tggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtc                                                                                                               >  5:1237549/1‑151 (MQ=255)
                                           gAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcat                                                                                                           <  2:1047108/151‑1 (MQ=255)
                                             aTAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                                                                                                          >  8:105126/1‑150 (MQ=255)
                                                       tcACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCa                                                                                                >  1:90740/1‑150 (MQ=255)
                                                        cACCAAGGTAAAGCTGAGTGGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                                                               <  7:1046186/150‑1 (MQ=255)
                                                        cACCAAGGTAAAGATGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAAAACAGTTCCTCTACGCCCAAATCCTTATCACTTGGTTCCAGTTCTACTAAGGCATTTTCGTCAAAAACATCAACTCCGATATCATAAAATACCGAAt                                                                                               >  6:100850/1‑150 (MQ=255)
                                                         aCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt                                                                                             <  2:516390/151‑1 (MQ=255)
                                                         aCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt                                                                                             <  3:1367020/151‑1 (MQ=255)
                                                         aCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATAACTGGGTTCCAGTTCTACTAAGGCCTTTTAGACAAAAACCTCAACTCCGTTCTAATCAAAATCCGCATCt                                                                                             >  1:1635956/1‑151 (MQ=255)
                                                          ccAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCtt                                                                                            >  8:207556/1‑151 (MQ=255)
                                                            aaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCa                                                                                          <  1:358851/151‑1 (MQ=255)
                                                             aGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAt                                                                                         >  8:1384128/1‑151 (MQ=255)
                                                             aGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAt                                                                                         <  3:926685/151‑1 (MQ=255)
                                                                 aaaGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa                                                                                      <  2:1542495/150‑1 (MQ=255)
                                                                    gCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                  <  8:1118506/151‑1 (MQ=255)
                                                                     cTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                  >  1:1091335/1‑150 (MQ=255)
                                                                     cTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                  >  7:766709/1‑150 (MQ=255)
                                                                      ggAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCAt                                                                                <  2:583227/150‑1 (MQ=255)
                                                                       gAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCa                                                                                 <  8:1463512/149‑1 (MQ=255)
                                                                       gAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                               <  4:234653/151‑1 (MQ=255)
                                                                       gAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCAt                                                                                <  3:1661678/150‑1 (MQ=255)
                                                                            gcgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACtt                                                                          <  6:1410785/151‑1 (MQ=255)
                                                                             cgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACttt                                                                         <  1:1441001/151‑1 (MQ=255)
                                                                              gTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATATTCATTTAAATAATGAACttt                                                                         >  6:1533028/1‑150 (MQ=255)
                                                                                ccAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCTAATTCCGCATCTTCATTTAAATCATGAACTTTc                                                                        >  3:1540681/1‑149 (MQ=255)
                                                                                ccAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAg                                                                      <  5:1655048/151‑1 (MQ=255)
                                                                                   acacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGt                                                                   >  1:76347/1‑151 (MQ=255)
                                                                                    cacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGt                                                                   >  6:1194966/1‑150 (MQ=255)
                                                                                     acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTAt                                                                 >  3:348262/1‑151 (MQ=255)
                                                                                       acGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTc                                                               >  5:251303/1‑151 (MQ=255)
                                                                                         gCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTg                                                             <  8:1610128/151‑1 (MQ=255)
                                                                                              gtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAATTCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTc                                                         >  5:1184997/1‑150 (MQ=255)
                                                                                              gtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTc                                                         <  6:1304292/150‑1 (MQ=255)
                                                                                               tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAt                                                       >  1:1387256/1‑151 (MQ=255)
                                                                                                  cTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATa                                                      >  6:78364/1‑149 (MQ=255)
                                                                                                    cccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACAAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAgg                                                   <  4:1395070/150‑1 (MQ=255)
                                                                                                     ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTg                                                 <  1:1192985/151‑1 (MQ=255)
                                                                                                      cTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCGAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTAAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTgg                                                >  8:31248/1‑151 (MQ=255)
                                                                                                      cTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTgg                                                >  5:1574751/1‑151 (MQ=255)
                                                                                                          gATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTGCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTcc                                            >  8:408872/1‑151 (MQ=255)
                                                                                                           aTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCt                                           >  7:1602836/1‑151 (MQ=255)
                                                                                                           aTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCt                                           >  7:1134774/1‑151 (MQ=255)
                                                                                                               cAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTAccc                                       >  1:1046333/1‑151 (MQ=255)
                                                                                                               cAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTAccc                                       <  3:502043/151‑1 (MQ=255)
                                                                                                               cAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTAccc                                       <  4:309554/151‑1 (MQ=255)
                                                                                                               cAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGCTGGCTCCTAccc                                       >  5:932216/1‑151 (MQ=255)
                                                                                                                aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCg                                      <  4:103532/151‑1 (MQ=255)
                                                                                                                  ttCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                    >  2:1368013/1‑151 (MQ=255)
                                                                                                                   tCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                    >  8:808037/1‑150 (MQ=255)
                                                                                                                    ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGa                                   <  8:77287/150‑1 (MQ=255)
                                                                                                                      tctTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATcc                                >  4:1208410/1‑151 (MQ=255)
                                                                                                                       ctTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATATTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATccc                               >  2:1588203/1‑151 (MQ=255)
                                                                                                                         tCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCtt                             >  3:578386/1‑151 (MQ=255)
                                                                                                                          gcgCCCAATCCTTAATACTGGGGTGCAGTTCCACTTACGCCCTTTCCGCCAAAACCCCCACCCCGTTCCCAACAAAATCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCtt                             <  3:1041811/148‑1 (MQ=255)
                                                                                                                             ccAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACg                         >  5:757824/1‑151 (MQ=255)
                                                                                                                                aaTCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGa                       >  3:784044/1‑150 (MQ=255)
                                                                                                                                 aTCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGaa                      >  5:1004866/1‑150 (MQ=255)
                                                                                                                                    gTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTATCCGTGATCCCTTGACGGAACAtt                  <  1:239121/151‑1 (MQ=255)
                                                                                                                                    gTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAtt                  <  2:881022/151‑1 (MQ=255)
                                                                                                                                    gTTATCACTGGGTTCCAGTTCTACTAAGGACTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTAATTTAAATCATTAACTTTCAGCGTATTATGACTCATAAGGTGACTCCTACCCGTGATCCCTTGACTGAACaa                   >  1:1327713/1‑149 (MQ=255)
                                                                                                                                     ttATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAtt                  >  2:1185057/1‑150 (MQ=255)
                                                                                                                                     ttATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAAATTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTc                 >  4:1126171/1‑151 (MQ=255)
                                                                                                                                      tATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTc                 <  5:358940/150‑1 (MQ=255)
                                                                                                                                      tATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTc                 <  7:1309404/150‑1 (MQ=255)
                                                                                                                                       aTCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCa                <  7:394011/150‑1 (MQ=255)
                                                                                                                                         cACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAg              <  8:1010467/150‑1 (MQ=255)
                                                                                                                                           cTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAACCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGCCTCATAAGGTGGCTCCTCCCCGTGCTCCCTTGACGGAACATTCAAGCa            >  4:591071/1‑150 (MQ=255)
                                                                                                                                               gTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGc       <  8:1251297/151‑1 (MQ=255)
                                                                                                                                                 tCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGcc      <  1:1285795/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGT  >  NC_000913/2867306‑2867604

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