Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,312,730 G→T 7.2% intergenic (+28/‑134) PP_2912 → / → hemB hypothetical protein/5‑aminolevulinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,312,7300GT 7.2% 215.6 / 3.7 68intergenic (+28/‑134)PP_2912/hemBhypothetical protein/5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base T (3/2);  total (37/31)
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.95e-01

CCGGGGTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATT  >  NC_002947/3312581‑3312877
                                                                                                                                                      |                                                                                                                                                    
ccGGGGTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCgggg                                                                                                                                                      <  1:1586181/150‑1 (MQ=255)
   gggTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTa                                                                                                                                                   >  2:1926593/1‑150 (MQ=255)
           cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGaa                                                                                                                                          <  1:213432/150‑1 (MQ=255)
           cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGaa                                                                                                                                          <  1:1658721/150‑1 (MQ=255)
           cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGaa                                                                                                                                          >  2:722288/1‑150 (MQ=255)
            gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAag                                                                                                                                         <  2:2388704/150‑1 (MQ=255)
             tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                        <  1:705507/150‑1 (MQ=255)
                    gACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAgg                                                                                                                                      >  2:711342/1‑145 (MQ=255)
                    gACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAgg                                                                                                                                      <  1:711342/145‑1 (MQ=255)
                        ccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAg                                                                                                                             <  2:2056059/150‑1 (MQ=255)
                        ccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAg                                                                                                                             <  2:1262004/150‑1 (MQ=255)
                         cgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAAGAGGcggt                                                                                                                                  <  1:575004/144‑4 (MQ=255)
                         cgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAAGAGGcggt                                                                                                                                  >  2:575004/1‑141 (MQ=255)
                            ccGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                        >  1:969934/1‑135 (MQ=255)
                            ccGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                        <  2:969934/135‑1 (MQ=255)
                                  ggACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACg                                                                                                                   <  1:134286/150‑1 (MQ=255)
                                  ggACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACg                                                                                                                   <  1:1926593/150‑1 (MQ=255)
                                      gacgCTCGTCACCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCg                                                                                                               >  1:1336893/1‑150 (MQ=255)
                                         gCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCCACACAGGCAAACCACGGCCGTTg                                                                                                            >  1:125929/1‑150 (MQ=255)
                                         gCTCGTCGCCGTGCAGCACGCGCCGAAGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTg                                                                                                            >  1:1644269/1‑150 (MQ=255)
                                           tcgtcgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAAcaca                                                                                                                              <  2:1097568/130‑1 (MQ=255)
                                           tcgtcgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAAcaca                                                                                                                              >  1:1097568/1‑130 (MQ=255)
                                              ggcgcGTGCTGCACGGCGACGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGGGGCCAACACAGGCAAACCACGGCCGTTGCCAt                                                                                                        >  2:1160120/5‑150 (MQ=255)
                                                               gccgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtc                                                                                       >  2:1588423/1‑150 (MQ=255)
                                                                 cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                     <  1:2237302/150‑1 (MQ=255)
                                                                 cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAAGAgg                                                                                                                                      <  2:1376750/101‑1 (MQ=38)
                                                                 cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGATGCAGGAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAAGAgg                                                                                                                                      >  1:1376750/1‑101 (MQ=38)
                                                                  gaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCg                                                                                    >  1:345211/1‑150 (MQ=255)
                                                                   aggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGt                                                                                   >  1:792059/1‑150 (MQ=255)
                                                                    ggccgagCATGGGCCGGCCAAGGCAGAGGAGGGAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCACCACAGTCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTc                                                                                  <  1:1160120/150‑1 (MQ=255)
                                                                       cgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCAtt                                                                               <  1:1588423/150‑1 (MQ=255)
                                                                        gagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCAttt                                                                              <  2:1336893/150‑1 (MQ=255)
                                                                                gccggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCtttt                                                                      >  1:2236513/1‑150 (MQ=255)
                                                                                  cggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGAc                                                                                                    >  2:153335/1‑118 (MQ=255)
                                                                                  cggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGAc                                                                                                    <  1:153335/118‑1 (MQ=255)
                                                                                    gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACt                                                                  <  2:125929/150‑1 (MQ=255)
                                                                                     ccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTg                                                                 <  2:2236513/150‑1 (MQ=255)
                                                                                      cAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTa                                                                             >  2:2059876/1‑137 (MQ=255)
                                                                                      cAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTa                                                                             <  1:2059876/137‑1 (MQ=255)
                                                                                         ggCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCAc                                                                   <  2:2178358/144‑1 (MQ=255)
                                                                                         ggCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCAGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTAACGCAAGACACATTCTCGTCATTTAGCCTTTTCAc                                                                   >  1:2178358/1‑144 (MQ=255)
                                                                                                   gTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTc                                                                                  >  1:385094/1‑119 (MQ=255)
                                                                                                   gTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTc                                                                                  <  2:385094/119‑1 (MQ=255)
                                                                                                   gTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGAAAGACACATTCTCGTCATTTAGCCTTGTCACTGGCTCCCGACAGGCa                                                   >  2:590595/1‑150 (MQ=255)
                                                                                                       cgccgcCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg                                               <  2:345211/150‑1 (MQ=255)
                                                                                                         ccgccCCGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTTGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGt                                             <  1:20964/150‑1 (MQ=255)
                                                                                                           gccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGtat                                           >  1:1156743/1‑150 (MQ=255)
                                                                                                             ccTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATAc                                         >  1:520544/1‑150 (MQ=255)
                                                                                                              cTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACg                                        >  1:2279267/1‑150 (MQ=255)
                                                                                                                 gcgcgagtagtTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGAAACATAATAGTACATGAGCCTATTCCAGGGAGCTCGACGGGCATAGGGTGTAGGCCg                                     >  1:1667276/6‑150 (MQ=255)
                                                                                                                   tGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTa                                                                             <  1:256869/108‑1 (MQ=255)
                                                                                                                   tGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTa                                                                             >  2:256869/1‑108 (MQ=255)
                                                                                                                     cAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCa                                 >  1:351969/1‑150 (MQ=255)
                                                                                                                        tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAAt                              >  1:960919/1‑150 (MQ=255)
                                                                                                                        tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAAAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAAt                              >  2:1677822/1‑150 (MQ=255)
                                                                                                                           tttATTCCTGTACCGCCTCTTCGCGGGGTTGCCCGCGAAGAGGCCAACACAGGCAAACCACGGTCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACa                                                       >  2:1514984/1‑122 (MQ=255)
                                                                                                                           tttATTCCTGTACCGCCTCTTCGCGGGGTTGCCCGCGAAGAGGCCAACACAGGCAAACCACGGTCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACa                                                       <  1:1514984/122‑1 (MQ=255)
                                                                                                                              aTTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTg                        <  2:960919/150‑1 (MQ=255)
                                                                                                                                  cTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGc                    >  2:1443550/1‑150 (MQ=255)
                                                                                                                                  cTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGc                    >  2:2062492/1‑150 (MQ=255)
                                                                                                                                   tGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCt                   <  2:1284235/150‑1 (MQ=255)
                                                                                                                                     tACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGt                 <  1:722288/150‑1 (MQ=255)
                                                                                                                                       ccGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCa               >  1:1522431/1‑150 (MQ=255)
                                                                                                                                           ctctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTc           >  1:1169475/1‑150 (MQ=255)
                                                                                                                                           ctctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTc           <  2:2277379/150‑1 (MQ=255)
                                                                                                                                              ttCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa        >  1:1349703/1‑150 (MQ=255)
                                                                                                                                                 gcgGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaacca     >  2:218769/1‑150 (MQ=255)
                                                                                                                                                 gcgGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaacca     >  2:1284847/1‑150 (MQ=255)
                                                                                                                                                 gcgGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaacca     >  2:1663484/1‑150 (MQ=255)
                                                                                                                                                    gggTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAAtt  <  2:792059/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCGGGGTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATT  >  NC_002947/3312581‑3312877

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: