Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,966,389 A→T 8.0% intergenic (‑402/+28) 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,3890AT 8.0% 210.9 / 5.6 71intergenic (‑402/+28)PP_5232/amtBhypothetical protein/ammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base A (25/40);  new base T (2/4);  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.97e-01

TCGGCGCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATG  >  NC_002947/5966241‑5966537
                                                                                                                                                    |                                                                                                                                                     
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:1293691/1‑150 (MQ=255)
                     ccGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGt                                                                                                                                 >  1:628551/1‑150 (MQ=255)
                        tGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGc                                                                                                                              >  1:1245295/1‑150 (MQ=255)
                        tGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGc                                                                                                                              <  2:655459/150‑1 (MQ=255)
                        tGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGc                                                                                                                              <  1:451205/150‑1 (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)
                                                             tAATTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑Ca                                                                                        <  1:1433511/150‑1 (MQ=255)
                                                             tAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑Ca                                                                                        >  1:203871/1‑150 (MQ=255)
                                                                gTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTCCCGCGTAGTACTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGt                                                                                     >  1:1452631/1‑150 (MQ=255)
                                                                 ttttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTc                                                                                    >  1:1054359/1‑150 (MQ=255)
                                                                   ttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGa                                                                                  <  1:723520/150‑1 (MQ=255)
                                                                    tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGGGCGCGTTGTGCTGAGC‑CAGGTCGAg                                                                                 <  2:1455423/150‑1 (MQ=255)
                                                                    tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAg                                                                                       <  1:931947/144‑1 (MQ=255)
                                                                    tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAg                                                                                       >  2:931947/1‑144 (MQ=255)
                                                                    tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAg                                                                                       <  2:1642683/144‑1 (MQ=255)
                                                                    tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAg                                                                                       >  1:1642683/1‑144 (MQ=255)
                                                                    tttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAg                                                                                 <  2:1245295/150‑1 (MQ=255)
                                                                              tCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGNGc                                                                                           <  1:1299528/131‑1 (MQ=255)
                                                                              tCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGc                                                                                           >  2:1299528/1‑131 (MQ=255)
                                                                              tCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCcgaccgac                                                                       <  2:628551/150‑1 (MQ=255)
                                                                                aCGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGActc                                                                     >  2:1466307/1‑150 (MQ=255)
                                                                                  gAACGTGTCGAATAAACCGTTTGCGGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACtctt                                                                   >  2:1284545/1‑150 (MQ=255)
                                                                                   aaCGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACtcttc                                                                  >  1:552287/1‑150 (MQ=255)
                                                                                             aTAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGC‑AAAAAATGCCCGGCAAGCTGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGa                                                       <  1:1731712/150‑1 (MQ=255)
                                                                                                 aCCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTATTCTTCTTCGTTGAccc                                                    >  1:1276108/1‑150 (MQ=255)
                                                                                                    gTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCa                                                 >  1:1733011/1‑150 (MQ=255)
                                                                                                         ccGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGccc                                            <  1:683813/150‑1 (MQ=255)
                                                                                                               gggAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATcacc                                      <  2:1750094/150‑1 (MQ=255)
                                                                                                                   aaTGGTGAAAAACTGGCACGCAGATGCGAGGCAAGCCGGCTTCTGTTTTAGCTCGTGGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACAt                                  <  2:1452631/150‑1 (MQ=255)
                                                                                                                       cTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAg                              <  1:1799293/150‑1 (MQ=255)
                                                                                                                        tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGc                             >  2:1213592/1‑150 (MQ=255)
                                                                                                                        tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGc                             >  2:378876/1‑150 (MQ=255)
                                                                                                                         gAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCa                            <  1:951726/150‑1 (MQ=255)
                                                                                                                         gAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGACAGCCTTTACAGCGTAGCCACTACAGGCGGTAGCCGGGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCa                            >  1:514559/1‑150 (MQ=255)
                                                                                                                            aaaCTGGCAAAAAAATGCCCGGTAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCt                         <  2:957105/150‑1 (MQ=255)
                                                                                                                                      aaaaaTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGATAATCAc               <  1:1138247/150‑1 (MQ=255)
                                                                                                                                      aaaaaTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCAc               <  2:517552/150‑1 (MQ=255)
                                                                                                                                         aacgCCCGGCATGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTGGCCGCGCTCGTTGTGCCCTGCTGCTGTCGA‑TCGGAGCGACTCCTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTa            <  2:1005293/147‑1 (MQ=255)
                                                                                                                                          ttGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAg           >  1:85922/2‑150 (MQ=255)
                                                                                                                                          ttGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAg           >  1:1262912/2‑150 (MQ=255)
                                                                                                                                          aTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAg           <  1:1213592/150‑1 (MQ=255)
                                                                                                                                           tGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAgg          <  2:85922/150‑1 (MQ=255)
                                                                                                                                              ccGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGa       >  1:1880336/1‑150 (MQ=255)
                                                                                                                                                ggCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCt                                                                                               <  2:1886243/61‑1 (MQ=255)
                                                                                                                                                ggCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCt                                                                                               >  1:1886243/1‑61 (MQ=255)
                                                                                                                                                   aaGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATg  <  2:552287/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCGGCGCGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGC‑CAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATG  >  NC_002947/5966241‑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: ‑