Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,442,673 C→G 5.4% intergenic (‑99/+39) nicP‑II ← / ← nicT porin‑like protein/metabolite transport protein NicT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,442,6730CG 5.4% 238.2 / 4.9 74intergenic (‑99/+39)nicP‑II/nicTporin‑like protein/metabolite transport protein NicT
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base G (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

CCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTGAC  >  NC_002947/4442529‑4442822
                                                                                                                                                |                                                                                                                                                     
ccAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGccc                                                                                                                                                  <  1:1531991/150‑1 (MQ=255)
    gCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCGGGGCCCGcct                                                                                                                                              <  2:165539/150‑1 (MQ=255)
     cGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcctg                                                                                                                                             <  2:1798321/150‑1 (MQ=255)
     cGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGGGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcctg                                                                                                                                             <  1:105754/150‑1 (MQ=255)
           gACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGGTTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcctgccctgc                                                                                                                                       >  1:1594232/1‑150 (MQ=255)
            aCGAAGAAAAAGGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCctgccctgcc                                                                                                                                      <  2:1562861/150‑1 (MQ=255)
              gaagaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGc                                                                                                                                                >  2:1185280/1‑138 (MQ=255)
              gaagaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGc                                                                                                                                                <  1:1185280/138‑1 (MQ=255)
              gaagaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGt                                                                                                                                    >  1:1642383/1‑150 (MQ=255)
               aagaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGtt                                                                                                                                   >  1:1096175/1‑150 (MQ=255)
                 gaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGc                                                                                                                                 <  2:1642383/150‑1 (MQ=255)
                 aaaaaaCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGc                                                                                                                                 <  2:1096175/149‑1 (MQ=255)
                                gTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGcc                                                                                                                  >  2:1835026/1‑150 (MQ=255)
                                     aCCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAg                                                                                                             >  2:1039644/1‑150 (MQ=255)
                                      ccTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcct                                                                                                                                              <  1:1730132/116‑1 (MQ=255)
                                      ccTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcct                                                                                                                                              >  2:1730132/1‑116 (MQ=255)
                                      ccTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCt                                                                                                                       >  1:624279/1‑139 (MQ=255)
                                      ccTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCt                                                                                                                       <  2:624279/139‑1 (MQ=255)
                                      ccTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGc                                                                                                            >  1:390159/1‑150 (MQ=255)
                                        ttGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCgg                                                                                                          <  2:1695132/150‑1 (MQ=255)
                                          gCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCGCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGtc                                                                                                        <  1:570516/150‑2 (MQ=255)
                                            aTAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCTAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGa                                                                                                      <  1:233266/150‑1 (MQ=255)
                                              aGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGAcc                                                                                                    >  2:1704344/1‑150 (MQ=255)
                                               gggCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCa                                                                                                   >  2:183759/1‑150 (MQ=255)
                                                 cgcaCTCCTGGATATGGGCAACAGCAACCAAGCGCGGCGGCACTGCATTGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGc                                                                                                 >  2:1530091/4‑150 (MQ=255)
                                                    aCTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGcc                                                                                                                                                   <  1:890316/97‑1 (MQ=255)
                                                    aCTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGcc                                                                                                                                                   >  2:890316/1‑97 (MQ=255)
                                                       ccTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCctgccctgcc                                                                                                                                      >  2:1671068/1‑107 (MQ=255)
                                                       ccTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCctgccctgcc                                                                                                                                      <  1:1671068/107‑1 (MQ=255)
                                                        cTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGcccc                                                                                                                                                       <  2:674689/89‑1 (MQ=255)
                                                        cTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGcccc                                                                                                                                                       >  1:674689/1‑89 (MQ=255)
                                                        cTGGATTTGGGCAACAGCAACCCGGCGCGGCGACACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCg                                                                                          >  2:671481/1‑150 (MQ=255)
                                                              tttggCAACAGCCACCCGGCGCGGCGGCCCCGCATGGCGACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAgg                                                                                    <  2:317392/147‑1 (MQ=255)
                                                               tGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGc                                                                                   <  1:1316794/150‑1 (MQ=255)
                                                                  gCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                                >  1:265301/1‑150 (MQ=255)
                                                                  gCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                                >  2:1098221/1‑150 (MQ=255)
                                                                          aaCCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCAGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGa                                                                        <  1:1096519/150‑1 (MQ=255)
                                                                           aCCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGaa                                                                       <  2:1657294/150‑1 (MQ=255)
                                                                                gcgcggcgGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAagcagc                                                                  >  2:44913/1‑150 (MQ=255)
                                                                                          aCTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAg                                                        >  2:1369444/1‑150 (MQ=255)
                                                                                              cATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAg                                                    >  2:35692/1‑150 (MQ=255)
                                                                                               aTGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGc                                                   >  2:1579899/1‑150 (MQ=255)
                                                                                                 ggCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCg                                                 >  1:239308/1‑150 (MQ=255)
                                                                                                 ggCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCg                                                 <  2:1368961/150‑1 (MQ=255)
                                                                                                  gCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCgg                                                <  1:1369444/150‑1 (MQ=255)
                                                                                                   ccACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGc                                               <  2:169056/150‑1 (MQ=255)
                                                                                                    cACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCa                                              <  1:1579899/150‑1 (MQ=255)
                                                                                                     aCCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                                <  1:1436899/115‑1 (MQ=255)
                                                                                                     aCCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                                >  2:1436899/1‑115 (MQ=255)
                                                                                                     aCCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCaa                                             >  2:148007/1‑150 (MQ=255)
                                                                                                                 gTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAagcag                                                                   >  2:290257/1‑116 (MQ=255)
                                                                                                                 gTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAagcag                                                                   <  1:290257/116‑1 (MQ=255)
                                                                                                                  ttGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGc                                                                                   >  1:106598/1‑99 (MQ=255)
                                                                                                                  ttGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGc                                                                                   <  2:106598/99‑1 (MQ=255)
                                                                                                                   tGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgca                               <  2:1594232/150‑1 (MQ=255)
                                                                                                                     cATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagc                             <  1:183759/150‑1 (MQ=255)
                                                                                                                      aTCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                                >  1:616673/1‑98 (MQ=255)
                                                                                                                      aTCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                                <  2:616673/98‑1 (MQ=255)
                                                                                                                         gagagaGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAg                                                    >  2:303618/1‑123 (MQ=255)
                                                                                                                         gagagaGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAg                                                    <  1:303618/123‑1 (MQ=255)
                                                                                                                          agagagGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGAAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCg                        >  1:1340055/1‑150 (MQ=255)
                                                                                                                                gTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACAAGGAagc                                                                     >  1:535902/1‑99 (MQ=255)
                                                                                                                                gTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACAAGGAagc                                                                     <  2:535902/99‑1 (MQ=255)
                                                                                                                                   ggCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCgg               >  1:1176461/1‑150 (MQ=255)
                                                                                                                                    gCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGt              <  1:1835026/150‑1 (MQ=255)
                                                                                                                                     cAGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTg             >  2:1898296/1‑150 (MQ=255)
                                                                                                                                      aGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGc                         >  2:1840227/1‑137 (MQ=255)
                                                                                                                                      aGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGc                         <  1:1840227/137‑1 (MQ=255)
                                                                                                                                      aGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa            <  1:1039644/150‑1 (MQ=255)
                                                                                                                                       ggCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAg           <  2:1176461/150‑1 (MQ=255)
                                                                                                                                       ggCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAg           <  1:1084402/150‑1 (MQ=255)
                                                                                                                                        gCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAgg          >  1:1061275/1‑150 (MQ=255)
                                                                                                                                          gggCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTc        >  2:351199/1‑150 (MQ=255)
                                                                                                                                          gggCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTc        <  1:44913/150‑1 (MQ=255)
                                                                                                                                           ggCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCg       <  2:1061275/150‑1 (MQ=255)
                                                                                                                                           ggCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCg       >  1:969279/1‑150 (MQ=255)
                                                                                                                                             ccccGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGCTGACACGGAAGTAGCTGAGAATCAGTTAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGt     <  1:35692/150‑1 (MQ=255)
                                                                                                                                             ccccGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGAAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGt     <  2:1340055/150‑1 (MQ=255)
                                                                                                                                             ccccGGCCCGCCTGCCCTGCCGTTGCGGGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAGGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGACGGt     >  2:1519422/1‑150 (MQ=255)
                                                                                                                                              cccGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTg    >  2:1937188/1‑150 (MQ=255)
                                                                                                                                              cccGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTg    <  1:85750/150‑1 (MQ=255)
                                                                                                                                               ccGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                                <  2:102785/73‑1 (MQ=255)
                                                                                                                                               ccGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                                >  1:102785/1‑73 (MQ=255)
                                                                                                                                               ccGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTGa   <  2:390159/150‑1 (MQ=255)
                                                                                                                                                cGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTGAc  >  1:1890324/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTGAC  >  NC_002947/4442529‑4442822

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: