Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,442,672 C→G 7.5% 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,6720CG 7.5% 290.2 / 12.9 95intergenic (‑98/+40)nicP‑II/nicTporin‑like protein/metabolite transport protein NicT
Reads supporting (aligned to +/- strand):  ref base C (42/45);  new base G (5/2);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

CGGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTGA  >  NC_002947/4442526‑4442821
                                                                                                                                                  |                                                                                                                                                     
cgGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCgg                                                                                                                                                    <  2:1599718/150‑1 (MQ=255)
 ggCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGc                                                                                                                                                   >  1:2221983/1‑150 (MQ=255)
 ggCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGc                                                                                                                                                   <  2:1683899/150‑1 (MQ=255)
     aGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGc                                                                                                                                               <  2:172744/150‑1 (MQ=255)
      ggCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcc                                                                                                                                              >  1:283488/1‑150 (MQ=255)
               aCGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGAGGGTTGCATCGAGAGAGGGCGGCAAGCGGGCCCCGGCCCGCctgccctgcc                                                                                                                                     >  2:519776/1‑150 (MQ=255)
                 gaagaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGt                                                                                                                                   >  2:433602/1‑150 (MQ=255)
                  aagaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGtt                                                                                                                                  >  1:777218/1‑150 (MQ=255)
                    gaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGc                                                                                                                                >  1:1663354/1‑150 (MQ=255)
                    gaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGc                                                                                                                                <  1:433602/150‑1 (MQ=255)
                      aaaaCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCCTTGCCACCGGCGTTGGGTTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGt                                                                                                                              <  1:2123570/150‑1 (MQ=255)
                             ggcgCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAgcgc                                                                                                                       >  2:1220632/1‑150 (MQ=255)
                             ggcgCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAgcgc                                                                                                                       <  2:2172849/150‑1 (MQ=255)
                                gcAGTGCACCCTTGCATAGGGCTACTCCTGGATTCGGGCAACGGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCCGCCGGCGGGCCCCGGCCCGTCTGCCCGGTCGTTGCGTGAAGCGCTgg                                                                                                                    <  2:145563/150‑1 (MQ=255)
                                    tGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCg                                                                                                                >  1:445065/1‑150 (MQ=255)
                                        aCCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAg                                                                                                            >  1:1432293/1‑150 (MQ=255)
                                          cTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCg                                                                                                          <  2:2269757/150‑1 (MQ=255)
                                           ttGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCgg                                                                                                         >  2:1417193/1‑150 (MQ=255)
                                               aTAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCTACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGt                                                                                                                              >  1:2142798/1‑125 (MQ=255)
                                               aTAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCTACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGt                                                                                                                              <  2:2142798/125‑1 (MQ=255)
                                               aTAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGa                                                                                                     >  1:2421829/1‑150 (MQ=255)
                                               aTAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGa                                                                                                     >  2:1608306/1‑150 (MQ=255)
                                               aTAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGACACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGa                                                                                                     >  2:1238405/1‑150 (MQ=255)
                                                   ggCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAg                                                                                                 >  2:1734002/1‑150 (MQ=255)
                                                       aCTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTg                                                                                             <  2:2117283/150‑1 (MQ=255)
                                                       aCTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTg                                                                                             <  1:1417193/150‑1 (MQ=255)
                                                          ccTGGATTTGGGCAACAGCAACCCGGCGCGGCCGAATAACATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCTCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGcc                                                                                          <  2:789954/150‑1 (MQ=255)
                                                           cTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCg                                                                                         >  1:1133540/1‑150 (MQ=255)
                                                           cTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCg                                                                                         >  2:290038/1‑150 (MQ=255)
                                                            tGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCgg                                                                                        >  1:1075090/1‑150 (MQ=255)
                                                            tGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGGCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCgg                                                                                        >  1:591097/1‑150 (MQ=255)
                                                               aTTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCaa                                                                                     <  1:214506/150‑1 (MQ=255)
                                                                 ttGGGCAACAGCAACCCGGCGCGGGGGCACGGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAgg                                                                                   <  2:990078/150‑1 (MQ=255)
                                                                     gCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGctgt                                                                               <  2:2221983/150‑4 (MQ=255)
                                                                       aaCAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGa                                                                             >  2:1435126/1‑150 (MQ=255)
                                                                          aGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCa                                                                          >  2:207033/1‑150 (MQ=255)
                                                                           gCAACCCGGCGCGGCGGCATTGCATGGAAACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAg                                                                         <  2:302088/150‑1 (MQ=255)
                                                                           gCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAg                                                                         <  2:1395809/150‑1 (MQ=255)
                                                                            cAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAgg                                                                        >  2:494943/1‑150 (MQ=255)
                                                                            cAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAgg                                                                        <  1:1778616/150‑1 (MQ=255)
                                                                             aaCCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCg                                                                                                          <  2:1631459/115‑1 (MQ=255)
                                                                             aaCCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCg                                                                                                          >  1:1631459/1‑115 (MQ=255)
                                                                                      gggcggcAGTGCATGGCCACCGGCGGGCGGGTTGCATCGAGAGAGGTCGGTAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                              <  2:667053/149‑1 (MQ=255)
                                                                                      cggcggCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                              <  1:103685/150‑1 (MQ=255)
                                                                                      cggcggCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                              <  2:1075090/150‑1 (MQ=255)
                                                                                        gcggcACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTAGGCAGGCGGGCCCCGGCCCGCCTGCCCTGTCGTTGCGTGAAGCGCTGGGCCTCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAgg                                                            >  1:762067/1‑150 (MQ=255)
                                                                                          ggcACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGAt                                                          >  2:2174209/1‑150 (MQ=255)
                                                                                          ggcACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGAt                                                          >  2:2174270/1‑150 (MQ=255)
                                                                                           gcACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATc                                                         <  2:1663354/150‑1 (MQ=255)
                                                                                            cACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGTCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCACGAAGCAGCTTAGGATCa                                                        <  1:519776/150‑1 (MQ=255)
                                                                                            cACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCa                                                        <  1:1245137/150‑1 (MQ=255)
                                                                                             aCTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAg                                                       <  1:2174209/150‑1 (MQ=255)
                                                                                             aCTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAg                                                       <  1:2174270/150‑1 (MQ=255)
                                                                                              cTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTg                                                                                                      >  1:866695/1‑102 (MQ=255)
                                                                                              cTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTg                                                                                                      <  2:866695/102‑1 (MQ=255)
                                                                                              cTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAgg                                                      <  1:207033/150‑1 (MQ=255)
                                                                                              cTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGACGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAgg                                                      <  1:153504/150‑1 (MQ=255)
                                                                                               tGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                              >  2:642455/1‑141 (MQ=255)
                                                                                               tGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                              <  1:642455/141‑1 (MQ=255)
                                                                                                gCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCGGGGCCCGcctgcc                                                                                                                                          >  2:466775/1‑64 (MQ=37)
                                                                                                gCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCGGGGCCCGcctgcc                                                                                                                                          <  1:466775/64‑1 (MQ=37)
                                                                                                 cATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAg                                                   >  2:950016/1‑150 (MQ=255)
                                                                                                  aTGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGc                                                  <  2:1162620/150‑1 (MQ=255)
                                                                                                  aTGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGc                                                  >  1:1811901/1‑150 (MQ=255)
                                                                                                  aTGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGGGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGc                                                  >  2:144045/1‑150 (MQ=255)
                                                                                                   tGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCCGCCTCCCCAAACTCATGGGCCACAAAGCGACAGGCCTCCAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGcc                                                 <  1:168654/150‑1 (MQ=255)
                                                                                                          cGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAac                                          >  1:1442159/1‑150 (MQ=255)
                                                                                                          cGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAac                                          <  2:1007543/150‑1 (MQ=255)
                                                                                                          cGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAac                                          >  1:1716000/1‑150 (MQ=255)
                                                                                                            gCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAacac                                        >  2:2341416/1‑150 (MQ=255)
                                                                                                             cGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACAcc                                       <  2:445065/150‑1 (MQ=255)
                                                                                                                  gggTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcag                                  <  1:1734002/150‑1 (MQ=255)
                                                                                                                      tGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgca                              <  2:777218/150‑1 (MQ=255)
                                                                                                                        cATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagc                            <  1:1238405/150‑1 (MQ=255)
                                                                                                                         aTCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCGAGGAAGGATGTGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagcg                           >  2:918815/1‑150 (MQ=255)
                                                                                                                         aTCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagcg                           >  1:631473/1‑150 (MQ=255)
                                                                                                                         aTCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagcg                           >  1:621981/1‑150 (MQ=255)
                                                                                                                           cGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTg                         >  2:1579987/1‑150 (MQ=255)
                                                                                                                             agagagGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGc                                                                                                               <  2:998753/62‑1 (MQ=255)
                                                                                                                             agagagGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGc                                                                                                               >  1:998753/1‑62 (MQ=255)
                                                                                                                              gagagGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGa                      <  2:621981/150‑1 (MQ=255)
                                                                                                                                 agGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGt                   <  2:631473/150‑1 (MQ=255)
                                                                                                                                   gTCGGAAGGCGGGCCCCGGCCCGCGTCCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTgg                 <  2:2421829/150‑2 (MQ=255)
                                                                                                                                    tCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCt                                                                                               <  1:1593270/71‑1 (MQ=255)
                                                                                                                                    tCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCt                                                                                               >  2:1593270/1‑71 (MQ=255)
                                                                                                                                    tCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGcc                >  1:841863/1‑150 (MQ=255)
                                                                                                                                      ggCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCgg              <  1:1579987/150‑1 (MQ=255)
                                                                                                                                        cAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCNGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGNTg            <  1:1220632/150‑1 (MQ=255)
                                                                                                                                         gggCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa           <  2:283488/149‑1 (MQ=255)
                                                                                                                                         aGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGCTCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa           >  2:105528/1‑150 (MQ=255)
                                                                                                                                         aGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATTAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa           >  2:2012265/1‑150 (MQ=255)
                                                                                                                                         aGGCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa           >  1:1621379/1‑150 (MQ=255)
                                                                                                                                         aGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa           >  1:1410611/1‑150 (MQ=255)
                                                                                                                                          ggCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAgg                                                                                   >  2:648714/1‑77 (MQ=39)
                                                                                                                                          ggCGGGCCGGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAgg                                                                                   <  1:648714/77‑1 (MQ=39)
                                                                                                                                                ccccGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGt    <  1:2341416/150‑1 (MQ=255)
                                                                                                                                                  ccGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTGa  >  2:1487328/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTGA  >  NC_002947/4442526‑4442821

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: