Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,501,774 C→A 5.9% 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.9% 1230.0 / 42.9 367intergenic (+50/+22)mdeA/PP_1309methionine gamma‑lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (196/148);  new base A (11/11);  total (207/160)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 371 total aligned reads displayed.

TCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTGCA  >  NC_002947/1501625‑1501922
                                                                                                                                                     |                                                                                                                                                    
tCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGc                                                                                                                                                      <  2:2678547/150‑1 (MQ=255)
  aCGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGa                                                                                                                                                    <  2:4145660/150‑1 (MQ=255)
  aCGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGa                                                                                                                                                    <  2:4594763/150‑1 (MQ=255)
  aCGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGa                                                                                                                                                    >  1:583224/1‑150 (MQ=255)
   cGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAg                                                                                                                                                   <  1:2990/150‑1 (MQ=255)
    gggATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAgg                                                                                                                                                  >  2:2380849/1‑150 (MQ=255)
     ggATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTCAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGc                                                                                                                                                 <  1:8605598/150‑1 (MQ=255)
         atCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCttt                                                                                                                                             >  2:3641644/1‑150 (MQ=255)
         atCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCttt                                                                                                                                             <  1:4722649/150‑1 (MQ=255)
         atCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCGCAAGAGGCCttt                                                                                                                                             <  1:1610569/150‑1 (MQ=255)
            ggAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTc                                                                                                                                          <  1:2534834/150‑1 (MQ=255)
                     ggTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCAGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGa                                                                                                                                 >  1:4470876/1‑150 (MQ=255)
                       tGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATc                                                                                                                               <  2:3776513/150‑1 (MQ=255)
                             tGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCt                                                                                                                         >  1:543252/1‑150 (MQ=255)
                              gTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCtt                                                                                                                        <  2:1151761/150‑1 (MQ=255)
                                cGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGa                                                                                                                      <  1:6414588/150‑1 (MQ=255)
                                  gTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATGGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGc                                                                                                                    <  2:8010065/150‑1 (MQ=255)
                                     gggCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGGTGCCCGCGAAGAGGCCCTTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCg                                                                                                                 <  2:7941998/150‑1 (MQ=255)
                                                tgtgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcgtc                                                                                                      <  1:1162775/150‑1 (MQ=255)
                                                   ggaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAg                                                                                                   >  2:5450356/1‑150 (MQ=255)
                                                   ggaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAg                                                                                                   >  2:2081198/1‑150 (MQ=255)
                                                   ggaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAg                                                                                                   >  2:744956/1‑150 (MQ=255)
                                                    gaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGc                                                                                                  >  1:7136846/1‑150 (MQ=255)
                                                         ccTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGcc                                                                                             >  1:1803137/1‑150 (MQ=255)
                                                          ctgctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGccc                                                                                            >  2:3325011/1‑150 (MQ=255)
                                                             ctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGc                                                                                         <  1:3235511/150‑1 (MQ=255)
                                                               ggCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGa                                                                                       <  2:5174165/150‑1 (MQ=255)
                                                               ggCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGa                                                                                       <  2:5687357/150‑1 (MQ=255)
                                                                gCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAg                                                                                      >  1:6738433/1‑150 (MQ=255)
                                                                  aGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGtt                                                                                    <  1:6643352/150‑1 (MQ=255)
                                                                   gATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTa                                                                                   <  2:262313/150‑1 (MQ=255)
                                                                    aTGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAg                                                                                  >  1:3869145/1‑150 (MQ=255)
                                                                    aTGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAg                                                                                  <  1:1929607/150‑1 (MQ=255)
                                                                     tGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAgg                                                                                 <  2:1496192/150‑1 (MQ=255)
                                                                      gTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGgc                                                                                >  1:940385/1‑150 (MQ=255)
                                                                         gAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCt                                                                             >  2:8189764/1‑150 (MQ=255)
                                                                         gAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCt                                                                             >  1:1229301/1‑150 (MQ=255)
                                                                          aaCAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCtt                                                                            >  1:6314918/1‑150 (MQ=255)
                                                                          aaCAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCtt                                                                            <  2:8362574/150‑1 (MQ=255)
                                                                           aCAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCttt                                                                           >  2:7901644/1‑150 (MQ=255)
                                                                             aGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCg                                                                         >  2:1243262/1‑150 (MQ=255)
                                                                              ggCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCgg                                                                        >  2:2036554/1‑150 (MQ=255)
                                                                               gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAATAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGt                                                                       >  2:507772/1‑150 (MQ=255)
                                                                               gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGt                                                                       >  2:4281674/1‑150 (MQ=255)
                                                                               gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGt                                                                       >  1:3753687/1‑150 (MQ=255)
                                                                               gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGAAGCTTTCGTCGAGCATGCCAGGCGAGTTAGGCGCTTTCGGt                                                                       >  1:108397/1‑150 (MQ=255)
                                                                               gCATTGCAGGCCTGCTGTGAGGTTGCCCGCGAAGAGGCCATTGACTCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGCTAGGCGCTTTCGGt                                                                       <  2:7475285/150‑1 (MQ=255)
                                                                                cATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTg                                                                      <  1:419655/150‑1 (MQ=255)
                                                                                  ttGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCa                                                                    >  2:4787770/1‑150 (MQ=255)
                                                                                    gCAGGCCTGCTTTTAGGTTGCCCGCGAAGAGGTCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATa                                                                  <  1:7720365/150‑1 (MQ=255)
                                                                                        gCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcg                                                              >  1:4910981/1‑150 (MQ=255)
                                                                                        gCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcg                                                              >  2:8604352/1‑150 (MQ=255)
                                                                                         ccTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgc                                                             <  1:5691669/150‑1 (MQ=255)
                                                                                          ctgctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcg                                                            <  2:7251361/150‑1 (MQ=255)
                                                                                            gctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgcg                                                          >  1:3540414/1‑150 (MQ=255)
                                                                                            gctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgcg                                                          >  1:7978508/1‑150 (MQ=255)
                                                                                             ctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTGTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCgg                                                         <  2:6296498/150‑1 (MQ=255)
                                                                                              tgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAATGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGc                                                        >  1:5159367/1‑150 (MQ=255)
                                                                                              tgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGc                                                        >  1:4662823/1‑150 (MQ=255)
                                                                                                ttAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCAc                                                      >  1:5916275/1‑150 (MQ=255)
                                                                                                ttAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCAc                                                      >  1:8595755/1‑150 (MQ=255)
                                                                                                 tAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAg                                                                                  >  2:3483990/1‑121 (MQ=255)
                                                                                                     ttGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcgg                                                 <  1:224311/150‑1 (MQ=255)
                                                                                                      tGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcggc                                                >  1:8277787/1‑150 (MQ=255)
                                                                                                      tGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcggc                                                <  2:7747138/150‑1 (MQ=255)
                                                                                                        cccGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcggcgg                                              >  1:5269167/1‑150 (MQ=255)
                                                                                                         ccGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgt                                             >  2:1822912/1‑150 (MQ=255)
                                                                                                           gcgAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgtgt                                           <  2:761647/150‑1 (MQ=255)
                                                                                                              aaGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCtt                                        >  2:2492783/1‑150 (MQ=255)
                                                                                                                gagTCCAGTGAATCCAGCCCACAAAACAGGCCTCGTGCGAGGCCTTTGTCACATTCAGATCAACTCTTGAGCCCGTAGCTTTCGTCGAGCAGGCCCGGCGAGTTAGGCGCTTTCGGTGCAGAGTCGCGCGGCACTTCGGCGGTGTCTTTa                                      <  2:6539636/150‑1 (MQ=255)
                                                                                                                    ccATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTg                                  <  1:2380849/150‑1 (MQ=255)
                                                                                                                     cATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgg                                 >  2:6874196/1‑150 (MQ=255)
                                                                                                                     cATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgg                                 >  2:6015238/1‑150 (MQ=255)
                                                                                                                      aTTGAATCCGGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTggg                                <  2:3132625/150‑1 (MQ=255)
                                                                                                                      aTTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTggg                                >  1:8183911/1‑150 (MQ=255)
                                                                                                                       ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgggg                               >  2:5570415/1‑150 (MQ=255)
                                                                                                                            tCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAgg                          >  1:162157/1‑150 (MQ=255)
                                                                                                                              cAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCg                        >  2:2989900/1‑150 (MQ=255)
                                                                                                                              cAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCg                        >  2:5009512/1‑150 (MQ=255)
                                                                                                                               aGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCgg                       >  2:3560245/1‑150 (MQ=255)
                                                                                                                               aGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCgg                       <  2:8469677/150‑1 (MQ=255)
                                                                                                                                 cccACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTc                     >  2:231793/1‑150 (MQ=255)
                                                                                                                                  ccACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCa                    >  1:6263239/1‑150 (MQ=255)
                                                                                                                                  ccACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCa                    >  2:2906806/1‑150 (MQ=255)
                                                                                                                                   cacaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAc                   >  2:6844558/1‑150 (MQ=255)
                                                                                                                                   cacaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAc                   >  1:6830220/1‑150 (MQ=255)
                                                                                                                                    acaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcg                  >  2:1725989/1‑150 (MQ=255)
                                                                                                                                    acaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTGTAGGTGGGGTCAGGCGGTCAcg                  >  2:6370315/1‑150 (MQ=255)
                                                                                                                                     caAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcgc                 <  1:2745580/150‑1 (MQ=255)
                                                                                                                                     caAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcgc                 >  1:1982751/1‑150 (MQ=255)
                                                                                                                                        caaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCggg              >  1:7237091/2‑150 (MQ=255)
                                                                                                                                        aaaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCggg              <  2:3958379/150‑1 (MQ=255)
                                                                                                                                         aaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc             >  2:1123758/1‑150 (MQ=255)
                                                                                                                                         aaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc             >  1:7381627/1‑150 (MQ=255)
                                                                                                                                         aaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc             >  1:5651814/1‑150 (MQ=255)
                                                                                                                                               ctctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCg       >  2:285772/1‑150 (MQ=255)
                                                                                                                                                 ctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGt     <  1:50263/150‑1 (MQ=255)
                                                                                                                                                 ctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGt     <  2:6637352/150‑1 (MQ=255)
                                                                                                                                                    gCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTGCa  >  2:4479025/1‑150 (MQ=255)
                                                                                                                                                    gCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTGCa  >  1:3432172/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTGCA  >  NC_002947/1501625‑1501922

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: