Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,442,628 C→G 6.9% intergenic (‑98/+40) 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,6280CG 6.9% 222.0 / 6.4 73intergenic (‑98/+40)nicP‑II/nicTporin‑like protein/metabolite transport protein NicT
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base G (2/3);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

CGGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCG  >  NC_002947/4442482‑4442773
                                                                                                                                                     |                                                                                                                                                  
cgGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCg                                                                                                                                                  >  2:120949/1‑149 (MQ=255)
 ggCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCgg                                                                                                                                                 >  3:280305/1‑149 (MQ=255)
  gCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGc                                                                                                                                                >  3:68795/1‑149 (MQ=255)
      tgtgTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGccc                                                                                                                                              <  3:276994/144‑1 (MQ=255)
        cGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCGGGGCCCGcct                                                                                                                                          <  8:48414/149‑1 (MQ=255)
        cGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCGGGGCCCGcct                                                                                                                                          <  4:141694/149‑1 (MQ=255)
          tGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTGCTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGcctgc                                                                                                                                        >  2:322308/1‑149 (MQ=255)
             aGACGAAGAAAAACGCGGCGCAGTGAAACCATGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGcctgccct                                                                                                                                     <  3:264750/149‑1 (MQ=255)
               aCGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGcctgccctgc                                                                                                                                   >  8:172221/1‑149 (MQ=255)
                cGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCCGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCctgccctgcc                                                                                                                                  <  1:120949/149‑1 (MQ=255)
                cGAACAAAAACGCCGCGCAGTGAAACCTTCCATACGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCCCGGCGGCACTGCATGGCCACCGGCGTGCTGGTTGCATCTACAGACG‑‑TCGGCAGGCGGGCCCCGGCCCGCctgccctgcc                                                                                                                                  <  5:360158/149‑1 (MQ=255)
                 gaagaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCg                                                                                                                                 >  8:243326/1‑149 (MQ=255)
                 gaagaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCg                                                                                                                                 >  3:174494/1‑149 (MQ=255)
                     aaaaaCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGc                                                                                                                             >  7:334237/1‑149 (MQ=255)
                       aaaCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGt                                                                                                                           <  2:29377/149‑1 (MQ=255)
                         aCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCGGGGCCCGc                                                                                                                                            >  2:45570/1‑130 (MQ=255)
                         aCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCGGGGCCCGc                                                                                                                                            <  1:45570/130‑1 (MQ=255)
                          cgcgGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGaa                                                                                                                        <  2:74398/149‑1 (MQ=255)
                                   gTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGc                                                                                                               >  4:47927/1‑149 (MQ=255)
                                    tGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCCGGCGGGCCCCCGCCCGCCTGGCCTGCCGGTTCGTGAAGCGCTGGGcc                                                                                                              >  1:360260/1‑149 (MQ=255)
                                    tGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGcc                                                                                                              >  5:3414/1‑149 (MQ=255)
                                        aCCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGc                                                                                                            >  6:345224/1‑147 (MQ=255)
                                         ccTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCAGCCCAGCCGTTGCGTGAAGTGCTGGGCCGCTAg                                                                                                         >  8:138313/1‑149 (MQ=255)
                                          cTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGc                                                                                                        >  7:242891/1‑149 (MQ=255)
                                              cATAGGGCTGCTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGtt                                                                                                    <  1:322308/149‑1 (MQ=255)
                                                tAGGGCTACTCCTTGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGa                                                                                                  <  4:280305/149‑1 (MQ=255)
                                                  gggCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGAc                                                                                                 <  3:115445/148‑1 (MQ=255)
                                                     caaCTCCTCCATTTGCGCAACAGCAACCCGGCGCGTCCGCACTCCATGGCCACCGGCCTGCGGGTAGCATCGAGAGAGG‑‑TCGGCAGCCCGGCCCCGGCCCCCCTGCCCAGCCGTTGCGTCAAGCGCAGGGCCGCAAGCGGTTGACCAGc                                                                                             >  2:67651/3‑149 (MQ=255)
                                                         tCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGAcc                                                                                                >  3:212975/1‑142 (MQ=255)
                                                         tCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGAcc                                                                                                <  4:212975/142‑1 (MQ=255)
                                                         tCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTgg                                                                                        >  6:307194/1‑149 (MQ=255)
                                                          ccTGGATTTCGCCCACAGCAACCCGGCGCGGCGGCACGGCCTGGCCACCGGCGGGCGGGTTGCATCGAGCGAGG‑‑CCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTCCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGc                                                                                       <  4:68795/149‑1 (MQ=255)
                                                           ttgtaTTGGGGAAACCGCACCCGGGCGGGGGGCCACTGCTGGCCCCCGGGGGTGGGGTTTGCATGTAGTCTGG‑‑GCGGCTGGGGGGCCCGGGCCGGCGTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGcc                                                                                      <  2:360260/145‑1 (MQ=255)
                                                           cTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCCGCCCCCCTGCCCAGCCCTTGCGTGAAGCGCTGGGGCGCGAGCGGTTGACCAGCGTgggc                                                                                      >  1:360476/1‑147 (MQ=255)
                                                            tGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTg                                                                                                                          >  7:63596/1‑113 (MQ=255)
                                                            tGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTg                                                                                                                          <  8:63596/113‑1 (MQ=255)
                                                                  tGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAg                                                                                <  7:172221/148‑1 (MQ=255)
                                                                     gCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGa                                                                            >  8:281001/1‑149 (MQ=255)
                                                                      cAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                           <  4:174494/149‑1 (MQ=255)
                                                                      cAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                           >  8:117036/1‑149 (MQ=255)
                                                                       aaCAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                           >  3:137369/1‑148 (MQ=255)
                                                                         cAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGAc                                                                        >  1:224944/1‑149 (MQ=255)
                                                                          agaaaCCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGAAGAGGCAGGCGCGCCCCGGCCCGCCTGCCATGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGa                                                                         >  1:107971/4‑149 (MQ=255)
                                                                              aCCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGa                                                                   <  6:3414/149‑1 (MQ=255)
                                                                               cccGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGaa                                                                  >  8:167969/1‑149 (MQ=255)
                                                                                      cggcggCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTg                                                           <  5:307194/149‑1 (MQ=255)
                                                                                       ggcggcACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                          <  3:47927/149‑1 (MQ=255)
                                                                                       ggcggcACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                          <  2:224944/149‑1 (MQ=255)
                                                                                        gcggcACTGATGGGCCCCGGGGGTGCGGTTTGAATCGAGTTCGG‑‑GCGGCTGGCGGGCCCCGGCCGGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAg                                                         <  2:360476/149‑1 (MQ=255)
                                                                                         cggcACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAgg                                                        <  7:117036/149‑1 (MQ=255)
                                                                                            cACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATc                                                     <  5:345224/149‑1 (MQ=255)
                                                                                                 cATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAgg                                                  <  7:281001/147‑1 (MQ=255)
                                                                                                 cATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGc                                                 <  7:138313/148‑1 (MQ=255)
                                                                                                  aTGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAg                                               >  6:24794/1‑149 (MQ=255)
                                                                                                             gCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAaca                                     >  5:332099/1‑149 (MQ=255)
                                                                                                             gCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCAGGCAAaca                                     >  7:353069/1‑149 (MQ=255)
                                                                                                                          aTCGTGAGTGG‑‑GCGGCAGACGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCTCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcag                              <  4:60519/143‑1 (MQ=255)
                                                                                                                          aTCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagc                        <  7:167969/149‑1 (MQ=255)
                                                                                                                          aTCGAGAGAGG‑‑TCGGCAGACGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcag                              >  3:60519/1‑143 (MQ=255)
                                                                                                                           tCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGCCCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagcg                       >  2:87128/1‑149 (MQ=255)
                                                                                                                                agagG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGc                    >  1:205492/1‑147 (MQ=255)
                                                                                                                                 gagG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAg                 <  7:243326/149‑1 (MQ=255)
                                                                                                                                    g‑‑TCGGCCGGCGGGCCCCGGCCCGCCGGCCCGGCCGGTGCGGGAAGCGCGGGGGCGCGAGCGGGTGGCCAGCGTGGGCGGCCAGGCGATGACCAGGAAGCAGCGGAGGATCAGGCAGCCGGCCAACCCCCACAGCCCCGCGTGCGAGGTg              >  8:20854/1‑149 (MQ=255)
                                                                                                                                    g‑‑TCGGCCGGCGGGCCCCCGCCCCCCGGCCCTGCCGTTGCGTGAAGCGCGGGGCCGCGAGCGGGTGGACAGCATGGGCCGGCAGGCGGGGGCCCCGGAGCAGCGGAGGGTGAGGCCGCCCGCCAACCCCCACAGCCCAGCGTGCggcgtg              >  8:221679/1‑144 (MQ=255)
                                                                                                                                        cGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGcc            >  2:363024/1‑149 (MQ=255)
                                                                                                                                         ggCAGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCg           >  3:154667/1‑149 (MQ=255)
                                                                                                                                         ggCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCg           <  7:221679/149‑1 (MQ=255)
                                                                                                                                          gCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCgg          >  5:268576/1‑149 (MQ=255)
                                                                                                                                          gCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCgg          >  3:331788/1‑149 (MQ=255)
                                                                                                                                            aGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTg        >  7:200315/1‑149 (MQ=255)
                                                                                                                                              gCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAg      <  8:200315/149‑1 (MQ=255)
                                                                                                                                                gggCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGt    <  4:331788/149‑1 (MQ=255)
                                                                                                                                                  gCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCg  <  4:137369/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGG‑‑TCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCG  >  NC_002947/4442482‑4442773

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