Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,256,890 A→T 5.1% noncoding (431/451 nt) PP_mr14 ← cspA regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,256,8900AT 5.1% 385.3 / 8.3 117noncoding (431/451 nt)PP_mr14cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base A (53/58);  new base T (2/4);  total (55/62)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 121 total aligned reads displayed.

TACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCTTTG  >  NC_002947/1256741‑1257038
                                                                                                                                                     |                                                                                                                                                    
tACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGaa                                                                                                                                                      <  2:652955/150‑1 (MQ=255)
tACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGaa                                                                                                                                                      <  1:1389155/150‑1 (MQ=255)
tACCAACAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGaa                                                                                                                                                      <  2:869776/150‑1 (MQ=255)
  ccATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAgg                                                                                                                                                    >  2:1350136/1‑150 (MQ=255)
   caacAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGc                                                                                                                                                   >  1:2407065/4‑150 (MQ=255)
    aTCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGAATAAAAAAACCCGCCAGAAGGCg                                                                                                                                                  >  2:490154/1‑150 (MQ=255)
     tCAGCGAATACTCCCCTGCCATACCGGGTTTGCCCGCACACAGCCAATCCGCCGTACTTGAAGCCTGTGCAGTCCCTGGGGGAGCGGGCGAGCCCGCCAAAAGGCCAGCCCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                 <  1:528877/150‑1 (MQ=255)
     tCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAGCCCGCCAGAAGGCgg                                                                                                                                                 >  2:1682598/1‑150 (MQ=255)
     tCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                 <  1:2141397/150‑1 (MQ=255)
     tCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                 <  2:1065793/150‑1 (MQ=255)
          gAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                 >  1:1149064/1‑145 (MQ=255)
          gAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                 <  2:1149064/145‑1 (MQ=255)
           aaGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGTCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttt                                                                                                                                           <  1:463327/150‑1 (MQ=255)
             gACTGACCTGCCAGACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttttt                                                                                                                                         <  1:649550/150‑1 (MQ=255)
              aCTGACCTGCCATACCGAGTTTGCCCGCACACAGCGAACCCGTCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGttttttta                                                                                                                                        <  1:633688/150‑2 (MQ=37)
                     tGCCATAGCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                                                                                                                                 <  2:1852587/150‑1 (MQ=255)
                        cATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGa                                                                                                                              >  2:1026777/1‑150 (MQ=255)
                         atcccGAGTTTGCCCGCACACAGCGACTGCGCCGTACTGGTGGCCTGTGCAGTCCTTGTGGGAGCGGGCGAGCCCGCAAAGAGGCCAGCGCCTATTTCTAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGaa                                                                                                                             <  1:140285/147‑1 (MQ=255)
                          tACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAg                                                                                                                            >  2:1223945/1‑150 (MQ=255)
                           aCCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                                                                                                                           <  1:416274/150‑1 (MQ=255)
                             cGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACACGCCAGAAGGGGGGTGTTTTTGAAAGCTAGGCAGCTa                                                                                                                         >  2:911087/1‑150 (MQ=255)
                                 tttGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATc                                                                                                                     <  2:1154166/150‑1 (MQ=255)
                                     cccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACt                                                                                                                 <  2:1503163/150‑1 (MQ=255)
                                      ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGtttttttat                                                                                                                                       >  1:1786512/1‑125 (MQ=37)
                                      ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACtt                                                                                                                >  2:1986611/1‑150 (MQ=255)
                                      ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACtt                                                                                                                >  2:659962/1‑150 (MQ=255)
                                      ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACtt                                                                                                                >  1:881461/1‑150 (MQ=255)
                                              aGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATc                                                                                                        <  1:1920275/150‑1 (MQ=255)
                                                cGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGCCTAAGTCGATCta                                                                                                      <  1:1008423/150‑2 (MQ=255)
                                                  aCTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGTGCGAGCCCGCGAAGAGGCCAGCGCATCACGCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGaa                                                                                                    <  1:761679/150‑1 (MQ=38)
                                                   cggcgcCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAAc                                                                                                   <  2:451049/148‑1 (MQ=255)
                                                          gTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCa                                                                                            <  1:770526/150‑1 (MQ=255)
                                                            aCTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGc                                                                                          >  1:1686143/1‑150 (MQ=255)
                                                             cTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATAGACTTAGGCGATCTGAACTTCTTCAGcc                                                                                         <  1:1644753/150‑1 (MQ=255)
                                                              ttGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                                                                                                                                 <  2:2190746/109‑1 (MQ=255)
                                                              ttGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATc                                                                                                                     <  2:532806/121‑1 (MQ=255)
                                                              ttGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAAACCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATc                                                                                                                     >  1:532806/1‑121 (MQ=255)
                                                                gAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGaa                                                                                                                             <  1:1371518/111‑1 (MQ=255)
                                                                gAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGc                                                                                      >  1:987712/1‑150 (MQ=255)
                                                                    ccTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGc                                                                                  >  1:28441/1‑150 (MQ=255)
                                                                     cTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGcc                                                                                 >  2:932196/1‑150 (MQ=255)
                                                                      tgtgCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCt                                                                                <  1:1026777/150‑1 (MQ=255)
                                                                      tgtgCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCt                                                                                >  2:830973/1‑150 (MQ=255)
                                                                      tgtgCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCt                                                                                >  1:1549093/1‑150 (MQ=255)
                                                                       gtgCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCtt                                                                               >  2:670959/1‑150 (MQ=255)
                                                                          cAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCt                                                                            >  1:1078277/1‑150 (MQ=255)
                                                                           aGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTg                                                                           <  2:2407065/150‑1 (MQ=255)
                                                                           aGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTg                                                                           >  2:2365017/1‑150 (MQ=255)
                                                                            gTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTgg                                                                          <  2:2073692/150‑1 (MQ=255)
                                                                             tCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTa                                                                                                               >  1:1112619/1‑112 (MQ=255)
                                                                             tCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTa                                                                                                               <  2:1112619/112‑1 (MQ=255)
                                                                                cTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                                                                                                                           >  2:1711661/1‑97 (MQ=255)
                                                                                cTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                                                                                                                           <  1:1711661/97‑1 (MQ=255)
                                                                                cTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcgc                                                                      <  2:28441/150‑1 (MQ=255)
                                                                                 tgtgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCttt                                                                              >  2:165152/1‑141 (MQ=255)
                                                                                 tgtgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCttt                                                                              <  1:165152/141‑1 (MQ=255)
                                                                                  gtgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCgg                                                                    <  1:659962/150‑1 (MQ=255)
                                                                                   tgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCggg                                                                   >  2:615221/1‑150 (MQ=255)
                                                                                     ttcgCGGGCTCGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTa                                                                 <  2:1125508/147‑1 (MQ=255)
                                                                                              gaccccGCGTGGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGaaa                                                        <  1:104362/147‑1 (MQ=255)
                                                                                                  ccGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                    >  2:1820615/1‑150 (MQ=255)
                                                                                                     cgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAg                                                 <  2:1453474/150‑1 (MQ=255)
                                                                                                         gagGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTg                                             <  2:2085716/150‑1 (MQ=255)
                                                                                                         gagGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTg                                             <  2:1549093/150‑1 (MQ=255)
                                                                                                              cAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                    <  1:1722458/138‑1 (MQ=255)
                                                                                                              cAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                    >  2:1722458/1‑138 (MQ=255)
                                                                                                                gcgcCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCcttct                                      >  1:775415/1‑150 (MQ=255)
                                                                                                                  gcCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCAt                                                                                    >  1:2007437/1‑102 (MQ=255)
                                                                                                                    cTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTg                                             >  2:1836005/1‑139 (MQ=255)
                                                                                                                     tATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCt                                                                                                       <  2:738526/80‑1 (MQ=255)
                                                                                                                     tATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCt                                                                                                       >  1:738526/1‑80 (MQ=255)
                                                                                                                     tATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAg                                 <  1:1820615/150‑1 (MQ=255)
                                                                                                                     tATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAg                                 <  1:1223945/150‑1 (MQ=255)
                                                                                                                      aTTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCctt                                        <  1:1540547/142‑1 (MQ=255)
                                                                                                                      aTTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCctt                                        >  2:1540547/1‑142 (MQ=255)
                                                                                                                      aTTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAgg                                <  1:1359519/150‑1 (MQ=255)
                                                                                                                       tttCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAggg                               >  2:1794891/1‑150 (MQ=255)
                                                                                                                        ttCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                    <  1:463297/128‑1 (MQ=255)
                                                                                                                              cAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGaa                        <  1:2365017/150‑1 (MQ=255)
                                                                                                                                 gCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCgg                                                                    <  1:126735/103‑1 (MQ=255)
                                                                                                                                 gCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCgg                                                                    >  2:126735/1‑103 (MQ=255)
                                                                                                                                 gCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACCAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGcc                     <  2:1390986/150‑1 (MQ=255)
                                                                                                                                   aTAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGt                   >  2:1272830/1‑150 (MQ=255)
                                                                                                                                     aaaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAgg                                >  1:623328/1‑135 (MQ=255)
                                                                                                                                     aaaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAgg                                <  2:623328/135‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAg                >  2:1245265/1‑150 (MQ=255)
                                                                                                                                       aaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGc               >  1:2031719/1‑150 (MQ=255)
                                                                                                                                        aaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCt              >  2:538900/1‑150 (MQ=255)
                                                                                                                                            cccGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTgg           >  1:692506/1‑149 (MQ=255)
                                                                                                                                            cccGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTgg           <  2:692506/149‑1 (MQ=255)
                                                                                                                                             ccGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAgg                                >  2:2192704/1‑127 (MQ=255)
                                                                                                                                             ccGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAgg                                <  1:2192704/127‑1 (MQ=255)
                                                                                                                                             ccGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGAt         >  1:1844374/1‑150 (MQ=255)
                                                                                                                                              cGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATg        >  1:1165262/1‑150 (MQ=255)
                                                                                                                                               gCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAgg                                <  1:311987/125‑1 (MQ=255)
                                                                                                                                               gCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAgg                                >  2:311987/1‑125 (MQ=255)
                                                                                                                                               gCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATgg       >  1:1500758/1‑150 (MQ=255)
                                                                                                                                               gCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATgg       >  2:1305229/1‑150 (MQ=255)
                                                                                                                                                 cAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCt     <  2:1688644/150‑1 (MQ=255)
                                                                                                                                                   gAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCttt   <  1:1272830/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCTTTG  >  NC_002947/1256741‑1257038

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑