Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,249,759 A→T 5.3% intergenic (+26/+45) 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,249,7590AT 5.3% 303.0 / 6.0 95intergenic (+26/+45)bvgA/PP_1091virulence factors transcription regulator BvgA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (41/48);  new base T (3/2);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 107 total aligned reads displayed.

CGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGATTTCAAC  >  NC_002947/1249613‑1249908
                                                                                                                                                   |                                                                                                                                                      
cGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTACTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTCAAGACTccgcccagta                                                                                                                                                     >  2:1118438/1‑148 (MQ=255)
 gACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCaggag                                                                                                                                                    <  1:1438845/150‑1 (MQ=255)
          tACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCa                                                                                                                                           >  1:2168047/1‑150 (MQ=255)
          tACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTAAAGCCTCCTGGCTGGAGGCttcaagg                                                                                                                                           <  1:2071689/150‑6 (MQ=255)
           aCCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAg                                                                                                                                          >  2:890032/1‑150 (MQ=255)
           aCCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTCGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAg                                                                                                                                          <  1:1031796/150‑1 (MQ=255)
             cTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAggg                                                                                                                                        <  2:1804472/150‑1 (MQ=255)
             cTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAggg                                                                                                                                        <  1:791814/150‑1 (MQ=255)
             cTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAggg                                                                                                                                        >  1:34706/1‑150 (MQ=255)
              tCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGc                                                                                                                                       >  1:1587086/1‑150 (MQ=255)
                     cAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGt                                                                                                                                <  1:214871/150‑1 (MQ=255)
                      aaGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCagg                                                                                                                                                      >  1:5469/1‑127 (MQ=255)
                      aaGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCagg                                                                                                                                                      <  2:5469/127‑1 (MQ=255)
                       aGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTcc                                                                                                                              >  2:677867/1‑150 (MQ=255)
                             tCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACt                                                                                                                        >  1:1905438/1‑150 (MQ=255)
                               aGTACCTCCCAGTCGCGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGAACTGCAGACGTGTGCTCTTCGGCAAGGGCGGACCTGATCGCCAAGCCTTGAAGCATCCATCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTAc                                                                                                                     <  1:2350110/150‑1 (MQ=255)
                                gTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCTGGCTGGAGGCttcaag                                                                                                                                            <  2:1257079/127‑5 (MQ=255)
                                gTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCTGGCTGGAGGCttcaag                                                                                                                                            >  1:1257079/1‑123 (MQ=255)
                                gTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACg                                                                                                                    <  1:842664/150‑1 (MQ=255)
                                  aCCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACggg                                                                                                                  <  2:1741587/150‑1 (MQ=255)
                                     tACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGAcc                                                                                                               >  2:1025960/1‑150 (MQ=255)
                                     tACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGAcc                                                                                                               >  1:13998/1‑150 (MQ=255)
                                      aCAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGAACTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGAccc                                                                                                              <  2:805503/150‑1 (MQ=255)
                                        aaGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGc                                                                                                            >  1:828712/1‑150 (MQ=255)
                                           tCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACt                                                                                                         >  1:989273/1‑150 (MQ=255)
                                                gctgctgccCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCACAGGCCACGTCATGACTGGTCGCCCACCATTGAAGGCCAAAGCCCCGAGGCTTTGCCATTGCTGTCGTTCCATGACTACGGGACCCGCACTACAgg                                                                                                    <  1:550447/150‑1 (MQ=255)
                                                 ctgctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGt                                                                                                   <  1:677867/150‑1 (MQ=255)
                                                 ctgctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGt                                                                                                   >  2:76556/1‑150 (MQ=255)
                                                  tgctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTc                                                                                                  <  1:1025960/150‑1 (MQ=255)
                                                      gctcaagctcaaTGCACTGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTaa                                                                                              <  2:1942637/150‑1 (MQ=255)
                                                      gctcaagctcaaTGCACGGCCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTaa                                                                                              <  2:1587086/150‑1 (MQ=255)
                                                       ctcaagctaaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAAt                                                                                             <  2:151757/150‑1 (MQ=255)
                                                        tcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATc                                                                                            >  1:692192/1‑150 (MQ=255)
                                                          aagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGa                                                                                          >  1:356069/1‑150 (MQ=255)
                                                          aagctcaaTGCACGGTCATGGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGCACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGa                                                                                          >  1:667186/1‑150 (MQ=255)
                                                             ctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAAt                                                                                       >  1:842468/1‑150 (MQ=255)
                                                              tcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATc                                                                                      <  1:929261/150‑1 (MQ=255)
                                                               gggggCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCg                                                                                     <  2:667186/146‑1 (MQ=255)
                                                                 aTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGCAGCCTCCAGCCAGGAGGCTTTGCAAGGGCTATCGGTCCATGACTACGGTACCCGCACTACAGGTCGTAATCGAAATCGGc                                                                                   <  2:884104/150‑1 (MQ=255)
                                                                   gCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCa                                                                                 >  1:1013426/1‑150 (MQ=255)
                                                                     aCGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGAAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAAc                                                                               >  1:2333151/1‑150 (MQ=255)
                                                                      cGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACt                                                                              <  1:2082318/150‑1 (MQ=255)
                                                                       ggTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAAGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTg                                                                             >  2:87586/1‑150 (MQ=255)
                                                                        gTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACgg                                                                                                                   >  1:1967611/1‑111 (MQ=255)
                                                                        gTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACgg                                                                                                                   <  2:1967611/111‑1 (MQ=255)
                                                                        gTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGc                                                                            <  2:2333151/150‑1 (MQ=255)
                                                                         tCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCt                                                                           >  2:2031411/1‑150 (MQ=255)
                                                                             tGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTg                                                                       <  2:2232066/150‑1 (MQ=255)
                                                                               gTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTcc                                                                                                                              <  1:1528595/94‑1 (MQ=255)
                                                                               gTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTcc                                                                                                                              >  2:1528595/1‑94 (MQ=255)
                                                                               gTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCa                                                                     <  1:2116173/150‑1 (MQ=255)
                                                                               gTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCa                                                                     >  2:2174696/1‑150 (MQ=255)
                                                                                       gATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCg                                                                                     <  1:2103488/126‑1 (MQ=255)
                                                                                       gATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCg                                                                                     >  2:2103488/1‑126 (MQ=255)
                                                                                       gATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGt                                                             <  2:1905438/150‑1 (MQ=255)
                                                                                           gAGCTGGCGCAACGGCAAGGTCTGACCTGCTGGAGCACGTTGGAAGGCGCTAGGCAGGGGGGTTTGCCGGGGCTAACGCCAGTGCCTGACCAAACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAAGCGCCGttctt                                                         >  1:1749770/1‑150 (MQ=255)
                                                                                            aGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                                        <  2:842468/150‑1 (MQ=255)
                                                                                            aGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                                        <  2:828712/150‑1 (MQ=255)
                                                                                              cTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                                      <  1:649588/150‑1 (MQ=255)
                                                                                              cTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                                      <  1:890032/150‑1 (MQ=255)
                                                                                                   cgcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGAACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAt                                                  <  1:1315112/150‑1 (MQ=255)
                                                                                                    gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtg                                                 >  1:935998/1‑150 (MQ=255)
                                                                                                      aaCGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                                        >  1:1452860/1‑141 (MQ=255)
                                                                                                      aaCGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                                        <  2:1452860/141‑1 (MQ=255)
                                                                                                       aCGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtgttg                                              <  1:1588150/150‑1 (MQ=255)
                                                                                                          gCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCg                                           >  1:1187584/1‑150 (MQ=255)
                                                                                                           cAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGTAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtgttgt                                             >  2:393504/1‑147 (MQ=255)
                                                                                                           cAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGTAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtgttgt                                             <  1:393504/147‑1 (MQ=255)
                                                                                                           cAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGa                                          >  1:2125153/1‑150 (MQ=255)
                                                                                                            aaGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGAt                                         <  2:935998/150‑1 (MQ=255)
                                                                                                             aGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATg                                        <  2:1017070/150‑1 (MQ=255)
                                                                                                                   gAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGc                                                                                   >  1:1671265/1‑101 (MQ=255)
                                                                                                                   gAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGc                                                                                   <  2:1671265/101‑1 (MQ=255)
                                                                                                                    aaCTGATCGCCAAGTCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGAc                                 <  2:692192/150‑1 (MQ=255)
                                                                                                                    aaCTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAAc                                                                               <  1:2010223/104‑1 (MQ=255)
                                                                                                                    aaCTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTg                                                                             <  1:2403546/106‑1 (MQ=255)
                                                                                                                    aaCTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTg                                                                             >  2:2403546/1‑106 (MQ=255)
                                                                                                                    aaCTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCGTGTTGTCGATGGTGCGAc                                 >  1:124656/1‑150 (MQ=255)
                                                                                                                       tGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                                        <  1:2076103/124‑1 (MQ=255)
                                                                                                                       tGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                                        >  2:2076103/1‑124 (MQ=255)
                                                                                                                        gATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGttt                             <  2:1501530/150‑1 (MQ=255)
                                                                                                                         aTCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGAccc                                                                                                              <  2:225029/68‑1 (MQ=255)
                                                                                                                         aTCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGAccc                                                                                                              >  1:225029/1‑68 (MQ=255)
                                                                                                                          tCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGc                           <  2:13998/150‑1 (MQ=255)
                                                                                                                            gCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGt                                                             >  1:898428/1‑114 (MQ=255)
                                                                                                                            gCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGt                                                             <  2:898428/114‑1 (MQ=255)
                                                                                                                              cAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAg                       <  2:2168047/150‑1 (MQ=255)
                                                                                                                                aGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGt                     <  2:1693342/150‑1 (MQ=255)
                                                                                                                                     tGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTc                <  2:1461993/150‑1 (MQ=255)
                                                                                                                                     caaaGCCTCCTGGCTGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTc                >  2:86124/3‑150 (MQ=255)
                                                                                                                                       aaGCCTCCTGGCTGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCTACGTTTGCTCAGGTTGGTCtt              >  1:2144478/1‑150 (MQ=255)
                                                                                                                                            tcaagccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGc                                                                                   <  1:170949/73‑1 (MQ=255)
                                                                                                                                            tCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGc                                                                                   >  2:170949/1‑76 (MQ=255)
                                                                                                                                              cagccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGAtt       >  2:2407440/1‑150 (MQ=255)
                                                                                                                                               agccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAgg                                                                                                    >  2:2389526/1‑56 (MQ=255)
                                                                                                                                               agccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAgg                                                                                                    <  1:2389526/56‑1 (MQ=255)
                                                                                                                                                 ccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAg                                                     >  1:1322663/1‑101 (MQ=255)
                                                                                                                                                 ccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAg                                                     <  2:1322663/101‑1 (MQ=255)
                                                                                                                                                  cagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGATTTCaa   >  2:769176/1‑150 (MQ=255)
                                                                                                                                                   aggaggCTTTGCCCGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGATTTCAAc  >  2:2394394/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGATTTCAAC  >  NC_002947/1249613‑1249908

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: