Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,986,096 C→G 7.2% 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,986,0960CG 7.2% 345.9 / 13.8 113P97R (CCG→CGG) flgAflagella basal body P‑ring formation protein
Reads supporting (aligned to +/- strand):  ref base C (54/50);  new base G (6/2);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 113 total aligned reads displayed.

TTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGG  >  NC_002947/4985955‑4986242
                                                                                                                                             |                                                                                                                                                   
ttGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCGGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGcc                                                                                                                                             >  2:1735485/1‑150 (MQ=255)
 tGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCg                                                                                                                                            >  1:996972/1‑150 (MQ=255)
  gAATTCACCGTTGAAGATTATCTGTCCATCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCgg                                                                                                                                           <  2:2034825/150‑2 (MQ=255)
  gAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCgt                                                                                                                                           <  2:1358656/150‑1 (MQ=255)
     ttCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGtgtg                                                                                                                                        >  2:232224/1‑150 (MQ=255)
      tCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGa                                                                                                                                       <  1:2123572/150‑1 (MQ=255)
      tCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGa                                                                                                                                       <  2:13438/150‑1 (MQ=255)
      tCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGa                                                                                                                                       <  2:258578/150‑1 (MQ=255)
        aCCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTg                                                                                                                                                 >  2:1656174/1‑138 (MQ=255)
        aCCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTg                                                                                                                                                 <  1:1656174/138‑1 (MQ=255)
        aCCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGGCg                                                                                                                                     >  2:48317/1‑150 (MQ=255)
        aCCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACg                                                                                                                                     >  2:589822/1‑150 (MQ=255)
              gAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACgg                                                                                                                               <  1:1628864/150‑1 (MQ=255)
                 gATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATa                                                                                                                           >  2:1032724/1‑150 (MQ=255)
                 gATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATa                                                                                                                           >  1:50598/1‑150 (MQ=255)
                 gATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATa                                                                                                                           >  2:483691/1‑150 (MQ=255)
                  aTTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATAc                                                                                                                          <  1:483691/150‑1 (MQ=255)
                    tATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGc                                                                                                                        >  1:1004486/1‑150 (MQ=255)
                        tGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCa                                                                                                                    >  1:2207164/1‑150 (MQ=255)
                            caccacCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCggg                                                                                                                                          <  2:1847826/125‑3 (MQ=255)
                            caccacCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCggg                                                                                                                                          >  1:1847826/1‑123 (MQ=255)
                              ccaccCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGgcg                                                                                                              <  1:795418/150‑1 (MQ=255)
                                 cccAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGCTACGCGGCGACGGCgcag                                                                                                           >  2:2243453/1‑150 (MQ=255)
                                  ccAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgcagc                                                                                                          >  2:2276765/1‑150 (MQ=255)
                                       cTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGt                                                                                                     <  1:1707642/150‑1 (MQ=255)
                                        tGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGTCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTg                                                                                                    <  1:1032724/150‑1 (MQ=255)
                                          ccggccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTTGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTg                                                                                                                      <  1:1641023/130‑1 (MQ=255)
                                          ccggccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTg                                                                                                                      >  2:1641023/1‑130 (MQ=255)
                                          ccggccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGgcg                                                                                                              <  1:2006311/138‑1 (MQ=255)
                                          ccggccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGgcg                                                                                                              >  2:2006311/1‑138 (MQ=255)
                                           cggccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTg                                                                                                                      <  1:1803766/129‑1 (MQ=255)
                                           cggccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTg                                                                                                                      >  2:1803766/1‑129 (MQ=255)
                                           cggccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGAc                                                                                                 >  1:1355781/1‑150 (MQ=255)
                                               cgCTATGAGATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCgt                                                                                                                                           <  2:1435618/105‑1 (MQ=255)
                                               cgCTATGAGATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCgt                                                                                                                                           >  1:1435618/1‑105 (MQ=255)
                                                  tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGCCGCATCCCTGGAAAGCACGCAAAAGCTGCTGGGCCGTGTGATTGTACGGTTTCTCTTGAACGGCGCATTACCGTGGACGGTATTc                                                                                          <  2:2081109/150‑1 (MQ=255)
                                                    tGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGt                                                                                        >  2:2160849/1‑150 (MQ=255)
                                                      aaaTCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGc                                                                                      <  1:1474146/150‑1 (MQ=255)
                                                         tCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCTTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGgc                                                                                                               <  1:1192534/122‑1 (MQ=255)
                                                         tCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCgg                                                                                   >  2:1783412/1‑150 (MQ=255)
                                                         tCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCgg                                                                                   >  1:575625/1‑150 (MQ=255)
                                                         tCCAGGTCAACACGCTGGACCCACGCCTGCGCATGCCTTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGgc                                                                                                               >  2:1192534/1‑122 (MQ=255)
                                                           cAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGcc                                                                                 >  2:388285/1‑150 (MQ=255)
                                                             ggTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCAc                                                                               >  1:386186/1‑150 (MQ=255)
                                                              gTCAACCCGCTGGCCCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGACGTGTGACGGGCGGAGTACGCTCCAGCGTCGTGGTATCATGGATGGTATGCTTTCCTGCCAcc                                                                              >  2:1129933/1‑150 (MQ=255)
                                                              gTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCAcc                                                                              >  1:1590736/1‑150 (MQ=255)
                                                                  aCCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGc                                                                          <  1:257269/150‑1 (MQ=255)
                                                                   cccGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTggg                                                                                                                                               <  1:2133412/81‑1 (MQ=255)
                                                                   cccGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCg                                                                         >  1:2133120/1‑150 (MQ=255)
                                                                          gACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTc                                                                  <  1:1322247/150‑1 (MQ=255)
                                                                           aCCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTcc                                                                 <  2:386186/150‑1 (MQ=255)
                                                                            cccACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCcg                                                                <  1:2276765/150‑1 (MQ=255)
                                                                                 gCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACg                                                           >  1:1980475/1‑150 (MQ=255)
                                                                                   cTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtg                                                         <  2:2207164/150‑1 (MQ=255)
                                                                                   cTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtg                                                         >  2:914062/1‑150 (MQ=255)
                                                                                    tGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtgg                                                        >  2:1995659/1‑150 (MQ=255)
                                                                                         aTGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggtg                                                   <  1:534925/150‑1 (MQ=255)
                                                                                          tGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGGCGGCCACCGTGCGACTTTTCGGCGTCgtggtggtgg                                                  >  2:384349/1‑150 (MQ=255)
                                                                                          tGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggtgg                                                  <  1:1451918/150‑1 (MQ=255)
                                                                                          tGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggtgg                                                  <  1:388285/150‑1 (MQ=255)
                                                                                          tGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCCAGTACTTTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggtgg                                                  <  1:2243453/150‑1 (MQ=255)
                                                                                           gCCGTTGTGCAGCCAGCAGCGGGACGCCTCACTGGAAAGCCCGGACCAGCCACGCGGCCGTGTGACGGTACGGATACGCTGCAACGGCGGAGCACCGTGGACGGTCTGCGTGCCGGCCATCGTGCGACTTTTCCGCGACgtggtggtggt                                                 <  2:677084/150‑1 (MQ=255)
                                                                                             cGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGa                                               <  1:2325155/150‑1 (MQ=255)
                                                                                              gTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAc                                              <  1:589822/150‑1 (MQ=255)
                                                                                                 gtgCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAcccg                                           <  1:1155033/150‑1 (MQ=255)
                                                                                                      gccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCt                                                                                                                       >  2:2247119/1‑69 (MQ=255)
                                                                                                      gccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCt                                                                                                                       <  1:2247119/69‑1 (MQ=255)
                                                                                                         agcagcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCTACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCa                                   >  2:221709/1‑150 (MQ=255)
                                                                                                              cTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggtgg                                                  <  2:75183/130‑1 (MQ=255)
                                                                                                              cTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtgggggtgg                                                  >  1:75183/1‑130 (MQ=255)
                                                                                                              cTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAACCGCCCGCTCAAgcgc                              <  2:1355781/150‑1 (MQ=255)
                                                                                                              cTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGGATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAgcgc                              <  1:554636/150‑1 (MQ=255)
                                                                                                                 gACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGac                           >  1:1670779/1‑150 (MQ=255)
                                                                                                                      cTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtg                                                         >  1:873723/1‑115 (MQ=255)
                                                                                                                      cTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtg                                                         <  2:873723/115‑1 (MQ=255)
                                                                                                                       tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                     >  2:325057/1‑150 (MQ=255)
                                                                                                                       tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                     >  2:844855/1‑150 (MQ=255)
                                                                                                                             gAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGc                                                                          >  1:1719191/1‑91 (MQ=255)
                                                                                                                             gAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGc                                                                          <  2:1719191/91‑1 (MQ=255)
                                                                                                                             gAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGc               <  1:1033030/150‑1 (MQ=255)
                                                                                                                             gAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACAGTCGGc               >  2:1670452/1‑150 (MQ=255)
                                                                                                                              aaaGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGacaac                        >  2:653262/1‑140 (MQ=255)
                                                                                                                              aaaGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGacaac                        <  1:653262/140‑1 (MQ=255)
                                                                                                                               aaGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGa             <  2:1980475/150‑1 (MQ=255)
                                                                                                                                 gCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg           >  1:2373672/1‑150 (MQ=255)
                                                                                                                                 acacGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg           <  1:1122555/147‑1 (MQ=255)
                                                                                                                                  ccgggCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg          >  1:2402049/4‑150 (MQ=255)
                                                                                                                                  cacGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCACCGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCAATCCCGCGCGGCAACACGGTCGGCGAAgg          >  2:1122555/3‑150 (MQ=255)
                                                                                                                                  cacGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg          >  2:2126734/3‑150 (MQ=255)
                                                                                                                                   ccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt         >  1:2128217/1‑150 (MQ=255)
                                                                                                                                   ccGGCACAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACGGTGCGACTTTTCCGCGTCGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt         <  1:821777/150‑1 (MQ=255)
                                                                                                                                   acGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt         >  2:2281566/2‑150 (MQ=255)
                                                                                                                                   acGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt         >  1:1115987/2‑150 (MQ=255)
                                                                                                                                   acGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt         >  1:1751711/2‑150 (MQ=255)
                                                                                                                                     ggCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGCTGa       <  2:34411/150‑1 (MQ=255)
                                                                                                                                        ccAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAccc                                            >  1:1092827/1‑110 (MQ=255)
                                                                                                                                        ccAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAccc                                            <  2:1092827/110‑1 (MQ=255)
                                                                                                                                         cAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGacaaca                       >  1:2088749/1‑130 (MQ=255)
                                                                                                                                         cAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGacaaca                       <  2:2088749/130‑1 (MQ=255)
                                                                                                                                          aGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTgg  <  1:1735485/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGG  >  NC_002947/4985955‑4986242

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: