Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,644,991 C→T 6.6% intergenic (‑43/+23) PP_3211 ← / ← PP_3212 ABC transporter ATP‑binding protein/Rieske 2Fe‑2S family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,644,9910CT 6.6% 268.0 / 7.0 91intergenic (‑43/+23)PP_3211/PP_3212ABC transporter ATP‑binding protein/Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base C (41/44);  new base T (3/3);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

CTGTAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAG  >  NC_002947/3644842‑3645140
                                                                                                                                                     |                                                                                                                                                     
cTGTAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGgc                                                                                                                                                       >  1:1272394/1‑150 (MQ=255)
 tGTAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGgcg                                                                                                                                                      >  1:1724913/1‑150 (MQ=255)
      gTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCg                                                                                                                                                 >  1:2084540/1‑150 (MQ=255)
           ttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATACCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGTGCGCCGGGGca                                                                                                                                            <  2:2182234/150‑2 (MQ=255)
             ccACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTACGCCCCGGTGCGCCGGGGcata                                                                                                                                          <  1:1413082/150‑4 (MQ=255)
              cACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTg                                                                                                                                         <  1:1050969/150‑1 (MQ=255)
               aCCTGATTGGCGAACAAATTGTTTTCCACTCCGATGGTGTTCGGCATGGCCTGCTGCTGTACTTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGTTCATGAGCGGTGGTGTATGCCCTGGCGCACCGTTGTGTTgg                                                                                                                                        <  2:481262/150‑1 (MQ=255)
                 cTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGc                                                                                                                                      >  2:508816/1‑150 (MQ=255)
                   gATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCaa                                                                                                                                    <  1:105624/150‑1 (MQ=255)
                   gATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCATGCGGCTGGACGCCGACCTGGGATAAGTTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCaa                                                                                                                                    <  1:863312/150‑1 (MQ=255)
                    aTTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAg                                                                                                                                   <  2:546327/150‑1 (MQ=255)
                    aTTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGGGTATGCCCCGGCGCACCGGGGCGTTGGGCGAg                                                                                                                                   >  1:215505/1‑150 (MQ=255)
                      tGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAgtg                                                                                                                                 <  2:1113988/150‑1 (MQ=255)
                       ggCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAgtgt                                                                                                                                <  2:1998039/150‑1 (MQ=255)
                         cGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCa                                                                                                                              <  1:973479/150‑1 (MQ=255)
                           aaCAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGt                                                                                                                            <  1:1119498/150‑1 (MQ=255)
                            aCAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTg                                                                                                                           <  1:832316/150‑1 (MQ=255)
                             cAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGTCAAGTGTCAGTgg                                                                                                                          >  2:79277/1‑150 (MQ=255)
                              aaaGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGc                                                                                                                         >  1:462853/1‑150 (MQ=255)
                                 ggggTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                      <  2:2292956/150‑1 (MQ=255)
                                        ccACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAacac                                                                                                               >  1:8302/1‑150 (MQ=255)
                                        ccACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAacac                                                                                                               >  2:2207998/1‑150 (MQ=255)
                                            acCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCg                                                                                                           >  1:360937/1‑150 (MQ=255)
                                            acCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCg                                                                                                           >  2:2071895/1‑150 (MQ=255)
                                                 tGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCtctc                                                                                                      <  1:508816/150‑1 (MQ=255)
                                                                     ggCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAc                                                                                  >  2:1941755/1‑150 (MQ=255)
                                                                      gCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAcc                                                                                 <  1:782078/150‑1 (MQ=255)
                                                                       cTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGc                                                                                      <  1:331688/144‑1 (MQ=255)
                                                                       cTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGc                                                                                      >  2:331688/1‑144 (MQ=255)
                                                                           aCGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCagta                                                                            <  1:1717994/150‑1 (MQ=255)
                                                                             gTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCagtagt                                                                          <  2:1724913/150‑1 (MQ=255)
                                                                                ggCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGc                                                                       <  1:899274/150‑1 (MQ=255)
                                                                                ggCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGc                                                                       <  2:8302/150‑1 (MQ=255)
                                                                                  ccTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAccc                                                                                >  1:2294147/1‑139 (MQ=255)
                                                                                  ccTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAccc                                                                                <  2:2294147/139‑1 (MQ=255)
                                                                                  ccTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCg                                                                     >  2:522082/1‑150 (MQ=255)
                                                                                  ccTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCg                                                                     >  1:1278607/1‑150 (MQ=255)
                                                                                   cTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAcc                                                                                 <  1:431254/137‑1 (MQ=255)
                                                                                   cTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAcc                                                                                 >  2:431254/1‑137 (MQ=255)
                                                                                   cTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACAATTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGa                                                                    <  1:2289626/150‑1 (MQ=255)
                                                                                        ttAGGGTCTGCATTAACATAGCGCTCTCCCCGTGGGTCATGGGGGGTGGTGGGTGCCCCGGCGCACCGGTGCGTGGGGCAATTGTCCGTGTCTGAGGGACACTTCGCTCGCGCCTTGGGCTTCAAGCAGAGCCGGTAGTTGCCGActcgt                                                               <  2:215505/149‑5 (MQ=255)
                                                                                        aTAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGAt                                                               >  1:358565/1‑150 (MQ=255)
                                                                                             gTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATgg                                                                                               >  2:1893300/1‑113 (MQ=255)
                                                                                             gTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATgg                                                                                               <  1:1893300/113‑1 (MQ=255)
                                                                                              tCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGGCCCAGTAGTTGCCGACGGATaaccaa                                                         >  1:1628578/1‑150 (MQ=255)
                                                                                                  cATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCg                                                     >  2:1825322/1‑150 (MQ=255)
                                                                                                   aTCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCAt                                                                                                 <  1:1708931/105‑1 (MQ=255)
                                                                                                   aTCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCAt                                                                                                 >  2:1708931/1‑105 (MQ=255)
                                                                                                   aTCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCgg                                                    >  2:1069840/1‑150 (MQ=255)
                                                                                                    tCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGa                                                                                   <  1:858139/118‑1 (MQ=255)
                                                                                                    tCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGa                                                                                   >  2:858139/1‑118 (MQ=255)
                                                                                                     cAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGc                                                  <  1:1941755/150‑1 (MQ=255)
                                                                                                       aCATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCtt                                                <  1:1174560/150‑1 (MQ=255)
                                                                                                       aCATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCtt                                                <  1:209524/150‑1 (MQ=255)
                                                                                                        cATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTg                                               >  2:912004/1‑150 (MQ=255)
                                                                                                          tAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTc                                             <  2:358565/150‑1 (MQ=255)
                                                                                                            gcgcTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCCCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCg                                           >  2:752235/1‑150 (MQ=255)
                                                                                                              gcTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGAt                                         <  2:462853/150‑1 (MQ=255)
                                                                                                               ctctCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATg                                        >  1:1029670/1‑150 (MQ=255)
                                                                                                                tctcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATgg                                       <  2:1278607/150‑1 (MQ=255)
                                                                                                                 ctcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGc                                      <  1:1069840/150‑1 (MQ=255)
                                                                                                                  tcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCa                                     >  2:702082/1‑150 (MQ=255)
                                                                                                                   ccccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATg                                        <  2:1745328/146‑1 (MQ=255)
                                                                                                                   ccccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATg                                        >  1:1745328/1‑146 (MQ=255)
                                                                                                                    cccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtg                                   >  1:198095/1‑150 (MQ=255)
                                                                                                                      cGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtgtg                                 >  2:2195586/1‑150 (MQ=255)
                                                                                                                       gTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGa                                                                                   <  2:89994/99‑1 (MQ=255)
                                                                                                                       gTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGa                                                                                   >  1:89994/1‑99 (MQ=255)
                                                                                                                            tCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCACCCCGGGCTTGTCCGATGGCATGTGCAGTTc                           >  1:475163/1‑150 (MQ=255)
                                                                                                                            tCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTc                           >  2:1888068/1‑150 (MQ=255)
                                                                                                                              aTGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACACAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCa                         >  1:803410/1‑150 (MQ=255)
                                                                                                                              aTGAGCGGTGGGGTAGGCCCCGGCGCACCTGGGCGTTGGCTAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAGGCAGACTCAGTAGTTTCAGGCGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCa                         >  2:119577/1‑150 (MQ=255)
                                                                                                                                      tggtgTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTgcgg                 >  2:2034385/1‑150 (MQ=255)
                                                                                                                                        gtgtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTc                                                                                          <  2:1645037/75‑1 (MQ=255)
                                                                                                                                        gtgtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTc                                                                                          >  1:1645037/1‑75 (MQ=255)
                                                                                                                                        gtgtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTgcggcg               <  2:509262/150‑1 (MQ=255)
                                                                                                                                        gtgtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTgcggcg               >  1:509262/1‑150 (MQ=255)
                                                                                                                                         tgtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg              <  2:475163/150‑1 (MQ=255)
                                                                                                                                         tgtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg              <  1:2207998/150‑1 (MQ=255)
                                                                                                                                         ccaacgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg              >  2:542142/6‑150 (MQ=255)
                                                                                                                                           aacgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGgt            >  1:1782565/4‑150 (MQ=255)
                                                                                                                                             tGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCAcc                       <  2:1362561/137‑1 (MQ=255)
                                                                                                                                             tGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCAcc                       >  1:1362561/1‑137 (MQ=255)
                                                                                                                                              gCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCgg                                                    <  1:2251560/107‑1 (MQ=255)
                                                                                                                                              gCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCgg                                                    >  2:2251560/1‑107 (MQ=255)
                                                                                                                                              gCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTc         <  1:912004/150‑1 (MQ=255)
                                                                                                                                              gCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTc         <  2:1029670/150‑1 (MQ=255)
                                                                                                                                               ccccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGCGGAACACTTCGCTCTCGCAATGGGATTCAAGCAGACGCAGTATTTGCCGACGGATCACCAACCCGGGCTTGACCGCTGGCATGTGCTGTTCCACCCTGCGGCGGGTGTCg        >  2:271687/1‑150 (MQ=255)
                                                                                                                                                cccGTCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCGCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGCCGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGa       <  1:271687/150‑1 (MQ=255)
                                                                                                                                                 ccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGac      >  1:1320995/1‑150 (MQ=255)
                                                                                                                                                 ccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGa       >  2:1050378/1‑149 (MQ=255)
                                                                                                                                                 ccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGa       <  1:1050378/149‑1 (MQ=255)
                                                                                                                                                  cGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAGCCCGGGCTTGTCCGCTGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGaca     <  1:1888068/150‑1 (MQ=255)
                                                                                                                                                    gcgcACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGAcaca   <  2:2084540/150‑1 (MQ=255)
                                                                                                                                                     cgcACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAg  <  1:2071895/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGTAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAG  >  NC_002947/3644842‑3645140

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