Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,988,018 C→G 7.0% P97R (CCG→CGG)  flgA → flagella basal body P‑ring formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,988,0180CG 7.0% 654.6 / 27.8 198P97R (CCG→CGG) flgAflagella basal body P‑ring formation protein
Reads supporting (aligned to +/- strand):  ref base C (113/71);  new base G (7/7);  total (120/78)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

TTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGCCC  >  NC_002947/4987873‑4988167
                                                                                                                                                 |                                                                                                                                                     
tttCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACAGGAAAGCCCGGCCCAGCGGCTg                                                                                                                                                   >  1:674085/1‑150 (MQ=255)
  tCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTggg                                                                                                                                                 >  2:4463627/1‑150 (MQ=255)
    ttGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGcc                                                                                                                                               >  1:398593/1‑150 (MQ=255)
      gAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCgg                                                                                                                                             <  1:3641419/150‑2 (MQ=255)
      gAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCgg                                                                                                                                             <  1:382781/150‑2 (MQ=255)
       aaTTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCgtg                                                                                                                                            <  2:2826737/150‑1 (MQ=255)
       aaTTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCACGGCCCAGCGGCTGGGCCggg                                                                                                                                            <  1:862068/150‑3 (MQ=255)
         ttCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGtgtg                                                                                                                                          >  2:2534206/1‑150 (MQ=255)
           cACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGAc                                                                                                                                        <  1:4723683/150‑1 (MQ=255)
                 tGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACg                                                                                                                                  >  1:3317529/1‑150 (MQ=255)
                  gAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACgg                                                                                                                                 >  2:1695279/1‑150 (MQ=255)
                   aaGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGa                                                                                                                                >  2:2247093/1‑150 (MQ=255)
                        tATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGc                                                                                                                           <  1:3697872/150‑1 (MQ=255)
                         aTCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTTACTGGAAAGCCCGGCCCAGcc                                                                                                                                                       <  2:298605/121‑1 (MQ=255)
                         aTCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTGACTGGAAAGCCCGGCCCAGcc                                                                                                                                                       >  1:298605/1‑121 (MQ=255)
                           cTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGc                                                                                                                        <  2:2454019/150‑1 (MQ=255)
                             gTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCaa                                                                                                                      >  2:3113857/1‑150 (MQ=255)
                                caccacCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACgg                                                                                                                   >  1:5358121/1‑150 (MQ=255)
                                caccacCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACgg                                                                                                                   >  2:2116480/1‑150 (MQ=255)
                                caccacCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACgg                                                                                                                   <  2:3640315/150‑1 (MQ=255)
                                caccacCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACgg                                                                                                                   >  1:2631046/1‑150 (MQ=255)
                                 caggccccGACTGCCGGCCGCTATGAAATCCCGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCCGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGgc                                                                                                                  >  2:526660/5‑150 (MQ=255)
                                 accaccCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCgg                                                                                                                                             >  1:4158323/1‑122 (MQ=255)
                                 accaccCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGgc                                                                                                                  <  1:5177279/150‑1 (MQ=255)
                                  ccaccCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGgcg                                                                                                                 <  2:3317529/150‑1 (MQ=255)
                                     cccAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgcag                                                                                                              >  1:959140/1‑150 (MQ=255)
                                       cAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgcagca                                                                                                            >  2:3390426/1‑150 (MQ=255)
                                            tGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTg                                                                                                       <  1:3583637/150‑1 (MQ=255)
                                                ggccgCGATGACATCCCGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACg                                                                                                   <  2:5348403/150‑1 (MQ=255)
                                                      tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                             >  1:3271105/1‑150 (MQ=255)
                                                      tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                             >  1:4824861/1‑150 (MQ=255)
                                                           aaTCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGcc                                                                                        <  2:293257/150‑1 (MQ=255)
                                                            aTCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCg                                                                                       >  1:3207861/1‑150 (MQ=255)
                                                            aTCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCg                                                                                       <  2:1001021/150‑1 (MQ=255)
                                                              ccAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGc                                                                                     >  1:4617542/1‑150 (MQ=255)
                                                                aGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCa                                                                                   >  1:3358263/1‑150 (MQ=255)
                                                                  gTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCAcc                                                                                 >  1:2920481/1‑150 (MQ=255)
                                                                  gTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCAcc                                                                                 >  2:4906785/1‑150 (MQ=255)
                                                                      aCCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGc                                                                             <  2:1598235/150‑1 (MQ=255)
                                                                       cccGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCg                                                                            >  1:2035134/1‑150 (MQ=255)
                                                                         cGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGAc                                                                          >  1:4940084/1‑150 (MQ=255)
                                                                          gCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACt                                                                         >  2:4297643/1‑150 (MQ=255)
                                                                            tGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACttt                                                                       <  1:3887194/150‑1 (MQ=255)
                                                                            tGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACttt                                                                       >  1:102906/1‑150 (MQ=255)
                                                                                   aCGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGa                                                                >  1:276408/1‑150 (MQ=255)
                                                                                      ccTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgt                                                             >  2:628376/1‑150 (MQ=255)
                                                                                       cTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtg                                                            <  2:2733698/150‑1 (MQ=255)
                                                                                           gcATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtgg                                                        <  1:784310/150‑1 (MQ=255)
                                                                                           gcATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtgg                                                        <  1:2389443/150‑1 (MQ=255)
                                                                                            cATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggt                                                       <  1:4045836/150‑1 (MQ=255)
                                                                                              tGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCACGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggtgg                                                     <  1:5257715/150‑1 (MQ=255)
                                                                                                 cGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGa                                                  >  2:4746440/1‑150 (MQ=255)
                                                                                                 cGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGa                                                  <  1:2525960/150‑1 (MQ=255)
                                                                                                 cGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGa                                                  >  2:2007248/1‑150 (MQ=255)
                                                                                                  gTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAc                                                 >  1:2672399/1‑150 (MQ=255)
                                                                                                        cagccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAcccgccc                                           <  1:3625311/150‑1 (MQ=255)
                                                                                                         agccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAcccgcccg                                          <  2:1850060/150‑1 (MQ=255)
                                                                                                           ccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtg                                                         <  2:3262386/133‑1 (MQ=255)
                                                                                                           ccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCt                                        <  2:3948944/150‑1 (MQ=255)
                                                                                                             agcagcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCa                                      >  1:1788847/1‑150 (MQ=255)
                                                                                                               cagcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAg                                    <  2:1203429/150‑1 (MQ=255)
                                                                                                                agcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAgc                                   <  2:1414606/150‑1 (MQ=255)
                                                                                                                 gcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAgcg                                  >  1:3938688/1‑150 (MQ=255)
                                                                                                                  cTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAgcgc                                 <  1:115133/150‑1 (MQ=255)
                                                                                                                   tGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGcgcg                                >  1:3455874/1‑150 (MQ=255)
                                                                                                                    ggACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGa                               <  2:5394400/150‑1 (MQ=255)
                                                                                                                    ggACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGa                               >  1:4454622/1‑150 (MQ=255)
                                                                                                                          cTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAacac                         >  2:3888457/1‑150 (MQ=255)
                                                                                                                          cTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAacac                         >  1:1849913/1‑150 (MQ=255)
                                                                                                                           tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggtg                                                      <  2:3159286/120‑1 (MQ=255)
                                                                                                                           tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                        >  1:2127389/1‑150 (MQ=255)
                                                                                                                           tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                        >  1:1424833/1‑150 (MQ=255)
                                                                                                                           tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                        >  2:1899895/1‑150 (MQ=255)
                                                                                                                           tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                        >  2:15576/1‑150 (MQ=255)
                                                                                                                           tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                        >  2:5109901/1‑150 (MQ=255)
                                                                                                                             aCTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGt                      >  2:264723/1‑150 (MQ=255)
                                                                                                                             aCTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGt                      >  1:830227/1‑150 (MQ=255)
                                                                                                                               tGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCg                    >  2:4038686/1‑150 (MQ=255)
                                                                                                                               tGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCg                    >  1:3414782/1‑150 (MQ=255)
                                                                                                                                ggAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCgg                   <  1:1574294/150‑1 (MQ=255)
                                                                                                                                 gAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGc                  >  1:4130258/1‑150 (MQ=255)
                                                                                                                                   tcacacGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGc                  <  1:1280172/143‑1 (MQ=255)
                                                                                                                                    aGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGaa               >  1:3156882/1‑150 (MQ=255)
                                                                                                                                     gCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg              <  2:3450022/150‑1 (MQ=255)
                                                                                                                                     gCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg              >  2:5002847/1‑150 (MQ=255)
                                                                                                                                     acacGGCCCAGCGGCTGGGCCGGGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg              <  1:1068690/147‑1 (MQ=255)
                                                                                                                                     acacGGCCCAGCGGCTGGGCCGGGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg              >  2:1068690/4‑150 (MQ=255)
                                                                                                                                      cccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg             >  1:5064670/1‑150 (MQ=255)
                                                                                                                                      cccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg             >  1:1601805/1‑150 (MQ=255)
                                                                                                                                       ccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt            >  1:4317621/1‑150 (MQ=255)
                                                                                                                                       acGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt            >  2:3460960/2‑150 (MQ=255)
                                                                                                                                       acGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGATGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGACCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt            >  1:244002/2‑150 (MQ=255)
                                                                                                                                         ggCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGa          <  1:4793339/150‑1 (MQ=255)
                                                                                                                                         ggCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGa          >  1:2759811/1‑150 (MQ=255)
                                                                                                                                          gCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGAc         <  1:2247093/150‑1 (MQ=255)
                                                                                                                                            ccAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGGATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGt       <  2:4945648/150‑1 (MQ=255)
                                                                                                                                             cAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTg      >  1:366328/1‑150 (MQ=255)
                                                                                                                                                ccGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGcc   >  1:2447553/1‑150 (MQ=255)
                                                                                                                                                 cGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGccc  >  1:2968279/1‑150 (MQ=255)
                                                                                                                                                 cGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGccc  >  1:197306/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGCCC  >  NC_002947/4987873‑4988167

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