Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,677,303 T→G 9.6% intergenic (+34/+70) PP_4137 → / ← chrR outer membrane siderophore receptor/chromate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,677,3030TG 9.6% 281.3 / 16.5 102intergenic (+34/+70)PP_4137/chrRouter membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand):  ref base T (21/68);  new base G (5/5);  total (29/73)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 102 total aligned reads displayed.

TGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTC  >  NC_002947/4677161‑4677442
                                                                                                                                                 |                                                                                                                                            
tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGTCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGCATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAAttt                                                                                                                                      <  2:2827482/150‑1 (MQ=255)
tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAAttt                                                                                                                                      <  1:265930/150‑1 (MQ=255)
tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAAttt                                                                                                                                      <  2:3213715/150‑1 (MQ=255)
tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAAttt                                                                                                                                      <  1:221941/150‑1 (MQ=255)
tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAAttt                                                                                                                                      <  1:1116646/150‑1 (MQ=255)
  cAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGc                                                                                                                                    <  1:546905/150‑1 (MQ=255)
      gACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGAc                                                                                                                                <  1:4779669/150‑1 (MQ=255)
       aCCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTTACCCCGAATTCGTCAGTGAATTTTCCGAcc                                                                                                                               >  1:437016/1‑149 (MQ=255)
       aCCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCCCGCATTCCGCATTTAAATTTCAGACg                                                                                                                               >  1:436750/1‑150 (MQ=255)
        ccGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGc                                                                                                                              <  2:1309216/150‑1 (MQ=255)
          gTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCAt                                                                                                                            <  2:3665596/150‑1 (MQ=255)
              cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTcgc                                                                                                                        <  1:4707577/150‑1 (MQ=255)
              cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTcgc                                                                                                                        <  1:3763054/150‑1 (MQ=255)
              cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTcgc                                                                                                                        <  1:3763219/150‑1 (MQ=255)
              cGCCGGCCGACCCCGACGCCGGCCAGGACCCGGCCCGTTCGACCCCGGGCGAAGGATGTGGCGTGTAGGCCGGGGTTTCCTACCGCTTCCGCTGACGGGTGGGCGCGGGTGTACCCGCGCATGCGTCATTGAATGTGCCGACGCATTcgc                                                                                                                        <  2:1556516/150‑1 (MQ=255)
               gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTcgcg                                                                                                                       <  2:1367917/150‑1 (MQ=37)
               gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGaa                                                                                                                                         >  2:4680640/1‑132 (MQ=255)
               gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGaa                                                                                                                                         <  1:4680640/132‑1 (MQ=255)
                ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTCGCgg                                                                                                                      <  2:3402003/150‑1 (MQ=37)
                ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTCGCgg                                                                                                                      <  2:2634065/150‑1 (MQ=37)
                ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTCGCgg                                                                                                                      <  2:443606/150‑1 (MQ=37)
                ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCgg                                                                                                                      <  2:3446300/150‑1 (MQ=255)
                ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCgg                                                                                                                      <  2:5300967/150‑1 (MQ=255)
                ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCgg                                                                                                                      <  2:4565245/150‑1 (MQ=255)
                ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCgg                                                                                                                      <  2:2052829/150‑1 (MQ=255)
                ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGAGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTCGCgg                                                                                                                      <  2:4215963/150‑1 (MQ=37)
                  ggCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGt                                                                                                                    >  2:5341912/1‑150 (MQ=255)
                     cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTAACCCCCAAAGCGTCAATGAATTTTACGACGCATTAGCCGGGaaa                                                                                                                 >  1:4484483/1‑150 (MQ=255)
                      tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAAc                                                                                                                <  1:1996649/150‑1 (MQ=255)
                       aCAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTTCCCGCGCATGCGTCATTGAATTTGCCGGCGCCTTAGCGGGTaaaac                                                                                                               >  2:1663163/1‑148 (MQ=255)
                        cAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAAccc                                                                                                              <  1:4063022/150‑1 (MQ=255)
                         aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGAc                                                                                                                                >  2:1989606/1‑131 (MQ=255)
                         aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGAc                                                                                                                                <  1:4728946/131‑1 (MQ=255)
                         aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGAc                                                                                                                                <  1:1989606/131‑1 (MQ=255)
                         aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGAc                                                                                                                                >  2:4728946/1‑131 (MQ=255)
                           cgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCt                                                                                                           <  2:4166837/150‑1 (MQ=255)
                            gccgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACCGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCCACGCATTCGCGGGTAAACCCGCTc                                                                                                          <  2:4597795/148‑1 (MQ=255)
                                      aGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAggg                                                                                                <  2:1748011/150‑1 (MQ=255)
                                           aaTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGAATTCGCGGGTAAACCCGCTCCTACACAGGGGCGtt                                                                                           <  2:1642017/150‑1 (MQ=255)
                                           aaTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTGCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGtt                                                                                           <  2:3268618/150‑1 (MQ=255)
                                            aTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTg                                                                                          <  2:1770174/150‑1 (MQ=255)
                                             tGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTCCCGCTGACGGGTGGGCGCGGTTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGc                                                                                         <  1:4932882/150‑1 (MQ=255)
                                              gCCCGTTCGACACCGGGCGAAGTTTTGGGGGTGTATGCGGGGGTTCCCTACAGCTTCCGGTGACGGGTGGGCGCGGTTTTACCCGCGAATGCGTATTTGAATTTGCCGACGAATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGca                                                                                        <  2:1430412/150‑1 (MQ=255)
                                              gCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGca                                                                                        <  2:4329158/150‑1 (MQ=255)
                                                ccGTTCGACACCGGGCGATGGTTTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGcaca                                                                                      <  1:4953213/150‑1 (MQ=255)
                                                 cGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGGTCCCCGCGCATGCGGCCTGGCATTTGCCCCCGCCGGCCCGGGGAAACCCGCGCCGACACCGGGGCGGTGCAcag                                                                                     >  2:3573696/1‑150 (MQ=255)
                                                  gTTCGACACCGGGCGATGGATTGGGCGTGTATGCCGGGGTTTCCTCCAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCAcagc                                                                                    <  2:3735245/150‑1 (MQ=255)
                                                   ttCGACACCGGGCTAAGTATTTGGCGTGTATGCCGGGGTTTCCTACAGCTCCCGCTGACGGGTGGGCGCGGGTTTACCCGCTAATGCTTCATTTAATTGGCCCACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCAcagca                                                                                   <  2:622517/150‑1 (MQ=38)
                                                        cacCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTTACCCCCGAATCGTCATTGAATTTTCCCACGCAGTCGCGGGGAAACCCCCGCCTACACAGGGGGGGTGGACAGCAGGGtt                                                                              >  1:3597797/1‑150 (MQ=255)
                                                         acCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTAACCCCGAATTCGTCATTGTATTTTGCGACCCATTCGCCGGTAAAACCCCTCCAAAAAAGGGGCGGTTCCCAGCAGGCTTa                                                                             >  2:3405312/1‑150 (MQ=255)
                                                         acCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTACCCGCGACAGCGGCCATGCAGTTTGCGCCGCAGTCGCGGGGGAACCCGCGCCGACACCGGGGCGGTGCCCAGCCggggtt                                                                             >  2:331027/1‑146 (MQ=255)
                                                          ccGGGCGATGGATTTGGCGTTTATGCCGGGTTTTCCTACAGCTTCCGCTGACGGGTGGGCGGGGTTTTACCCGCGAATGCTTCATTGAATTTGCCGACGCTTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAg                                                                            <  1:801518/150‑1 (MQ=255)
                                                            gggCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTAACCCCCCAATCGTCATTGAAATTTCCGACCCAATAGCGGGTAAAACCCCAACAACACACGGGCGGTGCACAGCCGGGATAggg                                                                          >  2:4205153/1‑150 (MQ=255)
                                                            gggCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTAACCCCCAAAGCGGGAATGTATTTTGCCACCCAATAGCGGGTAAAACCCCACCAAAACACGGGCGGTGCACAGCACGATTAggg                                                                          >  2:4205146/1‑150 (MQ=255)
                                                                gacggATTGGGCGTGTCTGCCGGGGTGTCCTCCCGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTc                                                                      <  2:583681/147‑1 (MQ=38)
                                                                 ttgtttttGGCGTTTTTGCGGGGTTTTCCTACTGCTTCCGGTGCGGGGGGGGGGGGGTTTTACCCGCTAATGCTTCTTTTATTTTGCCGACGAATGGGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCa                                                                     <  1:1119245/145‑1 (MQ=255)
                                                                      tttGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTTACCGCGAATTCGTCATTGAATTTTCCGACGCATTCGCGGGTAAACCCGCGCCTAAAAAGGGGCGTTGCACAGCAGGGTTAGGGGTTGAGACCg                                                                >  2:329190/1‑150 (MQ=255)
                                                                        tGGCGCGTAGGCCGGGGTTTCCAACCGCGCCCGCTGACGGGTGGGCGCGGGTGGCCCCGCGAATGCGACATGGAATTTGCCGACGCATGCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                              <  2:5097005/150‑1 (MQ=255)
                                                                         ggCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTTCCCCCGAAAGCGGCCTTGGAATTTGCCACGCAATAGCGGGTAAAACCCCGCCGACACAAGGGCGGTGCACCCCCGGATTAGGGGTTGAGAcccccc                                                             >  2:1187038/1‑150 (MQ=255)
                                                                          gCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTAACCCCGAAAGCGTCATTGAATTTTGCGACCCATTAGCGGGTAAAACCGCACCTAAAAAAGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcccccca                                                            >  2:3502255/1‑149 (MQ=255)
                                                                          gCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTAACCCCGAATTGGGGATTTAATTTTGCGACGCAATGGCGGGTTAACCCCCACCAAAACACGGGCGGTTCACAGCAGGATTATGGGTTAAGAACccccc                                                            >  2:3502333/1‑149 (MQ=255)
                                                                            gtgtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTACCCGCGAATGGGGCCTTTGGTTTGGCGGCGCATTCGCGGGGTAACCCCCTCCGACACAGGGGCGGTGCACCGCCGGGTTAGGGGTTGAGACCCCCCGGt                                                          >  2:4798119/1‑150 (MQ=255)
                                                                              gtATGCGGGGGTTTCCTACCGCTTCCGCTGTGGGGTGGGCGCGGGTTTACCCGCGATTGCTTATTTTATTTTGCCGACGCTTTGGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                        <  1:4205153/150‑1 (MQ=255)
                                                                              gtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                        <  2:4972664/150‑1 (MQ=255)
                                                                              gtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                        <  2:435105/150‑1 (MQ=255)
                                                                               tttggcGGGGTTTGCCTACCGCTGCGGCTGCGGGGGGGGCGCGGGTTTCCCCGGTATTGCTTATTTGATTTGGCCGACGATTTCGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                        <  1:4205146/145‑1 (MQ=255)
                                                                                ttGCGGGGTTTTCCTCCTGCTTCCGCTTCGGGGTGGGGGGGGTTTTCCCGGCTAGTGCTTCATTTATTTTGCCGCCGCTTTGGGGGTAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                      <  1:4674234/149‑1 (MQ=255)
                                                                                tggcgggggTTTCCTACCGCTTCCGCTGAGGGGTGGGCGGGGTTTTCCCCGCTAATGCGTATTTTAATTGGCCGACGCATTCGGGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                      <  1:3237869/145‑1 (MQ=255)
                                                                                aTGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                      <  1:4955390/150‑1 (MQ=255)
                                                                                 tGCCGGGGTTTCCTACAGCTGCCGCTGACGGGTGGGCGCGGGTTTACCGGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAAc                                                     <  2:4656945/150‑1 (MQ=255)
                                                                                       ggTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTTAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                               <  2:4707938/150‑1 (MQ=255)
                                                                                       ggTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                               <  2:4354594/150‑1 (MQ=255)
                                                                                          ttCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCAt                                            <  2:1590183/150‑1 (MQ=255)
                                                                                             ccTACTGCTTCCGCTGACGGGTGGGCGGGGGTTTACCCGCTATTGCGTAATTTAATTTGCCGACGCTTTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGa                                          <  1:3502255/150‑1 (MQ=255)
                                                                                             ccTACAGCTCCCGCTGACGGGTGGGCGCGGGTTTACCCGCTAATGCGTCATTTATTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGa                                          <  1:3502333/150‑1 (MQ=255)
                                                                                              atACTGCTTCGGGTGCGGGGTGGGGGGGGTTTTACCCGCTATTGCTTCTTTTATTTTGCCGCCGATTTCGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                         <  1:189902/149‑1 (MQ=255)
                                                                                                    cTTCCGCTGAGGGGTGGGCGCGGGTTTACCCGCGATTGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                   <  2:2017390/150‑1 (MQ=255)
                                                                                                      tCCGCTGCGGGGGGGGCGCGGGTTGACCCGCGAAGGCGTCATTGAATTTGCCGACGCTTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCt                                 <  1:1187038/150‑1 (MQ=255)
                                                                                                        cggggccGGGTGGGCGCGGGTTTACCCGCGATTGCGTCTTTGAATTTGCCGCCGCTTTGGGGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGt                               <  2:2274546/144‑1 (MQ=255)
                                                                                                           ggACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGAATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGa                            <  2:1587192/149‑1 (MQ=255)
                                                                                                                    ggCGCGGGTTTACCCGCGAATGCGTAATTTATTTTGCCGACGCTTTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCt                   <  1:3281242/150‑1 (MQ=255)
                                                                                                                     gcgcGGGTTTACCCGCTAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                  <  2:436750/150‑1 (MQ=255)
                                                                                                                     gcgcGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                  <  2:436768/150‑1 (MQ=255)
                                                                                                                     gcgcGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                  <  2:437016/150‑1 (MQ=255)
                                                                                                                       ggggggttttACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                  <  2:2341525/143‑1 (MQ=255)
                                                                                                                       gcgGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCa                >  1:2443995/1‑150 (MQ=255)
                                                                                                                            gTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                          <  2:1633197/104‑1 (MQ=255)
                                                                                                                            gTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                          >  1:1633197/1‑104 (MQ=255)
                                                                                                                            gTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa            <  2:1689060/150‑1 (MQ=255)
                                                                                                                                  cccGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc       <  2:2443995/150‑1 (MQ=255)
                                                                                                                                   ccGCGAATGCGTCGGCAAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg      >  1:2396191/1‑150 (MQ=37)
                                                                                                                                   ccGCGAATGCGTCGGCAAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg      >  2:237113/1‑150 (MQ=37)
                                                                                                                                   ccGCGAATGCGTCGGCAAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATCACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg      >  2:5112792/1‑150 (MQ=37)
                                                                                                                                   ccGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg      >  1:5369882/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg      >  2:2874174/1‑150 (MQ=255)
                                                                                                                                    cgcgAATGCGTCTTTGAATTGGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc     <  1:4296739/150‑1 (MQ=255)
                                                                                                                                    cgcgAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc     <  2:821310/150‑1 (MQ=255)
                                                                                                                                    cgcgAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc     >  1:5266800/1‑150 (MQ=255)
                                                                                                                                     gcgACTGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACgggcatt    <  2:2145972/150‑6 (MQ=255)
                                                                                                                                       gAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTc  >  2:1947230/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
TGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTC  >  NC_002947/4677161‑4677442

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: