Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,966,388 A→T 6.5% intergenic (‑401/+29) PP_5232 ← / ← amtB hypothetical protein/ammonium transporter AmtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,966,3880AT 6.5% 211.2 / 2.9 71intergenic (‑401/+29)PP_5232/amtBhypothetical protein/ammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base A (25/41);  new base T (2/3);  total (27/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCTTCGGCGCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATG  >  NC_002947/5966238‑5966537
                                                                                                                                                      |                                                                                                                                                      
ccTTCGGCGCGCATTCTAGAGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAATGCCCGGCa                                                                                                                                                        >  2:1731712/1‑150 (MQ=255)
 cTTCGGCGCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCa                                                                                                                                                        >  2:1309575/1‑150 (MQ=255)
   tCGGCGCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACCAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAg                                                                                                                                                      >  2:1799293/1‑150 (MQ=255)
      gcgcgcATTCTATCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTCTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCg                                                                                                                                                   <  2:1055629/150‑1 (MQ=255)
      gcgcgcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCg                                                                                                                                                   <  2:1251295/150‑1 (MQ=255)
       cgcgcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCgg                                                                                                                                                  <  1:913955/150‑1 (MQ=255)
           cATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCaa                                                                                                                                              <  1:209097/150‑1 (MQ=255)
            aTTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCTTGCCGGGCatt                                                                                                                                             <  1:693715/150‑3 (MQ=255)
            aTTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCGGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAg                                                                                                                                             >  1:957105/1‑150 (MQ=255)
              tCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCTTGCCGGGCatttt                                                                                                                                           <  1:823792/150‑5 (MQ=255)
               cTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGttt                                                                                                                                          <  1:1491072/150‑1 (MQ=255)
                   cgcgcCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAATCCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTAc                                                                                                                                      <  2:1088452/150‑1 (MQ=255)
                      gcCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACcgc                                                                                                                                   <  1:708829/150‑1 (MQ=255)
                        ccGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGt                                                                                                                                 >  1:628551/1‑150 (MQ=255)
                        ccGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGt                                                                                                                                 >  1:1293691/1‑150 (MQ=255)
                           tGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGc                                                                                                                              <  1:451205/150‑1 (MQ=255)
                           tGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGc                                                                                                                              <  2:655459/150‑1 (MQ=255)
                           tGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGc                                                                                                                              >  1:1245295/1‑150 (MQ=255)
                            ggAAAAAGCACACCCCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCa                                                                                                                             <  2:1480127/150‑1 (MQ=255)
                            ggAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCa                                                                                                                             <  1:1324104/150‑1 (MQ=255)
                             gAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAg                                                                                                                            >  2:723520/1‑150 (MQ=255)
                               aaaaGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTTCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGtt                                                                                                                          >  1:517552/1‑150 (MQ=255)
                                 aagcaaagcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTAc                                                                                                                        <  2:1138088/150‑1 (MQ=255)
                                   gcaaagcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGAACAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAg                                                                                                                      <  2:1215116/150‑1 (MQ=255)
                                    caaagcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAgg                                                                                                                     <  2:597768/150‑1 (MQ=255)
                                    caaagcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAgg                                                                                                                     <  1:1738002/150‑1 (MQ=255)
                                       agcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTAACTTCGGCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTg                                                                                                                  <  1:17700/150‑1 (MQ=255)
                                            cATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGcc                                                                                                             <  1:1529372/150‑1 (MQ=255)
                                                    cccTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGt                                                                                                     >  2:1138247/1‑150 (MQ=255)
                                                       tGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTtgt                                                                                                  <  1:1309575/150‑1 (MQ=255)
                                                         cGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGc                                                                                                <  1:883807/150‑1 (MQ=255)
                                                            cTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGa                                                                                             <  2:1575372/150‑1 (MQ=255)
                                                             ttGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAg                                                                                            >  2:1685802/1‑150 (MQ=255)
                                                              tGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGc                                                                                           <  1:501028/150‑1 (MQ=255)
                                                                tAATTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCa                                                                                        <  1:1433511/150‑1 (MQ=255)
                                                                tAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCa                                                                                        >  1:203871/1‑150 (MQ=255)
                                                                   gTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTCCCGCGTAGTACTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGt                                                                                     >  1:1452631/1‑150 (MQ=255)
                                                                    ttttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTc                                                                                    >  1:1054359/1‑150 (MQ=255)
                                                                      ttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGa                                                                                  <  1:723520/150‑1 (MQ=255)
                                                                       tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGGGCGCGTTGTGCTGAGCCAGGTCGAg                                                                                 <  2:1455423/150‑1 (MQ=255)
                                                                       tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAg                                                                                       >  1:1642683/1‑144 (MQ=255)
                                                                       tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAg                                                                                       <  1:931947/144‑1 (MQ=255)
                                                                       tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAg                                                                                       <  2:1642683/144‑1 (MQ=255)
                                                                       tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAg                                                                                       >  2:931947/1‑144 (MQ=255)
                                                                       tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAg                                                                                 <  2:1245295/150‑1 (MQ=255)
                                                                                 tCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGNGc                                                                                           <  1:1299528/131‑1 (MQ=255)
                                                                                 tCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGc                                                                                           >  2:1299528/1‑131 (MQ=255)
                                                                                 tCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCcgaccgac                                                                       <  2:628551/150‑1 (MQ=255)
                                                                                   aCGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGActc                                                                     >  2:1466307/1‑150 (MQ=255)
                                                                                     gAACGTGTCGAATAAACCGTTTGCGGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACtctt                                                                   >  2:1284545/1‑150 (MQ=255)
                                                                                      aaCGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACtcttc                                                                  >  1:552287/1‑150 (MQ=255)
                                                                                                aTAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAATGCCCGGCAAGCTGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGa                                                       <  1:1731712/150‑1 (MQ=255)
                                                                                                    aCCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTATTCTTCTTCGTTGAccc                                                    >  1:1276108/1‑150 (MQ=255)
                                                                                                       gTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCa                                                 >  1:1733011/1‑150 (MQ=255)
                                                                                                            ccGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGccc                                            <  1:683813/150‑1 (MQ=255)
                                                                                                                  gggAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATcacc                                      <  2:1750094/150‑1 (MQ=255)
                                                                                                                      aaTGGTGAAAAACTGGCACGCAGATGCGAGGCAAGCCGGCTTCTGTTTTAGCTCGTGGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACAt                                  <  2:1452631/150‑1 (MQ=255)
                                                                                                                          cTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAg                              <  1:1799293/150‑1 (MQ=255)
                                                                                                                           tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGc                             >  2:1213592/1‑150 (MQ=255)
                                                                                                                           tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGc                             >  2:378876/1‑150 (MQ=255)
                                                                                                                            gAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCa                            <  1:951726/150‑1 (MQ=255)
                                                                                                                            gAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGACAGCCTTTACAGCGTAGCCACTACAGGCGGTAGCCGGGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCa                            >  1:514559/1‑150 (MQ=255)
                                                                                                                               aaaCTGGCAAAAAAATGCCCGGTAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCt                         <  2:957105/150‑1 (MQ=255)
                                                                                                                                         aaaaaTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGATAATCAc               <  1:1138247/150‑1 (MQ=255)
                                                                                                                                         aaaaaTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCAc               <  2:517552/150‑1 (MQ=255)
                                                                                                                                            aacgCCCGGCATGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTGGCCGCGCTCGTTGTGCCCTGCTGCTGTCGATCGGAGCGACTCCTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTa            <  2:1005293/147‑1 (MQ=255)
                                                                                                                                             ttGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAg           >  1:85922/2‑150 (MQ=255)
                                                                                                                                             ttGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAg           >  1:1262912/2‑150 (MQ=255)
                                                                                                                                             aTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAg           <  1:1213592/150‑1 (MQ=255)
                                                                                                                                              tGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAgg          <  2:85922/150‑1 (MQ=255)
                                                                                                                                                 ccGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGa       >  1:1880336/1‑150 (MQ=255)
                                                                                                                                                   ggCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCt                                                                                               <  2:1886243/61‑1 (MQ=255)
                                                                                                                                                   ggCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCt                                                                                               >  1:1886243/1‑61 (MQ=255)
                                                                                                                                                      aaGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATg  <  2:552287/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CCTTCGGCGCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATG  >  NC_002947/5966238‑5966537

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: