Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,501,774 C→A 5.3% intergenic (+50/+22) 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,7740CA 5.3% 835.9 / 25.4 246intergenic (+50/+22)mdeA/PP_1309methionine gamma‑lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (132/101);  new base A (9/4);  total (141/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.15e-01

Only 100 of 251 total aligned reads displayed.

TCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTGC  >  NC_002947/1501625‑1501921
                                                                                                                                                     |                                                                                                                                                   
tCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGc                                                                                                                                                     >  2:5446788/1‑150 (MQ=255)
   cGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAg                                                                                                                                                  <  1:3861899/150‑1 (MQ=255)
     ggATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGc                                                                                                                                                <  1:662109/150‑1 (MQ=255)
      gATATCGGAGGGGTTGGTGCTTCTGTCGGTGGGGGTGGAGGATGTGGAGGACCTGCTGGGAGATGTCGAAGAGGGATTGCAGGCGTGGTGTTAGGTTGCCCGCGAAGAGGCGATTGAATCCAGCCGACTATAATGGCCTCTTGCGAgggg                                                                                                                                               >  2:3187159/1‑148 (MQ=255)
         atCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGATCAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCttt                                                                                                                                            <  1:1167950/150‑1 (MQ=255)
         atCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCttt                                                                                                                                            <  1:2515810/150‑1 (MQ=255)
          tCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCGCAAGAGGCCtttt                                                                                                                                           <  2:3349906/150‑2 (MQ=255)
          tCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCAATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTg                                                                                                                                           <  2:3917956/150‑1 (MQ=255)
            ggAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTc                                                                                                                                         <  2:2188285/150‑1 (MQ=255)
            ggAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCGCAAGAGGCCtttttt                                                                                                                                         <  2:2032873/150‑4 (MQ=255)
             gAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTca                                                                                                                                        >  2:4148605/1‑150 (MQ=255)
             gAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTca                                                                                                                                        <  2:3712224/150‑1 (MQ=255)
              aGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCGCAAGAGGCCttttttgt                                                                                                                                       <  2:3598340/150‑6 (MQ=255)
                  gTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACTGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTc                                                                                                                                   >  2:179798/1‑150 (MQ=255)
                  gTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTc                                                                                                                                   >  1:2594197/1‑150 (MQ=255)
                  gTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTc                                                                                                                                   <  2:2527185/150‑1 (MQ=255)
                   ttGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCa                                                                                                                                  >  1:1500866/1‑150 (MQ=255)
                      gggCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGAt                                                                                                                               >  1:3791608/3‑150 (MQ=255)
                      gTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGAt                                                                                                                               >  1:4573885/1‑150 (MQ=255)
                         cGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCaa                                                                                                                            <  1:2084088/150‑1 (MQ=255)
                           tCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACg                                                                                                                          >  2:1113810/1‑150 (MQ=255)
                             tGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGACCACAAAAAAGGCATCTTGCGAGGCCTTTGTCACATTCAGATAAACGCt                                                                                                                        >  1:5226043/1‑150 (MQ=255)
                              gTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGc                                                                                                                                                     >  2:4274337/1‑120 (MQ=255)
                              gTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGc                                                                                                                                                     <  1:4274337/120‑1 (MQ=255)
                                cGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGa                                                                                                                     <  1:3513346/150‑1 (MQ=255)
                                cGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGa                                                                                                                     >  1:1612957/1‑150 (MQ=255)
                                     gggCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCg                                                                                                                <  1:3286464/150‑1 (MQ=255)
                                            aggaTGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCttt                                                                                                         >  1:4430042/1‑150 (MQ=255)
                                              gaTGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcg                                                                                                       >  1:2984503/1‑150 (MQ=255)
                                                tgtgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGTGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcgtc                                                                                                     <  2:1816036/150‑1 (MQ=255)
                                                  tgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGa                                                                                                   <  1:1043539/150‑1 (MQ=255)
                                                  tgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGa                                                                                                   <  2:3260630/150‑1 (MQ=255)
                                                    gaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGc                                                                                                 >  1:3689678/1‑150 (MQ=255)
                                                    gaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGc                                                                                                 >  2:1289793/1‑150 (MQ=255)
                                                     aggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCa                                                                                                >  1:4357377/1‑150 (MQ=255)
                                                             ctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGc                                                                                        <  1:3394118/150‑1 (MQ=255)
                                                             ctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGc                                                                                        <  1:5203557/150‑1 (MQ=255)
                                                              tgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGa                                                                                                                     >  2:1272259/1‑120 (MQ=255)
                                                              tgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCg                                                                                       >  1:5109260/1‑150 (MQ=255)
                                                                   gATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTa                                                                                  <  2:1089926/150‑1 (MQ=255)
                                                                      gTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGgc                                                                               >  2:622537/1‑150 (MQ=255)
                                                                       tCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGgcg                                                                              >  1:2652417/1‑150 (MQ=255)
                                                                         gAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCt                                                                            >  1:1201758/1‑150 (MQ=255)
                                                                               gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGt                                                                      >  2:4750478/1‑150 (MQ=255)
                                                                               gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGt                                                                      >  2:4287597/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)
                                                                                      aGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGt                                                               >  1:4163338/1‑150 (MQ=255)
                                                                                          ctgctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcg                                                           >  1:5539823/1‑150 (MQ=255)
                                                                                           tgctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGa                                                                                                                     <  1:4779208/91‑1 (MQ=255)
                                                                                           tgctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgc                                                          >  2:5011627/1‑150 (MQ=255)
                                                                                            gctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgcg                                                         >  2:2905497/1‑150 (MQ=255)
                                                                                              tgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGc                                                       >  1:848856/1‑150 (MQ=255)
                                                                                               gTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAAc                                                                                                                           <  2:4696648/81‑1 (MQ=255)
                                                                                               gTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCa                                                      >  1:1433367/1‑150 (MQ=255)
                                                                                                ttAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCAc                                                     >  1:5132676/1‑150 (MQ=255)
                                                                                                ttAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCAc                                                     >  1:5446815/1‑150 (MQ=255)
                                                                                                ttAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCAc                                                     >  2:4500594/1‑150 (MQ=255)
                                                                                                  aGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACtt                                                   >  1:5518499/1‑150 (MQ=255)
                                                                                                      tGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcggc                                               <  2:3757979/150‑1 (MQ=255)
                                                                                                        cccGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcggcgg                                             >  1:3626913/1‑150 (MQ=255)
                                                                                                         ccGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgt                                            >  2:5353733/1‑150 (MQ=255)
                                                                                                         ccGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgt                                            <  2:3471722/150‑1 (MQ=255)
                                                                                                            cgAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTc                                         <  1:5364330/150‑1 (MQ=255)
                                                                                                               agagGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCttt                                      <  1:2905497/150‑1 (MQ=255)
                                                                                                                   gCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAggg                                  >  2:2488601/1‑149 (MQ=255)
                                                                                                                      aTTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTggg                               >  1:2305533/1‑150 (MQ=255)
                                                                                                                       ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgggg                              >  1:2204788/1‑150 (MQ=255)
                                                                                                                       ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgggg                              >  2:1729942/1‑150 (MQ=255)
                                                                                                                        tGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGt                             <  1:1302963/150‑1 (MQ=255)
                                                                                                                          aaTCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAg                                                                                 <  1:4947507/96‑1 (MQ=255)
                                                                                                                            tCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAgg                         >  1:3513042/1‑150 (MQ=255)
                                                                                                                             ccAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGc                        <  2:1418166/150‑1 (MQ=255)
                                                                                                                              cAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCg                       >  2:451917/1‑150 (MQ=255)
                                                                                                                               aGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCgg                      >  1:2427599/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)
                                                                                                                                    acaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcg                 <  1:3479196/150‑1 (MQ=255)
                                                                                                                                    acaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcg                 >  1:4115552/1‑150 (MQ=255)
                                                                                                                                      gacaaaGGCCTCGCAAGAGGCCTTTTTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcgcg               >  1:5158363/4‑150 (MQ=255)
                                                                                                                                      aaaaaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcgcg               >  1:80010/1‑150 (MQ=255)
                                                                                                                                      aaaaaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCGTTAGGTGGGGTCAGGCGGTCAcgcg               >  2:1993593/1‑150 (MQ=255)
                                                                                                                                        caaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCggg             >  2:3697414/2‑150 (MQ=255)
                                                                                                                                         aaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc            <  1:179798/150‑1 (MQ=255)
                                                                                                                                         aaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc            >  2:717570/1‑150 (MQ=255)
                                                                                                                                         aaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc            >  2:2261717/1‑150 (MQ=255)
                                                                                                                                             gCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTg        >  1:970151/1‑150 (MQ=255)
                                                                                                                                             gCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTg        >  1:611191/1‑150 (MQ=255)
                                                                                                                                               ctctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCg      <  2:3541625/150‑1 (MQ=255)
                                                                                                                                                tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCgg     >  2:383500/1‑150 (MQ=255)
                                                                                                                                                tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCgg     >  2:2956234/1‑150 (MQ=255)
                                                                                                                                                tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCgg     <  1:4287597/150‑1 (MQ=255)
                                                                                                                                                 ctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGt    <  2:2652417/150‑1 (MQ=255)
                                                                                                                                                   tGCGAGGACTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTAGAGCATGCCAGGCGAGTTAGGCGATTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTGc  <  1:3187159/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTGC  >  NC_002947/1501625‑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: