Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,585 A→G 5.2% L17P (CTC→CCC)  flhA ← flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,936,5850AG 5.2% 264.1 / 3.5 95L17P (CTC→CCC) flhAflagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base A (45/45);  new base G (3/2);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

AAGGCGTATACACAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATG  >  NC_002947/4936438‑4936734
                                                                                                                                                   |                                                                                                                                                      
aaGGCGTATACACAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCTGGCCGGGc                                                                                                                                                      <  1:1835464/150‑1 (MQ=255)
 aGGCGTATACACAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGcc                                                                                                                                                     <  1:1326121/150‑1 (MQ=38)
         acacAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAgg                                                                                                                                             >  1:223624/1‑150 (MQ=39)
           acaGACCAGCAGCACCACGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGgtt                                                                                                                                           >  1:1314138/1‑150 (MQ=39)
           acaGACCAGCAGCACCACGATCGACAGGGCGATGCTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGgtt                                                                                                                                           >  1:747030/1‑150 (MQ=255)
                ccAGCAGCACCACGATCGACAGGTCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGCTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCCCTAACAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGc                                                                                                                                      <  1:1124186/150‑1 (MQ=11)
                 cagcagCACCACGATAGCCAGGGCGCTGATGACCGAGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCg                                                                                                                                     <  2:593622/150‑1 (MQ=11)
                  agcagcCCCACGATCGACAGGGCGATGGTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCACCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCgg                                                                                                                                    <  2:643455/150‑1 (MQ=255)
                  agcagcACCACGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCgg                                                                                                                                    <  2:223624/150‑1 (MQ=39)
                   gcagcaCCACGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCggg                                                                                                                                   <  1:444216/150‑1 (MQ=38)
                    cagcaCCACGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGc                                                                                                                                  >  2:801001/1‑150 (MQ=39)
                          cacGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGGACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCt                                                                                                                            >  1:383859/1‑150 (MQ=38)
                          cacGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCt                                                                                                                            >  2:1904373/1‑150 (MQ=39)
                          cacGATCGACAGGGAGATCTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCt                                                                                                                            <  2:1556741/150‑1 (MQ=255)
                                cGACAGGGCGATGTTGAACGTGAAGAACACGTCGAACAGGAACGGCGGTATCGGCAACATCATCATCGCCAACATCACCAGCAGCAACAGCGGCACCCACAGGTTGCGCTGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTaa                                                                                                                      <  2:526579/150‑1 (MQ=11)
                                  aCAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACt                                                                                                                    >  2:1752732/1‑150 (MQ=39)
                                      ggCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGt                                                                                                                >  1:393575/1‑150 (MQ=39)
                                         gATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCg                                                                                                             >  1:1111804/1‑150 (MQ=38)
                                           tGTTGAACGTGACGAACCCGTCGCGAAGGAACGGCGGTATCGGCAACATCATCATGTCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGAt                                                                                                           >  2:512998/1‑150 (MQ=11)
                                                   gTGAAGAACACGTCGAGCCGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCCGCCGCAACCGCAGCACCCCCAGGTCGCCCCGGCCGAGCCCCACCAGGTTGTTGCGGGCGTTGCAGATAAACAGAATGCGATCCACcgcg                                                                                                   >  2:599639/1‑150 (MQ=11)
                                                   gTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCggg                                                                                                                                   >  2:438479/1‑118 (MQ=38)
                                                   gTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCggg                                                                                                                                   <  1:438479/118‑1 (MQ=38)
                                                   gTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCggg                                                                                                                                   <  1:1403191/118‑1 (MQ=38)
                                                   gTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCggg                                                                                                                                   >  2:1403191/1‑118 (MQ=38)
                                                   gTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCAAGGTTGCTGTGGGTGTCGCTGAGGAACGTAGCGCGATCCACcgcg                                                                                                   >  2:1666823/1‑150 (MQ=11)
                                                    tGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGa                                                                                                  <  1:1752732/150‑1 (MQ=35)
                                                    tGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGa                                                                                                  >  2:690826/1‑150 (MQ=34)
                                                    tGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAAAATCATAATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGa                                                                                                  >  2:710703/1‑150 (MQ=21)
                                                    accaggacacGTCGAGCAGGAACGGCGGTAGCGGAAACCTCATCATTGCAAACATCCCCAGGAGCATCAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTCGTCGCGGGAGTTGCTGATTAACTGAGTGCCATCCACCGCGa                                                                                                  <  2:932805/144‑1 (MQ=11)
                                                     gaagaaCACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAg                                                                                                 >  1:159682/1‑150 (MQ=35)
                                                      aagaaCACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGt                                                                                                >  1:478532/1‑150 (MQ=35)
                                                       agaaCACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTc                                                                                               <  2:393575/150‑1 (MQ=35)
                                                           cacGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGAt                                                                                                                         <  1:1441378/120‑1 (MQ=38)
                                                           cacGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACAAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGAt                                                                                                                         >  2:1441378/1‑120 (MQ=37)
                                                               tCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGAt                                                                                       <  2:1895849/150‑1 (MQ=35)
                                                                  aGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCa                                                                                    >  1:424201/1‑150 (MQ=35)
                                                                      ggAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAAc                                                                                >  1:591297/1‑150 (MQ=35)
                                                                                  tCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAg                                                                   <  1:237592/150‑1 (MQ=35)
                                                                                     gCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAgggg                                                                >  1:1805871/1‑150 (MQ=35)
                                                                                           tcatcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGa                                                                           <  1:472140/134‑1 (MQ=37)
                                                                                           tcatcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGa                                                                           >  2:472140/1‑134 (MQ=37)
                                                                                           tcatcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCt                                                          <  2:1314138/150‑1 (MQ=35)
                                                                                           tcatcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAACTTTGACGCCCAGGGGGCGCCTg                                                         >  2:1438083/1‑150 (MQ=25)
                                                                                            catcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGa                                                                                                                          >  2:1064897/1‑86 (MQ=40)
                                                                                            catcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGa                                                                                                                          <  1:1064897/86‑1 (MQ=40)
                                                                                            catcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTg                                                         >  1:742864/1‑150 (MQ=35)
                                                                                             atcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTgg                                                        >  1:1196589/1‑150 (MQ=35)
                                                                                              tcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCTGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTggg                                                       <  1:1666823/150‑1 (MQ=38)
                                                                                              tcatTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTggg                                                       >  1:1134531/1‑150 (MQ=38)
                                                                                                atTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCg                                                     >  1:850316/1‑150 (MQ=39)
                                                                                                 ttGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGt                                                    <  1:690826/150‑1 (MQ=38)
                                                                                                  tGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGCGCCCCCTGGGCGtc                                                   <  1:27637/150‑2 (MQ=255)
                                                                                                   gCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCACGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGGGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTg                                                   <  2:1292542/149‑1 (MQ=255)
                                                                                                   gCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCACGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTg                                                   >  1:1292542/1‑149 (MQ=255)
                                                                                                    ccAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGccc                                                                      <  2:1616712/130‑1 (MQ=37)
                                                                                                    ccAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGccc                                                                      >  1:1616712/1‑130 (MQ=37)
                                                                                                    ccAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGc                                                 >  1:1231613/1‑150 (MQ=255)
                                                                                                     cAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGtcaaa                                                <  1:710703/150‑5 (MQ=255)
                                                                                                         aTCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGc                                            <  2:1231613/150‑1 (MQ=255)
                                                                                                           cACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaa                                          >  2:373537/1‑150 (MQ=255)
                                                                                                           cACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGCGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaa                                          >  2:1005447/1‑150 (MQ=255)
                                                                                                            aCCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGGCGCCTg                                                         >  1:848370/1‑135 (MQ=255)
                                                                                                            aCCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGGCGCCTg                                                         <  2:848370/135‑1 (MQ=255)
                                                                                                             ccAGCAGCAACAGCGGCACCCCCGGGTTGCCCCGGCCGAGCCCGTCGAGGTTGTTGCGGGCGTTGCTGGTTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTATGACGCCCAGGAGGCGCCGATGCATCGCTTTGCcaga                                        >  2:316833/1‑150 (MQ=11)
                                                                                                              cagcagCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCaagaa                                       <  2:998187/150‑1 (MQ=255)
                                                                                                                  agcaACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAACCTTTGACGCCCAggg                                                                 <  2:382300/120‑1 (MQ=37)
                                                                                                                  agcaACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAggg                                                                 >  1:382300/1‑120 (MQ=37)
                                                                                                                  agcaACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCtt                                   <  1:512998/150‑1 (MQ=255)
                                                                                                                  agcaACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCttt                                  <  1:1438083/150‑1 (MQ=255)
                                                                                                                   gcaACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCttt                                  <  2:747030/150‑1 (MQ=255)
                                                                                                                    caACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTc                                 <  1:656117/150‑1 (MQ=255)
                                                                                                                     aaCAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGGTTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTcc                                >  2:299904/1‑150 (MQ=255)
                                                                                                                       cAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCaa                              <  2:850316/150‑1 (MQ=255)
                                                                                                                        aGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAAc                             <  1:801001/150‑1 (MQ=255)
                                                                                                                            gCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTg                         <  1:373537/150‑1 (MQ=255)
                                                                                                                               cccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCa                      <  2:1805871/150‑1 (MQ=255)
                                                                                                                               cccccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCa                      <  2:1134531/150‑1 (MQ=255)
                                                                                                                                 cccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGa                    >  2:1648550/1‑150 (MQ=255)
                                                                                                                                 cccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCACAGGCTGCCCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAgc                    <  1:1899464/150‑2 (MQ=255)
                                                                                                                                  ccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGac                   >  2:343460/1‑150 (MQ=255)
                                                                                                                                   cAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                        <  2:1940168/144‑1 (MQ=255)
                                                                                                                                   cAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTgg                        >  1:1940168/1‑144 (MQ=255)
                                                                                                                                    aGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGacac                 <  2:591297/150‑1 (MQ=255)
                                                                                                                                       ttGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCca              <  1:1648550/150‑1 (MQ=255)
                                                                                                                                         aacctggCCGGGCTCGGCCGGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGc                                                      <  2:509689/103‑1 (MQ=11)
                                                                                                                                         aacctggCCGGGCTCGGCCGGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGc                                                      >  1:509689/6‑108 (MQ=11)
                                                                                                                                          ccccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTACAACTTTGGCAGACACCCacac           >  1:1067650/1‑150 (MQ=255)
                                                                                                                                           cctggCCGGGCTCGGCCGGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACg          >  1:1452490/4‑150 (MQ=255)
                                                                                                                                           cctggCCGGGCTCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACg          >  1:1447928/4‑150 (MQ=255)
                                                                                                                                           cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCacac           <  2:1301745/149‑1 (MQ=255)
                                                                                                                                           cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCacac           >  1:1301745/1‑149 (MQ=255)
                                                                                                                                           cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACg          >  1:365732/1‑150 (MQ=255)
                                                                                                                                               gCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGaaaa      <  2:742864/150‑1 (MQ=255)
                                                                                                                                                 cGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGa    <  1:343460/150‑1 (MQ=255)
                                                                                                                                                 cGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGa    >  1:1782454/1‑150 (MQ=255)
                                                                                                                                                   aGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATg  <  2:383859/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
AAGGCGTATACACAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAACGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGCAACAGCGGCACCCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGGGCGTGGCTTTGCAAGAAGCTTTCCAACTTTGGCAGACACCCACACGAAAAGATG  >  NC_002947/4936438‑4936734

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: