Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,990,844:1 +C 5.8% coding (89/753 nt) 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,990,8441.C 5.8% 1440.5 / 75.4 328coding (89/753 nt)flgAflagella basal body P‑ring formation protein
Reads supporting (aligned to +/- strand):  ref base . (139/170);  new base C (7/12);  total (146/182)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 329 total aligned reads displayed.

TTCAGCCTTTCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCT  >  NC_002947/4990696‑4990988
                                                                                                                                                     |                                                                                                                                                
 tCAGCCTTTCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCa                                                                                                                                                >  1:3601186/1‑150 (MQ=255)
    gCCTTTCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACg                                                                                                                                              <  2:1370696/149‑1 (MQ=255)
     ccTTTCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGt                                                                                                                                            <  2:3758131/150‑1 (MQ=255)
        ttCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGc                                                                                                                                         <  2:5816768/150‑1 (MQ=255)
           gTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGaa                                                                                                                                      <  2:3412984/150‑1 (MQ=255)
            ttGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAAc                                                                                                                                     >  2:6791605/1‑150 (MQ=255)
            ttGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAAc                                                                                                                                     <  1:653955/150‑1 (MQ=255)
            ttGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAAc                                                                                                                                     >  2:4815599/1‑150 (MQ=255)
                aaCTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGTGCGAAcgccg                                                                                                                                 <  2:5179225/150‑3 (MQ=255)
                aaCTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGTGCGAAcgccg                                                                                                                                 <  1:275248/150‑3 (MQ=255)
                aaCTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAAcgccg                                                                                                                                 <  1:4371635/150‑3 (MQ=255)
                aaCTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGt                                                                                                                                 >  1:6174112/1‑150 (MQ=255)
                 aCTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGCGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGCCACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGtt                                                                                                                                <  2:3351194/150‑1 (MQ=255)
                 aCTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGTGCGAAcgccgg                                                                                                                                <  2:1806331/150‑4 (MQ=255)
                 aCTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGtt                                                                                                                                >  1:2272155/1‑150 (MQ=255)
                  cTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGTGCGAAcgccggg                                                                                                                               <  1:6859164/150‑5 (MQ=255)
                     aGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCAcc                                                                                                                            >  2:5294497/1‑150 (MQ=255)
                         aGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGc                                                                                                                        >  1:4855312/1‑150 (MQ=255)
                          gCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGcc                                                                                                                       >  1:692778/1‑150 (MQ=255)
                              gTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGaa                                                                                                                   <  1:347196/150‑1 (MQ=255)
                               tCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAAc                                                                                                                  <  2:3227738/150‑1 (MQ=255)
                               tCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAAc                                                                                                                  <  2:4140743/150‑1 (MQ=255)
                                  gCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGc                                                                                                               <  2:3471762/150‑1 (MQ=255)
                                   cACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCt                                                                                                              <  1:6380836/150‑1 (MQ=255)
                                    aCGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGTGCGAAcgccg                                                                                                                                 >  1:6311905/1‑128 (MQ=255)
                                     cGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTa                                                                                                            <  2:2503170/150‑1 (MQ=255)
                                          tttgctttTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGgt                                                                                                       <  1:674352/150‑1 (MQ=255)
                                             gctttTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTc                                                                                                    >  1:3529716/1‑150 (MQ=255)
                                               ttttttGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCAc                                                                                                  >  2:5651453/1‑150 (MQ=255)
                                                   ttGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCaa                                                                                              <  1:5294806/150‑1 (MQ=255)
                                                    tGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAg                                                                                             >  1:4061506/1‑150 (MQ=255)
                                                    tGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAg                                                                                             <  1:1835122/150‑1 (MQ=255)
                                                     gAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAgg                                                                                            <  1:2671744/150‑1 (MQ=255)
                                                        tgtCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGtt                                                                                         >  1:2006367/1‑150 (MQ=255)
                                                           cATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCt                                                                                      <  1:1171469/150‑1 (MQ=255)
                                                            aTGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCtt                                                                                     <  2:15509/150‑1 (MQ=255)
                                                            aTGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCtt                                                                                     <  1:3303179/150‑1 (MQ=255)
                                                            aTGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCtt                                                                                     <  1:5117229/150‑1 (MQ=255)
                                                            aTGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCtt                                                                                     <  2:6357089/150‑1 (MQ=255)
                                                             tGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTg                                                                                    >  1:1656919/1‑150 (MQ=255)
                                                              gCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGa                                                                                   <  1:5607258/150‑1 (MQ=255)
                                                               cacaCGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGaa                                                                                  <  1:77008/150‑1 (MQ=255)
                                                                acacGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAAt                                                                                 >  2:5750603/1‑150 (MQ=255)
                                                                acacGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAAt                                                                                 >  1:990958/1‑150 (MQ=255)
                                                                    gAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCCACCGTGCGAAcgccg                                                                                                                                 <  1:778520/98‑3 (MQ=255)
                                                                    gAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCAc                                                                             <  1:5526831/150‑1 (MQ=255)
                                                                        aCGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGtt                                                                         <  1:3795972/150‑1 (MQ=255)
                                                                          gACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGa                                                                       >  1:5407252/1‑150 (MQ=255)
                                                                            cATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAg                                                                     >  2:2285147/1‑150 (MQ=255)
                                                                             aTTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGa                                                                    <  1:4686170/150‑1 (MQ=255)
                                                                             aTTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGa                                                                    >  1:1017335/1‑150 (MQ=255)
                                                                              ttttCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGAt                                                                   >  1:2050699/1‑150 (MQ=255)
                                                                                 tCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTAt                                                                <  1:6702915/150‑1 (MQ=255)
                                                                                 tCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTAt                                                                <  2:841755/150‑1 (MQ=255)
                                                                                 tCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTAt                                                                <  1:308469/150‑1 (MQ=255)
                                                                                 tCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTAt                                                                <  2:5077070/150‑1 (MQ=255)
                                                                                  cccGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATc                                                               <  1:2543505/150‑1 (MQ=255)
                                                                                   ccGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCt                                                              <  1:5034063/150‑1 (MQ=255)
                                                                                   ccGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCt                                                              >  1:441853/1‑150 (MQ=255)
                                                                                    cgacgaTTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTg                                                             >  2:2420057/1‑150 (MQ=255)
                                                                                    cgacgaTTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCGTGAATTCACCGTTGAAGATTATCTg                                                             >  1:2857026/1‑150 (MQ=255)
                                                                                    cgacgaTTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTCCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTg                                                             >  2:3244282/1‑150 (MQ=255)
                                                                                     gacgaTTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGt                                                            <  2:4315024/150‑1 (MQ=255)
                                                                                     gacgaTTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGt                                                            <  1:5272580/150‑1 (MQ=255)
                                                                                      acgaTTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTc                                                           <  1:5758776/150‑1 (MQ=255)
                                                                                       cgaTTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTcc                                                          <  2:1656919/150‑1 (MQ=255)
                                                                                            gACACGCCTGATGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCAAGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCcaccac                                                     >  2:3925727/1‑150 (MQ=255)
                                                                                               acGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCAcccc                                                  >  2:607397/1‑149 (MQ=255)
                                                                                                cGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAg                                                 <  2:5881661/150‑1 (MQ=255)
                                                                                                 gCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTGTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGa                                                >  1:6127050/1‑150 (MQ=255)
                                                                                                  ccTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGAc                                               <  2:3008985/150‑1 (MQ=255)
                                                                                                      ctgACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTgcc                                           <  2:957696/150‑1 (MQ=255)
                                                                                                         aCTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGtt                                                                         >  2:5912173/1‑117 (MQ=255)
                                                                                                         aCTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGtt                                                                         <  1:5912173/117‑1 (MQ=255)
                                                                                                         aCTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTgccggc                                        <  1:2951825/150‑1 (MQ=255)
                                                                                                             gCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCt                                    >  1:4210893/1‑150 (MQ=255)
                                                                                                              cccGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTa                                   <  1:4983961/150‑1 (MQ=255)
                                                                                                               ccGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTAt                                  >  2:535359/1‑150 (MQ=255)
                                                                                                               ccGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTAt                                  <  1:5463819/150‑1 (MQ=255)
                                                                                                               ccGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCGAt                                  >  2:4272046/1‑150 (MQ=255)
                                                                                                                   tGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGaaa                              <  1:5239605/150‑1 (MQ=255)
                                                                                                                   tGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGaaa                              <  2:2801403/150‑1 (MQ=255)
                                                                                                                    ggCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAAt                             <  2:1229620/150‑1 (MQ=255)
                                                                                                                       cACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCa                          <  2:2950581/150‑1 (MQ=255)
                                                                                                                       cACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCa                          <  1:5837739/150‑1 (MQ=255)
                                                                                                                       cACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCa                          <  2:6875416/150‑1 (MQ=255)
                                                                                                                       cACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCa                          <  2:2076609/150‑1 (MQ=255)
                                                                                                                       cACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGACGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCa                          <  2:774983/150‑1 (MQ=255)
                                                                                                                        aCGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAg                         <  2:3491374/150‑1 (MQ=255)
                                                                                                                        aCGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAg                         <  1:6494854/150‑1 (MQ=255)
                                                                                                                        aCGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAg                         <  1:1398824/150‑1 (MQ=255)
                                                                                                                        aCGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAg                         <  2:537973/150‑1 (MQ=255)
                                                                                                                               gCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAAcc                  <  1:965974/150‑1 (MQ=255)
                                                                                                                                ccTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAAccc                 <  1:6518450/150‑1 (MQ=255)
                                                                                                                                 ctgctgATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCg                >  1:1508014/1‑150 (MQ=255)
                                                                                                                                  tgctgATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGc               >  1:1257018/1‑150 (MQ=255)
                                                                                                                                        tGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGAcc         <  2:6174112/150‑1 (MQ=255)
                                                                                                                                        tGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGAcc         <  2:6532583/150‑1 (MQ=255)
                                                                                                                                        gtgaacgCGTTCGCCACCGTGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTAGGAAATCCAGGTCAACCCGCTGGAcc         >  1:3986713/7‑150 (MQ=255)
                                                                                                                                               cGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCt  <  1:1460407/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TTCAGCCTTTCGTTGCAACTCAGGTAGCGAGTCGGCACGGGCTTTGCTTTTTTGAGTGTCATGCACACGAAAACGACATTTTCCCGACGATTGACACGCCTGCTGACTGGCCCGCTGGCCACGCTGTGCCTGCTGATGCCCGGCGTTCGCACGGTGGCGAACGCGTTCACCTTGCCTGAACAGCTTATCGGTGTCACCCAAGGGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCT  >  NC_002947/4990696‑4990988

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: