Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,254,706 C→G 6.8% intergenic (+24/+47) bvgA → / ← PP_1091 virulence factors transcription regulator BvgA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,254,7060CG 6.8% 2273.5 / 110.6 676intergenic (+24/+47)bvgA/PP_1091virulence factors transcription regulator BvgA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (353/277);  new base G (26/20);  total (379/297)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 708 total aligned reads displayed.

ATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGATTTCA  >  NC_002947/1254557‑1254855
                                                                                                                                                     |                                                                                                                                                     
  tGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTccagcca                                                                                                                                                     <  1:4797977/150‑1 (MQ=255)
   gCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTccagccag                                                                                                                                                    <  2:57427/150‑1 (MQ=255)
     cGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCagga                                                                                                                                                  >  1:3607508/1‑150 (MQ=255)
          gCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCtt                                                                                                                                             <  1:1572252/150‑1 (MQ=255)
           cATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCttt                                                                                                                                            <  2:792055/150‑1 (MQ=255)
           cATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCttt                                                                                                                                            <  1:3742823/150‑1 (MQ=255)
            aTGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTg                                                                                                                                           >  2:3064874/1‑150 (MQ=255)
            aTGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTg                                                                                                                                           >  1:7941651/1‑150 (MQ=255)
             tGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCTGGCTGGAGGCttcaa                                                                                                                                          <  2:4680140/150‑4 (MQ=255)
              gTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCttcaag                                                                                                                                         <  1:1514312/150‑5 (MQ=255)
               tACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCACTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCa                                                                                                                                        <  1:5689452/150‑1 (MQ=255)
                aCCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAg                                                                                                                                       <  2:7453874/150‑1 (MQ=255)
                     aGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTa                                                                                                                                  >  1:5044259/1‑150 (MQ=255)
                     aGTAACAAGACGGTCAGTACCCACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTa                                                                                                                                  <  2:3515885/150‑1 (MQ=255)
                      gTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTAt                                                                                                                                 <  1:1009179/150‑1 (MQ=255)
                       tAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATc                                                                                                                                >  1:5875256/1‑150 (MQ=255)
                         aCAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCATTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCgg                                                                                                                              <  1:7143334/150‑1 (MQ=255)
                         aCAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCgg                                                                                                                              >  2:4824054/1‑150 (MQ=255)
                           aaGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTc                                                                                                                            <  1:2551323/150‑1 (MQ=255)
                           aaGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTc                                                                                                                            <  1:5347103/150‑1 (MQ=255)
                              aCGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCAt                                                                                                                         <  2:4132583/150‑1 (MQ=255)
                              aCGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCAt                                                                                                                         >  1:428492/1‑150 (MQ=255)
                              aCGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCAt                                                                                                                         <  2:2084530/150‑1 (MQ=255)
                              aCGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCAt                                                                                                                         >  2:2675814/1‑150 (MQ=255)
                                      tACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACgg                                                                                                                 <  2:3006985/150‑1 (MQ=255)
                                         cTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGAc                                                                                                              >  1:3307271/1‑150 (MQ=255)
                                           aCAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGAccc                                                                                                            <  2:1831952/150‑1 (MQ=255)
                                           aCAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGAccc                                                                                                            >  1:3255820/1‑150 (MQ=255)
                                              aGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGACAAGCCATGAAGCATACAGCCAGGAGGCTTTGACAGGGCTATCGGTCCATGACTACGGGACCCGCa                                                                                                         >  1:6785399/1‑150 (MQ=255)
                                                  aCGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTAc                                                                                                     >  1:728849/1‑150 (MQ=255)
                                                  aCGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTAc                                                                                                     >  2:5325063/1‑150 (MQ=255)
                                                  aCGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTAc                                                                                                     >  2:2537917/1‑150 (MQ=255)
                                                   cGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACa                                                                                                    <  2:1534155/150‑1 (MQ=255)
                                                    ggCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAg                                                                                                   <  1:6937212/150‑1 (MQ=255)
                                                     gctgctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAgg                                                                                                  <  1:5624396/150‑1 (MQ=255)
                                                        gctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCg                                                                                               <  1:3120427/150‑1 (MQ=255)
                                                          tgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTa                                                                                             >  2:6312779/1‑150 (MQ=255)
                                                           gctcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTaa                                                                                            <  1:3052659/150‑1 (MQ=255)
                                                             tcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATc                                                                                          >  2:923388/1‑150 (MQ=255)
                                                             tcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATc                                                                                          >  1:7646820/1‑150 (MQ=255)
                                                               aagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGa                                                                                        <  2:5505516/150‑1 (MQ=255)
                                                                agctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGaa                                                                                       <  2:7662187/150‑1 (MQ=255)
                                                                  ctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAAt                                                                                     >  2:5605113/1‑150 (MQ=255)
                                                                   tcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATc                                                                                    >  1:8743282/1‑150 (MQ=255)
                                                                      aTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGc                                                                                 <  1:2471709/150‑1 (MQ=255)
                                                                      aTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGc                                                                                 >  2:5092125/1‑150 (MQ=255)
                                                                      aTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGc                                                                                 >  1:863987/1‑150 (MQ=255)
                                                                       tGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGcc                                                                                >  2:4507053/1‑150 (MQ=255)
                                                                       tGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGcc                                                                                >  1:3960061/1‑150 (MQ=255)
                                                                         cACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCaa                                                                              <  2:6502125/150‑1 (MQ=255)
                                                                              tCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCt                                                                         <  2:6120177/150‑1 (MQ=255)
                                                                                aTTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTc                                                                       >  2:735467/1‑150 (MQ=255)
                                                                                       gACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGgc                                                                >  2:4007202/1‑150 (MQ=255)
                                                                                         ccTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGgcgc                                                              >  2:8393741/1‑150 (MQ=255)
                                                                                          cTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGcc                                                             >  1:115198/1‑150 (MQ=255)
                                                                                           tGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCg                                                            <  2:5920851/150‑1 (MQ=255)
                                                                                            gATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGt                                                           <  2:8553313/150‑1 (MQ=255)
                                                                                               cGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttct                                                        <  1:938658/150‑1 (MQ=255)
                                                                                                gAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttctt                                                       >  1:2034022/1‑150 (MQ=255)
                                                                                                 aGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                                      <  2:8205643/150‑1 (MQ=255)
                                                                                                 aGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                                      >  1:4660663/1‑150 (MQ=255)
                                                                                                     ggCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGc                                                  <  1:2939776/150‑1 (MQ=255)
                                                                                                        gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtg                                               >  2:707349/1‑150 (MQ=255)
                                                                                                        gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtg                                               >  1:4076478/1‑150 (MQ=255)
                                                                                                        gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtg                                               >  1:3552661/1‑150 (MQ=255)
                                                                                                        gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtg                                               >  1:1814991/1‑150 (MQ=255)
                                                                                                        gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtg                                               >  2:8021600/1‑150 (MQ=255)
                                                                                                        gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtg                                               >  1:517991/1‑150 (MQ=255)
                                                                                                        gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtg                                               <  1:550172/150‑1 (MQ=255)
                                                                                                            cGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtgttgt                                           <  1:6457303/150‑1 (MQ=255)
                                                                                                               cAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGa                                        <  1:8255603/150‑1 (MQ=255)
                                                                                                               cAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGa                                        >  1:4438918/1‑150 (MQ=255)
                                                                                                                    tCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTg                                   <  2:7278613/150‑1 (MQ=255)
                                                                                                                    tCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTg                                   <  1:2259071/150‑1 (MQ=255)
                                                                                                                          cTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTAt                                                                                                                                 >  1:8434484/1‑50 (MQ=255)
                                                                                                                           tGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGtt                            <  1:8119818/150‑1 (MQ=255)
                                                                                                                           tGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATAGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGtt                            >  1:4822590/1‑150 (MQ=255)
                                                                                                                            gATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGttt                           <  1:6136793/150‑1 (MQ=255)
                                                                                                                             aTCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTg                          <  1:6618126/150‑1 (MQ=255)
                                                                                                                             aTCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTg                          <  2:8661847/150‑1 (MQ=255)
                                                                                                                              tCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGc                         >  2:1374550/1‑150 (MQ=255)
                                                                                                                               cGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCt                        <  2:4336345/150‑1 (MQ=255)
                                                                                                                                gCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTc                       >  1:1329293/1‑150 (MQ=255)
                                                                                                                                 ccAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCa                      <  2:4943532/150‑1 (MQ=255)
                                                                                                                                 ccAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCa                      <  2:3945323/150‑1 (MQ=255)
                                                                                                                                  cAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAg                     >  1:2789526/1‑150 (MQ=255)
                                                                                                                                  cAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAg                     <  2:6869100/150‑1 (MQ=255)
                                                                                                                                   aaGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAgg                    <  2:722058/150‑1 (MQ=255)
                                                                                                                                      tggcaaaGCCTCCTGGCTGGAGGCTTCGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTg                 >  2:3559608/6‑150 (MQ=255)
                                                                                                                                      ccTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTg                 <  2:2296222/150‑1 (MQ=255)
                                                                                                                                        ttGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGt               <  2:3374438/150‑1 (MQ=255)
                                                                                                                                         tGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTc              >  2:318739/1‑150 (MQ=255)
                                                                                                                                           aaGCCTCCTGGCTGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCtt            >  2:3756229/1‑150 (MQ=255)
                                                                                                                                            aGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCttt           <  1:7366115/150‑1 (MQ=255)
                                                                                                                                                tCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGTCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGa       >  2:3705660/1‑150 (MQ=255)
                                                                                                                                                tCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttct                                                        >  2:6687211/1‑101 (MQ=255)
                                                                                                                                                tCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGa       >  1:6283471/1‑150 (MQ=255)
                                                                                                                                                tCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGa       >  1:4533194/1‑150 (MQ=255)
                                                                                                                                                tCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGa       >  2:1240638/1‑150 (MQ=255)
                                                                                                                                                   agccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGAttt    >  2:2808501/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGATTTCA  >  NC_002947/1254557‑1254855

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: