Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,249,757 C→G 10.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,249,7570CG10.8% 270.5 / 19.3 90intergenic (+24/+47)bvgA/PP_1091virulence factors transcription regulator BvgA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (46/33);  new base G (6/4);  total (53/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGATTT  >  NC_002947/1249608‑1249904
                                                                                                                                                     |                                                                                                                                                   
aTTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGAGCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTccagc                                                                                                                                                     >  1:792295/1‑150 (MQ=255)
   gCCGACAGCATGTACTTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTccagccag                                                                                                                                                  <  2:767193/150‑1 (MQ=255)
               tACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCa                                                                                                                                      >  2:1749698/1‑150 (MQ=255)
                 ccTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCTGGCTGGAGGCttcaa                                                                                                                                        <  1:2401197/146‑4 (MQ=255)
                 ccTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCTGGCTGGAGGCttcaa                                                                                                                                        >  2:2401197/1‑143 (MQ=255)
                 ccTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAgg                                                                                                                                    >  2:1441871/1‑150 (MQ=255)
                   tCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAAGGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGc                                                                                                                                  >  1:1359888/1‑150 (MQ=255)
                      gTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTAt                                                                                                                               <  1:1441871/150‑1 (MQ=255)
                         aCAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCgg                                                                                                                            >  2:469451/1‑150 (MQ=255)
                          cAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGACCTGCACGCAAAGCATGTACGTCTGCGTCCAGGTGTCGGTGCCAGGGCTATCGGt                                                                                                                           >  2:623972/1‑150 (MQ=255)
                             gACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGGGTGTAGGCTTTGCCAGGGCTATCGGTCCa                                                                                                                        >  1:603293/1‑150 (MQ=255)
                                      tACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACgg                                                                                                               <  1:2106320/150‑1 (MQ=255)
                                        ccTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCACGGCTATCGGTCCATGACTACGGGa                                                                                                             <  2:1622372/150‑1 (MQ=255)
                                         cTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGAc                                                                                                            >  2:1898806/1‑150 (MQ=255)
                                         cTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTAAAGCCTCCTGGCTGGAGGCtt                                                                                                                                           >  1:2078968/1‑119 (MQ=255)
                                         cTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTAAAGCCTCCTGGCTGGAGGCtt                                                                                                                                           <  2:2078968/119‑1 (MQ=255)
                                           aCAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGCAGTGGGCGGGCATTGATGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGAccc                                                                                                          >  1:1124337/1‑150 (MQ=255)
                                             aaGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGc                                                                                                                                             <  1:314081/113‑1 (MQ=255)
                                             aaGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGc                                                                                                                                             >  2:314081/1‑113 (MQ=255)
                                                tCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACt                                                                                                     >  2:1235023/1‑150 (MQ=255)
                                                tCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACt                                                                                                     >  2:2323271/1‑150 (MQ=255)
                                                  aCGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTccagccag                                                                                                                                                  <  2:2414170/103‑1 (MQ=255)
                                                  aCGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTccagccag                                                                                                                                                  >  1:2414170/1‑103 (MQ=255)
                                                  aCGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCttt                                                                                                                                          >  1:881091/1‑111 (MQ=255)
                                                  aCGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCttt                                                                                                                                          <  2:881091/111‑1 (MQ=255)
                                                  aCGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTAc                                                                                                   >  2:1336991/1‑150 (MQ=255)
                                                    ggCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAg                                                                                                 <  1:720292/150‑1 (MQ=255)
                                                     gctgctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAgg                                                                                                <  2:792295/150‑1 (MQ=255)
                                                     gctgctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAgg                                                                                                >  2:2000209/1‑150 (MQ=255)
                                                         ctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGt                                                                                            >  2:815457/1‑150 (MQ=255)
                                                           gctcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTaa                                                                                          <  1:2323271/150‑1 (MQ=255)
                                                             tcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGACCTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATc                                                                                        <  2:74385/150‑1 (MQ=255)
                                                             tcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCg                                                                                                                             >  1:2079758/1‑113 (MQ=255)
                                                             tcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCg                                                                                                                             <  2:2079758/113‑1 (MQ=255)
                                                             tcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATc                                                                                        >  2:84133/1‑150 (MQ=255)
                                                             tcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATc                                                                                        >  1:74385/1‑150 (MQ=255)
                                                               aagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGa                                                                                      >  1:1557156/1‑150 (MQ=255)
                                                                 gctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGaaa                                                                                    >  1:2032590/1‑150 (MQ=255)
                                                                  ctcaaTGCACGGTCATTGGTTGAGCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAAt                                                                                   >  2:1193213/1‑150 (MQ=255)
                                                                     aaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAgg                                                                                                                                    <  2:1887188/98‑1 (MQ=255)
                                                                     aaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAgg                                                                                                                                    >  1:1887188/1‑98 (MQ=255)
                                                                        gctcgttCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCAGCACTACAGGTCGTAATCGAAATCGGCCa                                                                             <  2:1124337/144‑1 (MQ=255)
                                                                          aCGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAAc                                                                           >  2:1589026/1‑150 (MQ=255)
                                                                              tCATTGGTTGACCTGACCGAGCTGGCGCAACGGCAAGGTCTGAACTGCTCGCCAAGCCTTGAAGCCTCCCGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCt                                                                       >  2:1227006/1‑150 (MQ=255)
                                                                                aTTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGTAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTc                                                                     >  2:1887422/1‑150 (MQ=255)
                                                                                 ttGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCt                                                                    <  1:815457/150‑1 (MQ=255)
                                                                                  tGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTg                                                                   <  1:2000209/150‑1 (MQ=255)
                                                                                    gTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCa                                                                 <  2:2047556/150‑1 (MQ=255)
                                                                                    gTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCa                                                                 >  1:2259679/1‑150 (MQ=255)
                                                                                    gTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCa                                                                 >  2:1420584/1‑150 (MQ=255)
                                                                                    gTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCa                                                                 >  1:2258729/1‑150 (MQ=255)
                                                                                          cTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGcc                                                           >  1:1672698/1‑150 (MQ=255)
                                                                                           tGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCg                                                          >  2:913154/1‑150 (MQ=255)
                                                                                                gAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttctt                                                     >  1:1361372/1‑150 (MQ=255)
                                                                                                 aGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                                    <  2:1359888/150‑1 (MQ=255)
                                                                                                 aGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                                    <  1:84133/150‑1 (MQ=255)
                                                                                                   cTGGCGCAACGGCAAGGTCTGGACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGTCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAGCTGCTTCTGCAGGCGCCGTTCTTCCa                                                  <  2:1210075/150‑1 (MQ=255)
                                                                                                   cTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                                  <  1:1235023/150‑1 (MQ=255)
                                                                                                        gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtg                                             >  1:132346/1‑150 (MQ=255)
                                                                                                            cGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtgttgt                                         <  2:603293/150‑1 (MQ=255)
                                                                                                             ggCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTc                                        <  2:2258729/150‑1 (MQ=255)
                                                                                                             ggCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTc                                        <  2:2259679/150‑1 (MQ=255)
                                                                                                                aaGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTAGGGGCCCCGGACTCCAGGTCGTTATGTAAATCGGCCAACGGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGAt                                     >  1:765108/1‑150 (MQ=255)
                                                                                                                aaGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGAt                                     >  1:2035873/1‑150 (MQ=255)
                                                                                                                    tCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTg                                 <  1:1014208/150‑1 (MQ=255)
                                                                                                                      tGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCg                               >  1:2229759/1‑150 (MQ=255)
                                                                                                                        aaCTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGt                                                                                               <  2:650932/84‑1 (MQ=255)
                                                                                                                        aaCTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGt                                                                                               >  1:650932/1‑84 (MQ=255)
                                                                                                                           tGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                                  <  2:822379/126‑1 (MQ=255)
                                                                                                                           tGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                                  >  1:822379/1‑126 (MQ=255)
                                                                                                                           tGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGtt                          >  2:1618496/1‑150 (MQ=255)
                                                                                                                           tGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGtt                          >  1:1132842/1‑150 (MQ=255)
                                                                                                                           tGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGtt                          >  1:331888/1‑150 (MQ=255)
                                                                                                                                gCCAAGCCTTGAAGCCTACAGACAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGCTCTTCCAGCGTGTGGTGGATGGTGCGGCGTTTGCTc                     >  2:1418273/1‑150 (MQ=255)
                                                                                                                                 ccAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCa                    >  2:992872/1‑150 (MQ=255)
                                                                                                                                  cAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAg                   >  2:1367289/1‑150 (MQ=255)
                                                                                                                                  cAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAg                   <  1:1618496/150‑1 (MQ=255)
                                                                                                                                   aaGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAgg                  <  1:1336991/150‑1 (MQ=255)
                                                                                                                                    cctggcaaaGCCTCCTGGCTGGAGGCTTCGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAgg                  >  1:2326453/8‑149 (MQ=255)
                                                                                                                                    cctggcaaaGCCTCCTGGCTGGAGGCTTCGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAgg                  <  2:2326453/142‑1 (MQ=255)
                                                                                                                                     gCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGtt                >  2:1932664/1‑150 (MQ=255)
                                                                                                                                      ccTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTg               <  1:1898806/150‑1 (MQ=255)
                                                                                                                                        ttGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTCCAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATAGAGCAACGTTTGCTCAGGTTGGt             >  1:121514/1‑150 (MQ=255)
                                                                                                                                        ttGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGt             <  1:469451/150‑1 (MQ=255)
                                                                                                                                         tGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTc            <  1:1193213/150‑1 (MQ=255)
                                                                                                                                         caaaGCCTCCTGGCTGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                                  >  1:2429067/3‑112 (MQ=255)
                                                                                                                                         caaaGCCTCCTGGCTGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                                  <  2:2429067/110‑1 (MQ=255)
                                                                                                                                           aaGCCTCCTGGCTGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCtt          >  1:1953395/1‑150 (MQ=255)
                                                                                                                                            aGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCttt         <  2:2032590/150‑1 (MQ=255)
                                                                                                                                               cTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCg      >  1:1486011/1‑150 (MQ=255)
                                                                                                                                                tCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAgg                                                                                                <  2:1084596/59‑1 (MQ=255)
                                                                                                                                                tCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAgg                                                                                                >  1:1084596/1‑59 (MQ=255)
                                                                                                                                                tCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCt                      >  2:101368/1‑133 (MQ=255)
                                                                                                                                                tCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCt                      <  1:101368/133‑1 (MQ=255)
                                                                                                                                                  cagccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGAtt   <  1:1589026/150‑1 (MQ=255)
                                                                                                                                                  cagccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGAtt   <  1:1716567/150‑1 (MQ=255)
                                                                                                                                                   agccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGGttt  >  1:139182/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTTGCGATTT  >  NC_002947/1249608‑1249904

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