Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,314,655 A→C 11.6% intergenic (+31/‑131) 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,6550AC11.6% 340.4 / 34.6 121intergenic (+31/‑131)PP_2912/hemBhypothetical protein/5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (50/56);  new base C (3/11);  total (54/67)
Fisher's exact test for biased strand distribution p-value = 8.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

Only 100 of 122 total aligned reads displayed.

TGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTCCCCT  >  NC_002947/3314509‑3314804
                                                                                                                                                  |                                                                                                                                                     
tGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTAccc                                                                                                                                                    <  2:2699178/150‑1 (MQ=255)
     cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa                                                                                                                                               <  1:2676032/150‑1 (MQ=255)
     cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa                                                                                                                                               <  1:2593/150‑1 (MQ=255)
       tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                             <  1:497115/150‑1 (MQ=255)
       tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                             <  2:1050675/150‑1 (MQ=255)
       tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                             <  1:755328/150‑1 (MQ=255)
       tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                             <  1:2544481/150‑1 (MQ=255)
       tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                             <  1:496857/150‑1 (MQ=255)
       tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAaga                                                                                                                                             <  2:3628810/150‑1 (MQ=255)
       tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAaga                                                                                                                                             <  2:1707974/150‑1 (MQ=255)
       tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAaga                                                                                                                                             <  2:156406/150‑1 (MQ=255)
         agaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTATTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAgg                                                                                                                                           <  1:3919805/150‑1 (MQ=255)
         agaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAAGAgg                                                                                                                                           <  2:1083102/150‑1 (MQ=255)
          gaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAAGAGGc                                                                                                                                          <  1:2145143/150‑1 (MQ=255)
          gaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGc                                                                                                                                          <  2:4204008/150‑1 (MQ=255)
            tCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCa                                                                                                                                        <  1:3716083/150‑1 (MQ=255)
            tCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCa                                                                                                                                        <  1:3716143/150‑1 (MQ=255)
              gACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCCGCACGCGCCGAGGCCGAGCATGGGCCGGCCATGGCAGTGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCGTCGCGGGGTTACCCGCGAAGAGGCCAac                                                                                                                                      >  1:1191909/1‑150 (MQ=255)
                 gccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGATGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAAcaca                                                                                                                                   >  2:583242/1‑150 (MQ=255)
                 gccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAAcaca                                                                                                                                   >  2:3302596/1‑150 (MQ=255)
                      ccGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCaa                                                                                                                              <  2:2038877/150‑1 (MQ=255)
                          aTGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCATGGCTTAGGAGGTAACGCCGCCCTGGCTGCTGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCa                                                                                                                          <  1:418953/150‑1 (MQ=255)
                          aTGGACGACGCGCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCa                                                                                                                          <  1:3860058/150‑1 (MQ=255)
                            ggACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTAACCGCGAAGAGGCCAACACAGGCAAACCCCg                                                                                                                        >  1:4303564/1‑150 (MQ=255)
                             gacgacGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCACACCACgg                                                                                                                       >  1:3783583/1‑150 (MQ=255)
                              acgacgCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGc                                                                                                                      >  2:2068791/1‑150 (MQ=255)
                                  cgCTCGGCGCCGTGCAGCACGCGCCGAGGCCGAGCTTGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGtt                                                                                                                  <  1:423838/150‑1 (MQ=255)
                                   gCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTATTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTg                                                                                                                 <  1:3765827/150‑1 (MQ=255)
                                    cTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTCGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGc                                                                                                                <  1:845757/150‑1 (MQ=255)
                                     tcgtcgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGcc                                                                                                               <  1:367482/150‑1 (MQ=255)
                                         cgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGGTGCCATTg                                                                                                           >  2:2017198/1‑150 (MQ=255)
                                           ccGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGAc                                                                                                         <  2:3533306/150‑1 (MQ=255)
                                            cGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGTGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACg                                                                                                        <  1:116479/150‑1 (MQ=255)
                                               gcagcaCGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCaa                                                                                                     <  1:4121893/150‑1 (MQ=255)
                                                 agcaCGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGa                                                                                                   >  2:340806/1‑150 (MQ=255)
                                                 agcaCGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGa                                                                                                   <  2:943117/150‑1 (MQ=255)
                                                    ccgcgaagaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGacac                                                                                                >  2:2620605/8‑150 (MQ=255)
                                                    ccgcgaagaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGacac                                                                                                >  2:2620617/8‑150 (MQ=255)
                                                    aCGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGCAGAGGCCAACACAGGCAAACCACGGCCGTTGCCCTTGACGCAAGacac                                                                                                >  2:3882784/1‑150 (MQ=255)
                                                    aCGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCCACACAGGCAAACCCCGGCCGTTTCCATTTACGCAAGacac                                                                                                >  1:482128/1‑150 (MQ=255)
                                                       cgcCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTCGCCATTGACGCAAGACACAtt                                                                                             <  1:3120820/150‑1 (MQ=255)
                                                         ccgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtct                                                                                           <  2:3620508/150‑1 (MQ=255)
                                                          cgaggccgagCTTGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                          <  1:1996899/150‑1 (MQ=255)
                                                          cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                          <  1:3384813/150‑1 (MQ=255)
                                                          cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                          <  1:3381823/150‑1 (MQ=255)
                                                           gaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCGTCGCGGGGTTACCCGCGAAGAGGCCCACACAGGCAAACCCCGGCCGTTGCCATTTACGCAAGACACATTCTCg                                                                                         >  1:827495/1‑150 (MQ=255)
                                                              gccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAGGACACATTCTCGTCa                                                                                      >  2:2407540/1‑150 (MQ=255)
                                                                  agCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTa                                                                                  <  1:1168889/150‑1 (MQ=255)
                                                                   gCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGGCATTTAg                                                                                 >  1:707453/1‑150 (MQ=255)
                                                                    cATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGc                                                                                <  1:4406718/150‑1 (MQ=255)
                                                                        ggccggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCttt                                                                            <  2:91485/150‑1 (MQ=255)
                                                                            ggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCAc                                                                        <  2:3570527/150‑1 (MQ=255)
                                                                             gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACt                                                                       <  2:3648214/150‑1 (MQ=255)
                                                                                aaGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTTACTGgc                                                                    <  2:1143274/150‑1 (MQ=255)
                                                                                aaGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGgc                                                                    <  2:1143331/150‑1 (MQ=255)
                                                                                aaGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGgc                                                                    <  2:1143148/150‑1 (MQ=255)
                                                                                aaGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGgc                                                                    >  2:985419/1‑150 (MQ=255)
                                                                                aaGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGGTTCCCGCGAAGAGGCCCACACAGGCCAACCACCGCCGTTGCCAATGACGCAAGACACAGTCGCCGGATTTGGCCGTTTGACAGgc                                                                    >  1:3027231/1‑150 (MQ=255)
                                                                                   gCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGgcgcg                                                                 <  2:284858/150‑1 (MQ=255)
                                                                                     agagGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCCCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGGTGCCAGTGACGCAAGACACAGTCGCGACAATTAGCCTTTTCACAGGcgcgcg                                                               >  2:3250187/1‑150 (MQ=255)
                                                                                       aggaggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGAc                                                             <  2:3501336/150‑1 (MQ=255)
                                                                                        ggaggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTATTCGCGGGGTTACCCGCGGAGAGGCCAAAACAGGCAAACCACGGCCGGTGCCCTTGACGCAAGACCCAATCTCGTCATTTAGCCTTTTCACTGGCGCGCGACa                                                            >  2:798224/1‑150 (MQ=255)
                                                                                           ggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGc                                                         <  1:3606720/150‑1 (MQ=255)
                                                                                              aaCGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATa                                                      >  2:470999/1‑150 (MQ=255)
                                                                                               aCGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGTCGCAAGACACATTCTCGTCCTTTAGCCTTTTCACTGGCGCGCGCCAGGCATaa                                                     >  2:3939615/1‑150 (MQ=255)
                                                                                                cgccgcCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg                                                    <  2:3215959/150‑1 (MQ=255)
                                                                                                 gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGTGGCCAACACAGGCAAACCACGGCCGTTGCCAGTGACGCAAGACACCTTCTCGTCAGTGAGCCTTTGCACGGGCGCGCGACAGGccgcagg                                                   >  1:3697278/1‑144 (MQ=255)
                                                                                                 gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTAGCGGGGTTACCCCCGAAGAGGCCAACACAGGCAAACCACGGCCGGTGCCATTGACGCAAGAAACCGTCGCGGCAGTTAGCCTTTTCACGGGCGCGCGACAGGccggagg                                                   >  1:731721/1‑144 (MQ=255)
                                                                                                  ccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGt                                                  <  2:2787910/150‑1 (MQ=255)
                                                                                                  ccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGt                                                  <  2:929499/150‑1 (MQ=255)
                                                                                                   cgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGta                                                 <  1:710713/150‑1 (MQ=255)
                                                                                                   cgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGta                                                 <  1:1400080/150‑1 (MQ=255)
                                                                                                             gCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTc                                       <  1:1067807/150‑1 (MQ=255)
                                                                                                               agtagtTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCCAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCaa                                     >  2:2830290/1‑150 (MQ=255)
                                                                                                                      tATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTTCCATTGTCGCAACAAACATTCTCCTCATTTATCCGTTTAACAGGCGCGCGACAGGCAAAAGGTATACGCCGTTCAAATTCCCt                              >  2:1128952/1‑150 (MQ=255)
                                                                                                                       aTTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACCTTCTCGTCATTTAGCCATTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTg                             >  2:3500216/1‑150 (MQ=255)
                                                                                                                       aTTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTg                             <  2:3869115/150‑1 (MQ=255)
                                                                                                                         tCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGa                           <  1:4368799/150‑1 (MQ=255)
                                                                                                                         tCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGa                           <  2:2762728/150‑1 (MQ=255)
                                                                                                                           cTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGc                         <  1:1976256/150‑1 (MQ=255)
                                                                                                                            tGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGGCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCt                        >  2:3844416/1‑150 (MQ=255)
                                                                                                                             gTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCAGAAGGTATACGCCGGTCAAATTACCTGGAGCTg                       >  1:706468/1‑150 (MQ=255)
                                                                                                                              tACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGt                      >  2:2232107/1‑150 (MQ=255)
                                                                                                                                ccGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCa                    >  2:1887861/1‑150 (MQ=255)
                                                                                                                                ccGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCa                    >  2:2750421/1‑150 (MQ=255)
                                                                                                                                 cGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCAt                   <  2:1689257/150‑1 (MQ=255)
                                                                                                                                 cGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCAt                   >  1:906816/1‑150 (MQ=255)
                                                                                                                                 cGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCAt                   >  2:3178018/1‑150 (MQ=255)
                                                                                                                                 cGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTACACGCCGTTCAAATTCCCTGGAGCTGTCAt                   <  2:4272089/150‑1 (MQ=255)
                                                                                                                                   ccTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGt                 >  2:2974779/1‑150 (MQ=255)
                                                                                                                                    ctctTCGCGGGTAACCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTc                >  1:3700730/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTc                                 <  2:3389355/132‑1 (MQ=255)
                                                                                                                                     tctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTcc               >  1:1857118/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa             <  1:3464257/150‑1 (MQ=255)
                                                                                                                                       ttCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa             >  1:359881/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGGCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa             >  1:1268946/1‑150 (MQ=255)
                                                                                                                                        tCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaac            <  1:1521599/150‑1 (MQ=255)
                                                                                                                                          gcgGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaacca          >  1:510939/1‑150 (MQ=255)
                                                                                                                                                 tACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTcccc   >  1:2601086/1‑150 (MQ=255)
                                                                                                                                                  aCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTCACCt  >  2:3126597/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTCCCCT  >  NC_002947/3314509‑3314804

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: