Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,299,973 T→A 8.2% intergenic (+34/+36) slyB → / ← nhaA‑I outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,299,9730TA 8.2% 1695.6 / 97.0 547intergenic (+34/+36)slyB/nhaA‑Iouter membrane lipoprotein/Na(+)/H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (266/235);  new base A (27/18);  total (293/254)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 586 total aligned reads displayed.

AAGATGATGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCT  >  NC_002947/1299824‑1300122
                                                                                                                                                     |                                                                                                                                                       
  gatgatGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGgctgc                                                                                                                                                       >  1:2180913/1‑150 (MQ=255)
     gatgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCtgctgc                                                                                                                                                    <  1:4791545/150‑1 (MQ=255)
     gatgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCtgctgc                                                                                                                                                    <  2:2908857/150‑1 (MQ=255)
      atgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGcc                                                                                                                                                   <  1:5460040/150‑1 (MQ=255)
       tgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCa                                                                                                                                                  >  1:7100432/1‑150 (MQ=255)
          cagcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTc                                                                                                                                              >  2:4011428/1‑150 (MQ=255)
             caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCggg                                                                                                                                           <  1:1323886/150‑1 (MQ=255)
             caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCggg                                                                                                                                           <  2:6155388/150‑1 (MQ=255)
             caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCggg                                                                                                                                           <  1:6484549/150‑1 (MQ=255)
             caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCggg                                                                                                                                           <  2:332317/150‑1 (MQ=255)
             caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCggg                                                                                                                                           <  1:8650969/150‑1 (MQ=255)
                 cgcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGttt                                                                                                                                       <  2:7071447/150‑1 (MQ=255)
                 cgcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGttt                                                                                                                                       <  1:5914035/150‑1 (MQ=255)
                      cTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCa                                                                                                                                  >  1:5161034/1‑150 (MQ=255)
                       tATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCAt                                                                                                                                 <  1:976301/150‑1 (MQ=255)
                          gtgCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAAATGGCTGCTGCCAGTCGGGGTTTTTTCATGcc                                                                                                                              >  1:8398325/1‑150 (MQ=255)
                             cAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCa                                                                                                                           >  2:1448635/1‑150 (MQ=255)
                              aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAg                                                                                                                          <  1:7811305/150‑1 (MQ=255)
                                  cGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCaa                                                                                                                      <  1:4618031/150‑1 (MQ=255)
                                      cAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCa                                                                                                                  <  2:528745/150‑1 (MQ=255)
                                        ggAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAgg                                                                                                                >  2:2794672/1‑150 (MQ=255)
                                         gagaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGc                                                                                                               <  1:8041696/150‑1 (MQ=255)
                                             aCGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATg                                                                                                           <  1:6602918/150‑1 (MQ=255)
                                             aCGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATg                                                                                                           >  1:6975098/1‑150 (MQ=255)
                                              cGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTCATGCCTCAGCCAAGTCAGGCAATGCt                                                                                                         >  2:5201169/1‑150 (MQ=255)
                                               gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCt                                                                                                         >  2:2543803/1‑150 (MQ=255)
                                               gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCt                                                                                                         >  2:556832/1‑150 (MQ=255)
                                                aaaTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTg                                                                                                        <  2:5789339/150‑1 (MQ=255)
                                                aaaTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTg                                                                                                        <  2:7390080/150‑1 (MQ=255)
                                                 aaTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTgg                                                                                                       >  1:4323415/1‑150 (MQ=255)
                                                  aTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGt                                                                                                      >  2:7940369/1‑150 (MQ=255)
                                                   tCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGtt                                                                                                     >  2:1725652/1‑150 (MQ=255)
                                                     ttCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGtttt                                                                                                   >  1:5365443/1‑150 (MQ=255)
                                                     ttCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGtttt                                                                                                   >  1:3332148/1‑150 (MQ=255)
                                                       ccGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGc                                                                                                 >  2:4829836/1‑150 (MQ=255)
                                                        cGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCg                                                                                                >  2:2686717/1‑150 (MQ=255)
                                                                 gACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCAt                                                                                       <  2:3165791/150‑1 (MQ=255)
                                                                 gACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCAt                                                                                       <  2:5671413/150‑1 (MQ=255)
                                                                    cGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGc                                                                                    >  1:6422480/1‑150 (MQ=255)
                                                                     gAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCg                                                                                   >  1:4882428/1‑150 (MQ=255)
                                                                     gAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCg                                                                                   <  2:7219840/150‑1 (MQ=255)
                                                                       gTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGt                                                                                 <  1:2743955/150‑1 (MQ=255)
                                                                           gcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACg                                                                            >  1:3538948/1‑150 (MQ=255)
                                                                           gcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACg                                                                            <  2:6290129/150‑1 (MQ=255)
                                                                             atcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGc                                                                          >  1:3984436/1‑150 (MQ=255)
                                                                               catGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGt                                                                        <  1:7533084/150‑1 (MQ=255)
                                                                                atGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTa                                                                       <  1:7276689/150‑1 (MQ=255)
                                                                                atGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTa                                                                       <  2:7630900/150‑1 (MQ=255)
                                                                                atGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTa                                                                       >  1:8215408/1‑150 (MQ=255)
                                                                                atGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTa                                                                       <  1:2289587/150‑1 (MQ=255)
                                                                                 tGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTg                                                                                            >  2:3391904/1‑129 (MQ=255)
                                                                                 tGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAg                                                                      <  1:3802708/150‑1 (MQ=255)
                                                                                  gACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGc                                                                     >  2:8577344/1‑150 (MQ=255)
                                                                                  gACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGc                                                                     >  2:105932/1‑150 (MQ=255)
                                                                                     tGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGa                                                                  <  1:2940362/150‑1 (MQ=255)
                                                                                       tAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATc                                                                <  1:2661336/150‑1 (MQ=255)
                                                                                       tAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATc                                                                >  1:3023400/1‑150 (MQ=255)
                                                                                       tAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATc                                                                >  1:3440033/1‑150 (MQ=255)
                                                                                       tAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATc                                                                >  2:2422883/1‑150 (MQ=255)
                                                                                       tAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATc                                                                >  1:6328712/1‑150 (MQ=255)
                                                                                        aGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCa                                                               >  1:1008599/1‑150 (MQ=255)
                                                                                          aCGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc                                                             >  2:8472362/1‑150 (MQ=255)
                                                                                          aCGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc                                                             >  1:8772717/1‑150 (MQ=255)
                                                                                           cGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCg                                                            <  2:6438619/150‑1 (MQ=255)
                                                                                            ggTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgc                                                           <  1:5609688/150‑1 (MQ=255)
                                                                                             gTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgca                                                          <  2:6781954/150‑1 (MQ=255)
                                                                                              tACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGATCGGCTACGCCGTCACCgcaa                                                         <  2:543528/150‑2 (MQ=255)
                                                                                                  aGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCaaa                                                     >  2:3483555/1‑150 (MQ=255)
                                                                                                  aGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCaaa                                                     <  1:2488186/150‑1 (MQ=255)
                                                                                                    tCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAga                                                   >  2:3760830/1‑150 (MQ=255)
                                                                                                    tCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAga                                                   >  1:1531984/1‑150 (MQ=255)
                                                                                                         tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGTCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCTATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCATCAAAgaacgta                                              <  2:3908867/150‑3 (MQ=255)
                                                                                                         tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaa                                              >  2:8074852/1‑150 (MQ=255)
                                                                                                              cACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGt                                         <  1:5763429/150‑1 (MQ=255)
                                                                                                              cACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGNGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGt                                         <  2:2234652/150‑1 (MQ=255)
                                                                                                                cTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTaa                                       <  2:1921693/150‑1 (MQ=255)
                                                                                                                    ccGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAAt                                   <  2:6654298/150‑1 (MQ=255)
                                                                                                                     cGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATg                                  >  1:3131437/1‑150 (MQ=255)
                                                                                                                      ctAGTGCAGTAATAAAAAACCCCGACATGGCTGCGGCCAGTCGGGGTTTTTTCATGCCTCCGCCAGGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGc                                 <  2:5999277/149‑1 (MQ=255)
                                                                                                                           gCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATg                            <  2:1247354/150‑1 (MQ=255)
                                                                                                                           gCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATg                            <  2:7460531/150‑1 (MQ=255)
                                                                                                                            cAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGc                           >  1:1854181/1‑150 (MQ=255)
                                                                                                                             aGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacga                                               <  2:4395893/129‑1 (MQ=255)
                                                                                                                             aGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCg                          >  1:2039323/1‑150 (MQ=255)
                                                                                                                               taataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCACTGCAGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGt                        >  1:4212207/1‑150 (MQ=255)
                                                                                                                               taataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGt                        <  2:6926745/150‑1 (MQ=255)
                                                                                                                               taataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGt                        >  2:8320774/1‑150 (MQ=255)
                                                                                                                                aataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTc                       <  1:1924247/150‑1 (MQ=255)
                                                                                                                                   aaaaaaCCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCGTTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCAt                    <  2:4763259/150‑1 (MQ=255)
                                                                                                                                      aaaCCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGcc                 >  2:6534204/1‑150 (MQ=255)
                                                                                                                                       aaCCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCTGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCg                <  1:6290570/150‑1 (MQ=255)
                                                                                                                                        aCCCCGACTGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCgg               >  1:5240366/1‑150 (MQ=255)
                                                                                                                                         ccccGACTGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc              >  2:4997476/1‑150 (MQ=255)
                                                                                                                                          cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCg             >  1:6765523/1‑150 (MQ=255)
                                                                                                                                          cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCg             >  2:7555914/1‑150 (MQ=255)
                                                                                                                                          cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCg             >  1:7355952/1‑150 (MQ=255)
                                                                                                                                              aCATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCaaa                                                     <  1:7599845/106‑1 (MQ=255)
                                                                                                                                                aTGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCa       >  2:87627/1‑150 (MQ=255)
                                                                                                                                                aTGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCa       >  1:6090216/1‑150 (MQ=255)
                                                                                                                                                 tGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCAc      <  1:1209227/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AAGATGATGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCT  >  NC_002947/1299824‑1300122

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