Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,314,654 T→A 7.0% intergenic (+30/‑132) 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,314,6540TA 7.0% 337.0 / 9.0 108intergenic (+30/‑132)PP_2912/hemBhypothetical protein/5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (44/56);  new base A (3/5);  total (47/61)
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

Only 100 of 108 total aligned reads displayed.

GGTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTC  >  NC_002947/3314507‑3314800
                                                                                                                                                   |                                                                                                                                                  
ggTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTAc                                                                                                                                                  <  2:3253162/150‑1 (MQ=255)
     tGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCcgcg                                                                                                                                             <  1:664666/150‑1 (MQ=255)
     tGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCcgcg                                                                                                                                             >  1:867883/1‑150 (MQ=255)
       cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa                                                                                                                                           <  2:2846240/150‑1 (MQ=255)
       cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGaa                                                                                                                                           <  2:4937400/150‑1 (MQ=255)
       cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGaa                                                                                                                                           <  1:3979426/150‑1 (MQ=255)
       cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGaa                                                                                                                                           <  2:2257076/150‑1 (MQ=255)
         tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTGGCGGGGTTAACCGCGAaga                                                                                                                                         >  1:908280/1‑150 (MQ=255)
         tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                         <  2:2379761/150‑1 (MQ=255)
          gagaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGTGGAGGTAACGCCGCCCTGGCTGCAGGAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAagag                                                                                                                                        <  1:3208499/150‑1 (MQ=255)
          gagaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAagag                                                                                                                                        <  2:1068665/150‑1 (MQ=255)
            gaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGGTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGc                                                                                                                                      <  2:4166447/150‑1 (MQ=255)
             agcGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAAGAGGcg                                                                                                                                     <  2:3747468/148‑2 (MQ=255)
               cGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCaa                                                                                                                                   <  2:2100558/150‑1 (MQ=255)
                  cgccgcCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGCAGAGGCcaccac                                                                                                                                >  1:3851798/1‑150 (MQ=255)
                   gccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAAcaca                                                                                                                               >  1:5392592/1‑150 (MQ=255)
                   gccgccCGGATGGACGACGCCCGTCGCCGTGCAGCACGCGCCGAGGCCGCGCATGGGCCGGCCACGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAAGAGGcggtacag                                                                                                                               <  2:5443203/150‑8 (MQ=38)
                    ccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAg                                                                                                                              <  2:1154250/150‑1 (MQ=255)
                    ccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAg                                                                                                                              <  2:2371953/150‑1 (MQ=255)
                     cgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAgg                                                                                                                             >  1:4799832/1‑150 (MQ=255)
                        ccGGATGGACGACGCTCGTCGCCGTGCATCACGCGCCGAGGCCGAGAATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGTCTGCAGAAGTTTATTCCTGTACCGGCTCTTCGCCGGGGTACCCCCGGACTCGCCAAAACAGGCaa                                                                                                                          >  1:532559/1‑150 (MQ=255)
                        ccGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCAGGGGCCGGCCAAGGCAGTGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCaa                                                                                                                          <  1:246322/150‑1 (MQ=255)
                         cGGATGGCCGACGCTCGCCGCCGTGCAGCCCGCGCCGAGGCCGAGCATGGGCCGGCCACGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTCCCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCaaa                                                                                                                         <  2:5000259/150‑1 (MQ=255)
                          ggATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAAc                                                                                                                        <  1:3263406/150‑1 (MQ=255)
                          ggATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAAc                                                                                                                        <  1:3263861/150‑1 (MQ=255)
                                     gCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTg                                                                                                             <  1:1683225/150‑1 (MQ=255)
                                      cTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTTTTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGTAGTGGCCAACTCAGGCAATCCACGGCCGTTGc                                                                                                            <  1:4990648/150‑1 (MQ=255)
                                            gCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGa                                                                                                      <  2:2258914/150‑1 (MQ=255)
                                            gCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGa                                                                                                      >  2:5439572/1‑150 (MQ=255)
                                                  cagctCGCGCGGAGGCGGTGCTTGGGCCGGCCTATGCATTGGTGTTAACGCCGGCCTGGGTGCAGTTGTTTTTTCCTGTACCGCCTCTTGGCGTGGTTACCCGCGAAGAGGCCATCACAGGCATACCTCGTCCGTTGCCATTGACGCAAg                                                                                                <  2:532559/150‑1 (MQ=255)
                                                     caCGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGaca                                                                                             >  2:3665216/1‑150 (MQ=255)
                                                        gcgcCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGGAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGacacct                                                                                          >  1:3184171/1‑148 (MQ=255)
                                                        gcgcCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACAt                                                                                          >  2:3925088/1‑150 (MQ=255)
                                                        gcgcCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGGTACCCGCGGAGAGGCCAACAAAGGCCAACCACGGCCGTTGCCAATGACGCAAGacacca                                                                                          >  1:3183920/1‑148 (MQ=255)
                                                          gccgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGAc                                                                                                     <  1:4788361/137‑1 (MQ=255)
                                                          gccgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGAc                                                                                                     >  2:4788361/1‑137 (MQ=255)
                                                          gccgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtc                                                                                        >  2:491474/1‑150 (MQ=255)
                                                            cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                      <  1:893731/150‑1 (MQ=255)
                                                            cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                      <  2:1301485/150‑1 (MQ=255)
                                                               ggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTc                                                                                   <  2:3054227/150‑1 (MQ=255)
                                                                  cgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCCGTGACGCAACACACATTCTCGTCAtt                                                                                >  2:1951514/1‑150 (MQ=255)
                                                                  cgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCAATGACGCAAAACACATTATCGGCAtt                                                                                >  1:4403959/1‑150 (MQ=255)
                                                                   gagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGGCAttt                                                                               >  1:5396173/1‑150 (MQ=255)
                                                                   gagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGGTACCCGCGAAGAGGCCACCACAGGCCAACCCCGGCCGGTGCCCTTGACGCAAGACACAGTCTCGTcctgt                                                                               >  2:4292844/1‑146 (MQ=255)
                                                                      cATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGc                                                                            <  1:3665216/150‑1 (MQ=255)
                                                                         agctcgcccAAGGCAGAGGAGGAAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCACATCGCGGGGTTACCCGCGAAGAGGCCAACACACGCAAACCACGGCCGTTGCCATTGACCCAACACACATTCTCCTCATTTAGCCtt                                                                         <  2:5171233/143‑1 (MQ=255)
                                                                           gccggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCtttt                                                                       <  2:2478855/150‑1 (MQ=255)
                                                                            ccggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCCTTGCCCATGACGCAAGACACATTCTCGTCATTTAGCCGTTTc                                                                      >  1:4837486/1‑150 (MQ=255)
                                                                             cggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCa                                                                     <  1:1664036/150‑1 (MQ=255)
                                                                             cggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCa                                                                     <  2:635772/150‑1 (MQ=255)
                                                                              ggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTCATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCAc                                                                    <  2:5153017/150‑1 (MQ=255)
                                                                               gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACt                                                                   <  2:730419/150‑1 (MQ=255)
                                                                                 cAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTAGCGGGGTTACCCGCGAATATGCCAAAAAAGGCAAAACACGGCCGTTTCCATTGAAGCAAGACAAATTCTCGTCATTTAGACTTTTCACTgg                                                                 >  2:188540/1‑150 (MQ=255)
                                                                                      cAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTTGCCTTTTCACGGgcgcgc                                                            >  2:42786/1‑150 (MQ=255)
                                                                                           gaggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAg                                                       >  2:4811548/1‑150 (MQ=255)
                                                                                           gaggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAg                                                       <  1:1998995/150‑1 (MQ=255)
                                                                                             ggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGc                                                     >  2:833426/1‑150 (MQ=255)
                                                                                                  cgccgcCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg                                                <  1:1042697/150‑1 (MQ=255)
                                                                                                  cgccgcCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg                                                <  2:2636077/150‑1 (MQ=255)
                                                                                                  cgccgcCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg                                                <  2:2259959/150‑1 (MQ=255)
                                                                                                   gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAgg                                               >  2:2962904/1‑150 (MQ=255)
                                                                                                    ccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGt                                              <  1:2290718/150‑1 (MQ=255)
                                                                                                              tGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTAACCGCGAAGAGGCCAAAACACGCAAACCACCGCCGGTGCCCTTGACCCAAGACACAAACTCGGCCATTTGCCTTTTAACTTGCCCGCGACAGGCATAAGGTATACGCCGtt                                    >  1:2166143/1‑150 (MQ=255)
                                                                                                               gCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTc                                   <  1:2985563/150‑1 (MQ=255)
                                                                                                                      tttATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGGAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTcc                            >  2:866981/1‑150 (MQ=255)
                                                                                                                      tttATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTcc                            <  1:3184012/150‑1 (MQ=255)
                                                                                                                        tATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCt                          <  2:4799832/150‑1 (MQ=255)
                                                                                                                              tGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCt                    >  2:2017559/1‑150 (MQ=255)
                                                                                                                              tGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCt                    <  2:1650025/150‑1 (MQ=255)
                                                                                                                               gTACCGCCTCTTCGCGGGGTTACCCCCGAAGAGGCCCACACAGGCAAACCACGGCCCTTGGCATTTAACCAAGAAACATTATCGTGAGTTAGCCTTTTGACGGGCGCGCGACAGGCATAAGGTATAAGCCCTTAAAAATACCAGGAGCTg                   >  1:2317909/1‑150 (MQ=255)
                                                                                                                                tACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGt                  <  2:4378929/150‑1 (MQ=255)
                                                                                                                                 aCCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTc                 <  2:867883/150‑1 (MQ=255)
                                                                                                                                 aCCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCAGAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTc                 <  2:3302599/150‑1 (MQ=255)
                                                                                                                                   cGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCt                    >  1:5130277/1‑145 (MQ=255)
                                                                                                                                   cGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCt                    <  2:5130277/145‑1 (MQ=255)
                                                                                                                                   cGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCAGTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAAATCCCTGGAGCTGTCAt               >  2:243848/1‑150 (MQ=255)
                                                                                                                                    gCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATg              >  1:3915508/1‑150 (MQ=255)
                                                                                                                                     ccTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACAGATTCTGGTCATTTTGCCTTTTGAGTGGGGCGCGTGTGGCATTAGGTATACGCGGTTGAAATTCCGTGGAGCTGTGATGt             >  1:18957/1‑150 (MQ=255)
                                                                                                                                     ccTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGt             >  1:18945/1‑150 (MQ=255)
                                                                                                                                      ctctTCGCGGGTAACCCCGCGAAGAGG‑CGGTACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACGGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTcc           >  1:1559134/1‑150 (MQ=38)
                                                                                                                                       tctTCGCGGGTAACCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGa                       >  1:4320563/1‑138 (MQ=255)
                                                                                                                                       tctTCGCGGGTAACCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGa                       <  2:4320563/138‑1 (MQ=255)
                                                                                                                                       tctTCGCGGGTAACCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTcc           >  2:3568508/1‑150 (MQ=255)
                                                                                                                                       tctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCtttt                                                                       >  2:1009934/1‑90 (MQ=255)
                                                                                                                                       tctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCtttt                                                                       <  1:1009934/90‑1 (MQ=255)
                                                                                                                                       tctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTcc           >  2:329664/1‑150 (MQ=255)
                                                                                                                                       tctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTcc           >  2:110490/1‑150 (MQ=255)
                                                                                                                                       tctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTcc           >  2:3059911/1‑150 (MQ=255)
                                                                                                                                       tctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTcc           >  1:4307212/1‑150 (MQ=255)
                                                                                                                                        ctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTcca          <  2:2487503/150‑1 (MQ=255)
                                                                                                                                         ttCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa         <  2:908280/150‑1 (MQ=255)
                                                                                                                                         ttCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa         >  1:3177519/1‑150 (MQ=255)
                                                                                                                                         ttCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa         >  2:3233450/1‑150 (MQ=255)
                                                                                                                                          tCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaac        >  1:4681746/1‑150 (MQ=255)
                                                                                                                                          tCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaac        <  1:4811548/150‑1 (MQ=255)
                                                                                                                                          tCGCGGGCTTACCCGCGAAGAGCCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAACACACATACTCCTCATTTACCCTATCCACTCGCGCGCGACACGCATCACGTATACGCCGTTCAAATTCCCAGGAGCTCACATCTccaac        <  2:4915722/150‑1 (MQ=255)
                                                                                                                                           cgcgGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaacc       <  2:3851798/150‑1 (MQ=255)
                                                                                                                                             cgGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaaccaa     <  1:2438730/150‑1 (MQ=255)
                                                                                                                                                ggTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTc  <  2:3374648/150‑1 (MQ=255)
                                                                                                                                                ggTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTc  >  2:2533431/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTC  >  NC_002947/3314507‑3314800

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