Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,501,772 T→C 5.2% intergenic (+48/+24) mdeA → / ← PP_1309 methionine gamma‑lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,501,7720TC 5.2% 810.9 / 20.0 249intergenic (+48/+24)mdeA/PP_1309methionine gamma‑lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (134/101);  new base C (9/4);  total (144/105)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 253 total aligned reads displayed.

CATCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTGC  >  NC_002947/1501623‑1501921
                                                                                                                                                     |                                                                                                                                                     
cATCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGACATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCtt                                                                                                                                                       <  1:3436486/150‑1 (MQ=255)
  tCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGc                                                                                                                                                     >  2:5446788/1‑150 (MQ=255)
    aCGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGa                                                                                                                                                   >  1:3541625/1‑150 (MQ=255)
    aCGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGa                                                                                                                                                   <  1:2897220/150‑1 (MQ=255)
     cGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAg                                                                                                                                                  <  1:3861773/150‑1 (MQ=255)
       ggATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGc                                                                                                                                                <  1:662617/150‑1 (MQ=255)
        gATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGcc                                                                                                                                               <  1:5446788/150‑1 (MQ=255)
           atCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGATCAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCttt                                                                                                                                            <  1:1167950/150‑1 (MQ=255)
             cGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTATGTTGCCCGCGACGAGGCCATTTAATCCAGCCCACAAAAATGTCCTCTTGCGAGGCCTTTGt                                                                                                                                          <  2:1280897/150‑1 (MQ=255)
               gAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTca                                                                                                                                        >  2:4148605/1‑150 (MQ=255)
               gAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTca                                                                                                                                        <  2:3854845/150‑1 (MQ=255)
               gAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTca                                                                                                                                        <  2:3712224/150‑1 (MQ=255)
                aGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAATAAAGGCCTCTTGCGAGGCCTTTGTcac                                                                                                                                       >  1:5196902/1‑150 (MQ=255)
                    gTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTc                                                                                                                                   >  1:4329866/1‑150 (MQ=255)
                    gTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTc                                                                                                                                   <  2:2527085/150‑1 (MQ=255)
                     ttGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCa                                                                                                                                  >  1:1500866/1‑150 (MQ=255)
                     ttGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCa                                                                                                                                  >  2:1173388/1‑150 (MQ=255)
                       ggTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGa                                                                                                                                <  1:2494865/150‑1 (MQ=255)
                           cGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCaa                                                                                                                            <  1:2084088/150‑1 (MQ=255)
                            gTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAAc                                                                                                                           >  2:2417323/1‑150 (MQ=255)
                             tCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACg                                                                                                                          >  2:1113761/1‑150 (MQ=255)
                                gTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGc                                                                                                                                                     >  2:4274337/1‑120 (MQ=255)
                                  cGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGa                                                                                                                     >  1:1612957/1‑150 (MQ=255)
                                      ggggCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGccc                                                                                                                 >  1:4819676/1‑150 (MQ=255)
                                       gggCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCg                                                                                                                <  1:5051291/150‑1 (MQ=255)
                                       gggCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCg                                                                                                                <  1:3286464/150‑1 (MQ=255)
                                                gaTGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcg                                                                                                       >  1:2984503/1‑150 (MQ=255)
                                                 aTGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcgt                                                                                                      <  1:4148605/150‑1 (MQ=255)
                                                 aTGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCGGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCACCCCACACAAAAGGCCTCTTGCGAGGCCATTGTCACATACAGATCAACGCTTGAGCCCGTAGCTTtcgt                                                                                                      <  2:1958975/150‑1 (MQ=255)
                                                    tgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGa                                                                                                   <  2:3260630/150‑1 (MQ=255)
                                                    tgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGa                                                                                                   <  1:1043539/150‑1 (MQ=255)
                                                     ggaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAg                                                                                                  >  2:3558436/1‑150 (MQ=255)
                                                      gaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGc                                                                                                 >  1:3657954/1‑150 (MQ=255)
                                                      gaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGc                                                                                                 <  1:828305/150‑1 (MQ=255)
                                                             tgctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCg                                                                                          <  1:3375331/150‑1 (MQ=255)
                                                                 ggCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGa                                                                                      <  1:3481564/150‑1 (MQ=255)
                                                                    aGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGtt                                                                                   >  2:2645100/1‑150 (MQ=255)
                                                                    aGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGtt                                                                                   >  1:2409503/1‑150 (MQ=255)
                                                                       tGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAgg                                                                                >  1:5127090/1‑150 (MQ=255)
                                                                        gTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGgc                                                                               >  2:622537/1‑150 (MQ=255)
                                                                         tCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGgcg                                                                              >  1:2652417/1‑150 (MQ=255)
                                                                           gAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCt                                                                            >  2:1804317/1‑150 (MQ=255)
                                                                           gAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCt                                                                            >  1:1201758/1‑150 (MQ=255)
                                                                            aaCAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCtt                                                                           >  2:2201042/1‑150 (MQ=255)
                                                                               aGGCATTGCAGGCCTGCGGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCg                                                                        >  2:334703/1‑150 (MQ=255)
                                                                                 gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGt                                                                      >  2:4750478/1‑150 (MQ=255)
                                                                                 gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGt                                                                      >  2:4287597/1‑150 (MQ=255)
                                                                                   aTTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAAAGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGc                                                                    >  1:5133642/1‑150 (MQ=255)
                                                                                    ttGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCa                                                                   <  2:4794362/150‑1 (MQ=255)
                                                                                      gCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATa                                                                 >  2:3613932/1‑150 (MQ=255)
                                                                                       cAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAg                                                                >  2:521346/1‑150 (MQ=255)
                                                                                        aGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGt                                                               >  1:1908856/1‑150 (MQ=255)
                                                                                        aGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGt                                                               >  2:355249/1‑150 (MQ=255)
                                                                                             tgctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgc                                                          >  2:5011627/1‑150 (MQ=255)
                                                                                                tgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGc                                                       >  1:848929/1‑150 (MQ=255)
                                                                                                tgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGc                                                       >  2:2012848/1‑150 (MQ=255)
                                                                                                 gTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAAc                                                                                                                           >  1:4696648/1‑81 (MQ=255)
                                                                                                 gTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAAc                                                                                                                           <  2:4696648/81‑1 (MQ=255)
                                                                                                 gTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCtt                                                                                                                       <  2:372332/85‑1 (MQ=255)
                                                                                                 gTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCtt                                                                                                                       >  1:372332/1‑85 (MQ=255)
                                                                                                 gTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCa                                                      >  1:584810/1‑150 (MQ=255)
                                                                                                 gTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCa                                                      >  1:1433367/1‑150 (MQ=255)
                                                                                                 gTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCa                                                      >  1:4219615/1‑150 (MQ=255)
                                                                                                 gTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCa                                                      >  2:1794034/1‑150 (MQ=255)
                                                                                                  ttAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCAc                                                     >  2:4500594/1‑150 (MQ=255)
                                                                                                         gCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcggcg                                              <  1:3041878/150‑1 (MQ=255)
                                                                                                           ccGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgt                                            >  2:2879336/1‑150 (MQ=255)
                                                                                                           ccGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgt                                            <  2:3471722/150‑1 (MQ=255)
                                                                                                           ccGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgt                                            >  2:5353733/1‑150 (MQ=255)
                                                                                                            cgcgAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgtg                                           >  1:2294963/1‑150 (MQ=255)
                                                                                                            cgcgAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgtg                                           >  2:1386512/1‑150 (MQ=255)
                                                                                                               gAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCt                                        <  1:4016141/150‑1 (MQ=255)
                                                                                                                 agagGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCttt                                      <  1:4373245/150‑1 (MQ=255)
                                                                                                                 agagGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCttt                                      <  2:351353/150‑1 (MQ=255)
                                                                                                                   agGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAg                                    <  2:3150415/150‑1 (MQ=255)
                                                                                                                   agGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAg                                    <  2:3689678/150‑1 (MQ=255)
                                                                                                                       cATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCACGTAGCTTTCGTCGAGCATGACAGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgg                                >  1:4465341/1‑150 (MQ=255)
                                                                                                                         ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgggg                              >  1:5334627/1‑150 (MQ=255)
                                                                                                                         ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgggg                              >  1:2204788/1‑150 (MQ=255)
                                                                                                                          tGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGt                             <  1:1302963/150‑1 (MQ=255)
                                                                                                                            aaTCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCa                           >  1:258493/1‑150 (MQ=255)
                                                                                                                              tCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcggcg                                              <  1:5522880/129‑1 (MQ=255)
                                                                                                                              tCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAgg                         >  2:1070275/1‑150 (MQ=255)
                                                                                                                              tCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGCTGGGGTCAgg                         >  1:3513122/1‑150 (MQ=255)
                                                                                                                               ccAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGc                        <  2:1418243/150‑1 (MQ=255)
                                                                                                                                cAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCg                       >  2:451917/1‑150 (MQ=255)
                                                                                                                                 aGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCgg                      <  1:1113810/150‑1 (MQ=255)
                                                                                                                                    ccACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCa                   >  1:1568000/1‑150 (MQ=255)
                                                                                                                                     cacaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAc                  >  2:4847680/1‑150 (MQ=255)
                                                                                                                                      acaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAgg                                                                                <  2:2744081/87‑1 (MQ=255)
                                                                                                                                      acaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCagt                 >  2:569490/1‑148 (MQ=255)
                                                                                                                                        gacaaaGGCCTCGCAAGAGGCCTTTTTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcgcg               >  2:2112203/4‑150 (MQ=255)
                                                                                                                                           aaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc            <  1:179798/150‑1 (MQ=255)
                                                                                                                                           aaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc            >  2:717570/1‑150 (MQ=255)
                                                                                                                                           aaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGTGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc            >  1:5213720/1‑150 (MQ=255)
                                                                                                                                           aaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc            >  1:617194/1‑150 (MQ=255)
                                                                                                                                                  tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCgg     >  2:383500/1‑150 (MQ=255)
                                                                                                                                                  tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCgg     >  2:2956234/1‑150 (MQ=255)
                                                                                                                                                  tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCgg     <  1:4287477/150‑1 (MQ=255)
                                                                                                                                                   ctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGt    <  1:5511414/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTGC  >  NC_002947/1501623‑1501921

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: