Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,781,365 T→G 16.2% intergenic (+55/‑88) PP_3338 → / → PP_3339 cytochrome bd‑type quinol oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,781,3650TG16.2% 395.8 / 13.8 187intergenic (+55/‑88)PP_3338/PP_3339cytochrome bd‑type quinol oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (54/97);  new base G (21/15);  total (75/112)
Fisher's exact test for biased strand distribution p-value = 2.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

CGTCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCT  >  NC_002947/3781216‑3781514
                                                                                                                                                     |                                                                                                                                                      
cgTCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCgtggt                                                                                                                                                        <  1:3673932/150‑1 (MQ=255)
  tCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTgc                                                                                                                                                      <  1:8421420/150‑1 (MQ=255)
  tCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTgc                                                                                                                                                      <  1:3988450/150‑1 (MQ=255)
      tCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGccg                                                                                                                                                  >  2:3575885/1‑149 (MQ=255)
       ttggctCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCACTCGCCGTGGTGccccag                                                                                                                                                 >  1:5488331/6‑148 (MQ=255)
       cAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCtt                                                                                                                                                 <  2:5684584/150‑2 (MQ=255)
        aGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGc                                                                                                                                                <  2:4799905/150‑1 (MQ=255)
          cATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGa                                                                                                                                              <  1:124974/150‑1 (MQ=255)
          cATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGa                                                                                                                                              >  1:4971630/1‑150 (MQ=255)
           aTCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCgcc                                                                                                                                             >  1:2925838/1‑148 (MQ=255)
            tCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACa                                                                                                                                            >  2:2312397/1‑150 (MQ=255)
            tCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGAcc                                                                                                                                            >  2:1181764/1‑149 (MQ=255)
            tCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGAcc                                                                                                                                            >  1:5914470/1‑149 (MQ=255)
              tGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGGCTGCGAaaaa                                                                                                                                          >  1:5694180/1‑150 (MQ=255)
                 gAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAg                                                                                                                                       <  1:6442104/150‑1 (MQ=255)
                  aaGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAgg                                                                                                                                      <  1:6927842/150‑1 (MQ=255)
                  aaGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAgg                                                                                                                                      <  1:2799885/150‑1 (MQ=255)
                   aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                     <  2:3124499/150‑1 (MQ=255)
                   aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                     <  1:1675109/150‑1 (MQ=255)
                     ccGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTCGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcgc                                                                                                                                   <  2:6816463/150‑1 (MQ=255)
                      cGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGca                                                                                                                                  <  2:8341214/150‑1 (MQ=255)
                       gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcac                                                                                                                                 <  2:2093202/150‑1 (MQ=255)
                         aGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGCCAAAAAGGCGcacca                                                                                                                              >  1:2023354/1‑150 (MQ=255)
                          gcgcCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcaccac                                                                                                                             <  2:3395209/150‑1 (MQ=255)
                           cgccgcCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACg                                                                                                                            >  2:8818870/1‑150 (MQ=255)
                             ccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACCAAAAGGCGCCCCACGGa                                                                                                                          >  1:8816601/1‑150 (MQ=255)
                              cgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCGGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGTGCCTGCGACAAAAAGGCGCACCACGGaa                                                                                                                         <  2:1397975/150‑1 (MQ=255)
                               gccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGacc                                                                                                                        >  2:6281509/1‑148 (MQ=255)
                                    agagCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAAGGGCACCCCCGAAcagtga                                                                                                                   >  1:6390679/1‑145 (MQ=255)
                                        ccAAGGCTTCATGCTGTTTGGAACACTTTTCATCCTGCCCTGTTGGGTTGGCTGATCCGACTTTTCATGTTGGCCGCCGGAATGCGCAGCAATTGAGCCCTCGCCTTGGTGCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAGc                                                                                                               <  2:2728425/150‑1 (MQ=255)
                                        ccAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGc                                                                                                               <  2:586291/150‑1 (MQ=255)
                                          aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGAACATGTAAAGccc                                                                                                             >  2:5602920/1‑150 (MQ=255)
                                           aGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGATTGTTCATGATGGCCGCCGGCATGCGCAGCCTTTGAGCCCCCGCGTTGGTGCGCCTGCGAAAAAAAGGCGAACCACGGAACTTGTACAGCCCg                                                                                                            <  2:5424211/150‑1 (MQ=255)
                                           aGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCCACAAAAAAGCGCCCCAACCAACAGGTACAGcccc                                                                                                            >  1:1361388/1‑149 (MQ=255)
                                                 cATGCTGGTAGGCACACTGTTCATCCGGCCCAGTTGGCTTGGCTGATCCGAGTGTTCATGTGGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACgg                                                                                                      <  2:6316406/150‑1 (MQ=255)
                                                  aTGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAAGCGCACCACCGAACATGTACAGCCCCCCACGGa                                                                                                     >  1:6555160/1‑150 (MQ=255)
                                                   tgggggTCGGCACCCTGTTCTCCCGGCCCTGTGCGGTTGGCTGACCCGAGTTTTCATGTTGGCCGCCGGCATGCGCCGCATTTGCGCCCGGGCGGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGAt                                                                                                    <  2:2949979/146‑1 (MQ=255)
                                                     gtgttcggCACACTGTTCCGCCTGCCCTTTTCGGTTGGCTTATCCGACTTTTCTTGAGGGCCGCGGGCTGGCGCAGCCTTTGCGCCCGGGCGGTGGTGCGCTTGCTAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                  <  2:1635910/144‑1 (MQ=255)
                                                     cTGGTTGGCACACTGTTCTGCCTGCCCTGTTGGCTGGGCTGATCCGACTGTTCTTGATGGCCGCCGGCCGGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                  <  2:4700073/150‑1 (MQ=255)
                                                       ggTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCGGCGACAAAAAAGCGCACCACCGAACATGTAAAGCCCCCCAACGATGCgg                                                                                                >  2:5968179/1‑150 (MQ=255)
                                                          aGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCtt                                                                                             <  2:7216589/150‑1 (MQ=255)
                                                          aGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCCCCACGGAACATGTACAGCCCCCCACGGGTGCGGCtt                                                                                             >  2:4291762/1‑150 (MQ=255)
                                                          aGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGCCCAAAAGGCGCCCCAACGGAAATGTAAACCCCCCCACCGATGAGGCtt                                                                                             >  2:4049110/1‑150 (MQ=255)
                                                            gCACACTGTACATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACCAAAAGGCGCACCACCGCAATTGTACACCCCGCCACGGATGCGGCTTaa                                                                                           >  2:4005363/1‑150 (MQ=255)
                                                              acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                     >  2:8349968/1‑107 (MQ=255)
                                                               cacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGGCAAAAAGGCGCCCCCCGGCCCGTGTTCAGCCCCCCACGGGGGCGGGGTAAttc                                                                                        >  2:3885041/1‑148 (MQ=255)
                                                                 cgggccctCCGGCCCAGTTCGGTTGGCTGAACCGCCTGTTCAGGTGGGCCGCGGGCATGCGCAGCCATTGCGCCCTCGCCGTGGGGCGCCGGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgc                                                                                      <  2:2552724/143‑1 (MQ=255)
                                                                 cTGTTCAGCCTGCCCTGTTCGCTTGGCTGATCCGAGTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgc                                                                                      <  1:7204466/150‑1 (MQ=255)
                                                                   gTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGTCAAAAAGGCGCACCCCCGGAAGTGGAAAGCCCGCAAACGGTGCGGCGTAATgcgcgc                                                                                    >  2:2252242/1‑150 (MQ=255)
                                                                       aTCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGGCAAAAAGGCGCACCACCGGACGTGTACCGCCCGCAACGGCTGCGGCGTAATGCGCGCCcta                                                                                >  2:8423389/1‑150 (MQ=255)
                                                                        tCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGACAAAAAGGCGCACCCCCGCACATTTAAAGCCCCCAACCGCTGCGGCGTAATGCGCGCCaaaa                                                                               >  2:3933145/1‑147 (MQ=255)
                                                                        tCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTTCGACAAAAAGGCGCACCACGGAACTTGTAAAGCCCGCAAAGGATGCGGCGTAATTCGCGCCatat                                                                               >  2:7803027/1‑150 (MQ=255)
                                                                           tGCCCTGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCg                                                                            <  2:4900234/150‑1 (MQ=255)
                                                                           tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCCCGGGACTTGTACAGCCCGCAACCGATGCGGCTTAATGCGCGCCATATTCg                                                                            >  2:2821352/1‑150 (MQ=255)
                                                                           tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATAttag                                                                            >  2:3504731/1‑148 (MQ=255)
                                                                             cccAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCAAAAAAGCGCCCCACCGAACGTGTTAAGCCCCCAACCGGTGGGGCGGAATGCGCGCCAAATACGcg                                                                          >  1:2805998/1‑149 (MQ=255)
                                                                                cgTTCGCTTGGCTGATCCTACTGTGCATGATGGCCGCCGGCTTGCGCAGCCATTGCGCCCTCGCCGGGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCgg                                                                       <  2:4423822/149‑1 (MQ=255)
                                                                                 gTTCGCTTGGCTGATCCGACTGTTCTTGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCggg                                                                      <  2:8620690/150‑1 (MQ=255)
                                                                                    ggCTTGGATGATCCGACTGTTCTTGATGGCCGCCGGCATGCGCAGCCATTGAGCCCCCGCGGTGGTGCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGa                                                                   <  2:2378405/149‑1 (MQ=255)
                                                                                     ggttGGGTTATCGGACTGTTCTTGTTGGCGGGGGGCTTGCGCTGCCTTTGAGCCCTGGGCGTGGTGCGCCTGCGAGTAAAAGGCGCAGCACGGTACTTGTACAGCCCGCAACGGATGCGGCTTTATGCGCGCCATATTGGCTCGGGTgtc                                                                  <  2:664372/148‑3 (MQ=255)
                                                                                     gCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCCCCACCCAACATGTACAGCCCCCAACCCATGACGCTTAATTCGCGCCAAATTAGCACggggtac                                                                  >  2:5222622/1‑146 (MQ=255)
                                                                                       ttGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCTGCGCCCAAAAGGCGCACCCCCGAACTTGTACAGCCCCCAACCGAGGCGGCTTAATTCGCGCCAAATTCGCTCGGGTGAcac                                                                >  2:6298855/1‑148 (MQ=255)
                                                                                         ggCTGAGCCGACTGTTCATGTTGGCCGCCGGCATGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCAt                                                              <  2:6721002/150‑1 (MQ=255)
                                                                                           cTGTTCCGAGTGTTCTTGCTGCCCGCCGGCTTGCGCAGCAATTGAGCCCTGGCGGTGGTGCGCCTGCCAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGc                                                            <  1:8520239/150‑1 (MQ=255)
                                                                                               tCCGACTGTTGGTGTTGGCGGCCGGCATGCGCAGCATTTGCGCCGTCGCCGTGGTGCGCCTGCCAAAAAACGGGCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTc                                                        <  2:2995027/150‑1 (MQ=255)
                                                                                               cccGACTGTTCCTGATGGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTc                                                        <  1:7937902/149‑1 (MQ=255)
                                                                                                ccGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGCCAAAAAAGCGCACCAACGAACATGTAAAGCCCCCAACCGATGCGGCTTAATTCGCGCCAAAATAGCGCCGGTGAATAATGCTCTca                                                       >  2:3264572/1‑150 (MQ=255)
                                                                                                ccGACTGTTCATGAGGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTca                                                       <  1:8779519/150‑1 (MQ=255)
                                                                                                   aCTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGGCAAAAAAGCGCACCACCGAACGTTTACAGCCCCCCACCGATTCGGCTTAATTCGCGCCCAAATAGCTCGGGTGTCTCATGCGTTcacgt                                                    >  1:7346958/1‑150 (MQ=255)
                                                                                                          atgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGGCAAAAAGGCGCAACACCGCACATGTAAAGCCCCCAAACGATGCGGCTTAATGAGCGCCAAATTCGCTCGGGTGGCGAAAGCTTTcacgtgcacggg                                             >  2:2635763/1‑148 (MQ=255)
                                                                                                           tgtgggCCGCGGGCTGGCGCCGCCCTGGAGCCCTCGCGGGGGGGCGCCGGCCAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCacgtgcacgtgc                                            <  1:8255281/146‑1 (MQ=255)
                                                                                                           tgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGACAAAAAGGCGCACCACGGCACTTTTACAGCCCGCCACGGCTGCGGCGTAATGCGCGCCCAATTCGCTCGGGGGGCGCATGGTTTGacgggcccgggc                                            >  2:5498245/1‑150 (MQ=255)
                                                                                                           tgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGCGGTGCGCCTGCGTCAAAAAGGCGCACCACGGAACTTTTACATCCCGCAACGGATTCGGCTTAATGCGCGCCATATTAGCTAGGGTGAATAATGCTTTAacgtgcacgtgc                                            >  2:3260750/1‑150 (MQ=255)
                                                                                                             atgGCCGCCGGCATGCGCCGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACgtgcacgtgcac                                          <  2:2801877/150‑1 (MQ=255)
                                                                                                                   gcgggattGCGCAGCATTTGCGCCCTGGCGGGGGGGGGCCTGCTAAAAAATGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  1:3270603/143‑1 (MQ=255)
                                                                                                                   gcgggCTTGCGCAGCCATTGCGCCCTCGCCTTGTTGCGCGTGCGACAAAATGGGGCACCTCGGAACTTGTACAGCCCGCAACGGATGGGGCTTAATGCGCTCCATTTTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  1:4276918/147‑1 (MQ=255)
                                                                                                                   gcgggCATGCGCAGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  1:2894267/147‑1 (MQ=255)
                                                                                                                   gccgGCATGCGCAGCCATTGCGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  2:8400561/150‑1 (MQ=255)
                                                                                                                   gccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  2:7147625/150‑1 (MQ=255)
                                                                                                                   gccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  2:7068118/150‑1 (MQ=255)
                                                                                                                   gccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  2:4250374/150‑1 (MQ=255)
                                                                                                                   gccgGCAGGCGCCGCCTTTGAGCCCGCGCCGGGGTGCGCCTGCGACACAACGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                    <  1:8480961/150‑1 (MQ=255)
                                                                                                                    ccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCCCCGGCGACCAAAAGGCGCACCACCGCACGTGTACAGCCCGCAACGGATGCGGCTTTAAGCGCGCCATATTCGCGCGGGTGACTCATGCTTTCACCTGCACGGGCACATCTTTg                                   >  2:363927/1‑150 (MQ=255)
                                                                                                                     cgtcATGCGCAGCAATTGAGCCCTCGCCTTGGTGCGCCTGCTAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                  <  1:5066442/147‑1 (MQ=255)
                                                                                                                     cgGCATGCGCAGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                  <  2:1236402/150‑1 (MQ=255)
                                                                                                                      ggCAGGCGCAGCCATTGCGCCCGCGCGGTGGGGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTc                                 <  2:1039621/150‑1 (MQ=255)
                                                                                                                             gctgCCATTGCGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgc                          <  2:1276899/147‑1 (MQ=255)
                                                                                                                              cAGCCATTGCGCCCTCGCCGGGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcc                         <  2:2468970/150‑1 (MQ=255)
                                                                                                                               aGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcct                        <  1:3575885/150‑1 (MQ=255)
                                                                                                                                 ccATTGAGCCCTCGCCGTGGGGGGCCTGCGCCAAAAAGGCGCACCACCGAACCTGTAAAGCCCCCCACCGCTGCGGCTTAATGGGCGCCAACTTAGCGCGGGGGGCTCAATCGTTAGACTGCACCTGCACAGCGTTTTGTGTTAgcctgc                      >  2:3785510/1‑150 (MQ=255)
                                                                                                                                          ccgcGCCGTGGTGCGCCGGCCAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTg             <  1:7003556/147‑1 (MQ=255)
                                                                                                                                              gCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGt         <  1:1372469/150‑1 (MQ=255)
                                                                                                                                                 gtggtgCGCCTGCGACAAAATGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCa      <  1:3437492/150‑1 (MQ=255)
                                                                                                                                                 gtggtgCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCa      <  2:2905789/150‑1 (MQ=255)
                                                                                                                                                 gtggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGTTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCa      >  1:2529262/1‑150 (MQ=255)
                                                                                                                                                 gtggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCa      >  1:3085771/1‑150 (MQ=255)
                                                                                                                                                  tggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCaa     <  2:8220907/150‑1 (MQ=255)
                                                                                                                                                   ggtgCGCGTGCCACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc    <  1:4788164/150‑1 (MQ=255)
                                                                                                                                                   ggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc    <  2:2046238/150‑1 (MQ=255)
                                                                                                                                                    gtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAcc   >  2:6179531/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CGTCGATCAGCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCT  >  NC_002947/3781216‑3781514

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