Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,571,149 G→A 7.1% intergenic (‑24/+20) PP_2254 ← / ← dctA‑III GntR family transcriptional regulator/C4‑dicarboxylate transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,571,1490GA 7.1% 285.3 / 8.7 98intergenic (‑24/+20)PP_2254/dctA‑IIIGntR family transcriptional regulator/C4‑dicarboxylate transport protein
Reads supporting (aligned to +/- strand):  ref base G (39/52);  new base A (4/3);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

TGCCCAGTTGTTCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCACGGTG  >  NC_002947/2571002‑2571295
                                                                                                                                                   |                                                                                                                                                  
tGCCCAGTTGTTCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTc                                                                                                                                                  <  2:413109/150‑1 (MQ=255)
  cccAGTTGTTCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCAt                                                                                                                                                <  1:69339/150‑1 (MQ=255)
   ccAGTTGTTCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATg                                                                                                                                               >  2:851579/1‑150 (MQ=255)
    cAGTTGTTCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGc                                                                                                                                              >  2:2232266/1‑150 (MQ=255)
     aGTTGTTCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGCCGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGcc                                                                                                                                             >  1:1507879/1‑150 (MQ=255)
      gttgttCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCAGGCGTCATGCCg                                                                                                                                            <  1:1177624/150‑1 (MQ=255)
       ttgttCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCgg                                                                                                                                           <  2:1500017/150‑1 (MQ=255)
       ttgttCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCgg                                                                                                                                           <  1:851579/150‑1 (MQ=255)
        tgttCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGc                                                                                                                                          <  2:2062604/150‑1 (MQ=255)
         gttCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTTCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCAtgcctgcc                                                                                                                                         >  1:70545/1‑150 (MQ=255)
             ggCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTCGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCGTGCCGGCCACGt                                                                                                                                     >  1:2296623/1‑150 (MQ=255)
              gctcgTTCGCGTGGCCCTGGCAGCCTGTGGCCGGGGCGCTATTGGCCGCGGGTCAGGGCATGGCGGAGTTTTTCCACCACTGAGTTGACGTCGGTTGATGAGNTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGtt                                                                                                                                    <  1:899157/146‑1 (MQ=255)
              gCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGtt                                                                                                                                    >  1:2092779/1‑150 (MQ=255)
               ccAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTc                                                                                                                                   >  1:126870/1‑150 (MQ=255)
               ccAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTc                                                                                                                                   <  2:426492/150‑1 (MQ=255)
               ccAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGGTCGCTCCTGGAAGGGCCGGCGTCATGCCGGACACGttt                                                                                                                                   >  1:1611826/1‑149 (MQ=255)
                 aGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGc                                                                                                                                 <  1:1715286/150‑1 (MQ=255)
                          tggcctggcAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCg                                                                                                                                  >  1:288811/1‑140 (MQ=255)
                          tggcctggcAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCg                                                                                                                                  <  2:288811/140‑1 (MQ=255)
                          tggcctggcAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTc                                                                                                                        <  2:841964/150‑1 (MQ=255)
                          tggcctggcAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCATGACGCCGGccc                                                                                                                                        >  2:537893/1‑133 (MQ=255)
                          tggcctggcAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCATGACGCCGGccc                                                                                                                                        <  1:537893/134‑2 (MQ=255)
                                  cagcaTGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGAt                                                                                                                <  2:126870/150‑1 (MQ=255)
                                    gcaTGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATgg                                                                                                              <  2:197654/150‑1 (MQ=255)
                                          cGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCACGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTg                                                                                                        <  1:734537/150‑1 (MQ=255)
                                           gCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTgg                                                                                                       <  1:1747696/150‑1 (MQ=255)
                                              ggTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGtt                                                                                                    <  1:1734439/150‑1 (MQ=255)
                                                   gcCATTGGCCGCGGGCCAGGGCCTGGCGGAGGTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGttgcttg                                                                                               >  1:1923232/1‑150 (MQ=255)
                                                    ccATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGttgcttgc                                                                                              >  1:582249/1‑150 (MQ=255)
                                                    ccATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGAGCTGGGttgcttgc                                                                                              >  2:2337374/1‑150 (MQ=255)
                                                      aTTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgt                                                                                            <  2:1550787/150‑1 (MQ=255)
                                                       ttGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgtg                                                                                           <  1:562121/150‑1 (MQ=255)
                                                           ccGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGATGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgtggcgt                                                                                       >  2:1488686/1‑150 (MQ=255)
                                                           ccGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgtggcgt                                                                                       >  1:138849/1‑150 (MQ=255)
                                                            cgcgGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgtggcgtg                                                                                      >  1:371885/1‑150 (MQ=255)
                                                              cgGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTggg                                                                                    <  1:327587/150‑1 (MQ=255)
                                                               gggCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCt                                                                                                                      <  2:789225/115‑1 (MQ=255)
                                                               gggCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCt                                                                                                                      >  1:789225/1‑115 (MQ=255)
                                                                 gCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCTGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgtg                                                                                           <  1:1301362/140‑1 (MQ=255)
                                                                 gCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCTGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgtg                                                                                           >  2:1301362/1‑140 (MQ=255)
                                                                    aGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTg                                                                                                        <  2:721435/124‑1 (MQ=255)
                                                                    aGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTg                                                                                                        >  1:721435/1‑124 (MQ=255)
                                                                      ggCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGc                                                                            >  2:1514140/1‑150 (MQ=255)
                                                                         cTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgt                                                                                            <  1:1761181/131‑1 (MQ=255)
                                                                         cTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgt                                                                                            >  2:1761181/1‑131 (MQ=255)
                                                                          tGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAgg                                                                        <  2:404655/150‑1 (MQ=255)
                                                                          cggCGGACTTCTGCCCCCGCTTTGTTGATGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAgg                                                                        <  2:138849/149‑1 (MQ=255)
                                                                           ggCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGc                                                                       <  2:1164419/150‑1 (MQ=255)
                                                                                    tttCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGGTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCg                                                              >  1:1428097/1‑150 (MQ=255)
                                                                                     ttCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGa                                                             <  1:2284846/150‑1 (MQ=255)
                                                                                     ttCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGa                                                             <  2:2398213/150‑1 (MQ=255)
                                                                                     ttCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGAGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGa                                                             <  1:2337374/150‑1 (MQ=255)
                                                                                        caccacTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCa                                                          >  1:2272534/1‑150 (MQ=255)
                                                                                         accacTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGc                                                                                                                                 <  2:1599094/78‑1 (MQ=255)
                                                                                         accacTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGc                                                                                                                                 <  2:1599766/78‑1 (MQ=255)
                                                                                         accacTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGc                                                                                                                                 >  1:1599094/1‑78 (MQ=255)
                                                                                         accacTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGc                                                                                                                                 >  1:1599766/1‑78 (MQ=255)
                                                                                         accacTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCATCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCAc                                                         >  2:1959730/1‑150 (MQ=255)
                                                                                           cacTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGt                                                       >  1:1367704/1‑150 (MQ=255)
                                                                                                 gTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTgcg                                                 >  1:1743662/1‑150 (MQ=255)
                                                                                                  ttGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGAt                                                                                                                >  1:293828/1‑86 (MQ=255)
                                                                                                  ttGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGAt                                                                                                                <  2:293828/86‑1 (MQ=255)
                                                                                                  ttGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCATCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAgg                                                                        <  1:46550/126‑1 (MQ=255)
                                                                                                  ttGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCATCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAgg                                                                        >  2:46550/1‑126 (MQ=255)
                                                                                                       gacgGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCt                                                                                                                      >  1:1338295/1‑75 (MQ=255)
                                                                                                       gacgGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCt                                                                                                                      <  2:1338295/75‑1 (MQ=255)
                                                                                                       gacgGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGt                                                                                                                    <  1:1323468/77‑1 (MQ=255)
                                                                                                       gacgGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGt                                                                                                                    >  2:1323468/1‑77 (MQ=255)
                                                                                                           gTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGACCACGTTCGCTCAGAGTTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAg                                       <  1:286050/150‑1 (MQ=255)
                                                                                                             tgatgaGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGt                                                                                                                                     <  1:1119333/54‑1 (MQ=255)
                                                                                                             tgatgaGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGt                                                                                                                                     >  2:1119333/1‑54 (MQ=255)
                                                                                                              gatgaGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCa                                        >  2:2282444/1‑146 (MQ=255)
                                                                                                              gatgaGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCa                                        <  1:2282444/146‑1 (MQ=255)
                                                                                                              gatgaGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTc                                    >  2:1576402/1‑150 (MQ=255)
                                                                                                                  aGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATg                                <  2:2086090/150‑1 (MQ=255)
                                                                                                                   gCTAATCACGGTTCGCTCCTGGAAGGTCAGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGAGGCAAAGCCAGGCTTGCAGTCAAACACGTTCTGCGCCCGTTCCGGGTCAATGt                               >  2:1040913/1‑150 (MQ=255)
                                                                                                                    cTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCTAGCACGTTCTGCGCACGTTCCAGGTCAATGTc                              >  1:1479942/1‑150 (MQ=255)
                                                                                                                    cTGATCACGGTTCGCTCCTCCAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTCTGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTc                              <  1:1931613/150‑1 (MQ=255)
                                                                                                                      gATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCt                            >  1:1697407/1‑150 (MQ=255)
                                                                                                                       aTCACGGTTCGCTCCTGGAAGGGCCGGGGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCtt                           <  2:726425/150‑1 (MQ=255)
                                                                                                                          aCGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTtctc                        >  1:408219/1‑150 (MQ=255)
                                                                                                                          aCGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTtctc                        <  2:582249/150‑1 (MQ=255)
                                                                                                                          aCGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTtctc                        <  2:2092779/150‑1 (MQ=255)
                                                                                                                           cGGTTCGCTCCTGGAAGGGCCGGCATCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTcc                       <  1:1959730/150‑1 (MQ=255)
                                                                                                                                 gCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACgg                 <  2:1697407/150‑1 (MQ=255)
                                                                                                                                    ccTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCa              <  1:790907/150‑1 (MQ=255)
                                                                                                                                      tGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAAt            >  1:1239198/1‑150 (MQ=255)
                                                                                                                                        gAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATgg          <  2:371885/150‑1 (MQ=255)
                                                                                                                                        gAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATgg          <  2:1923232/150‑1 (MQ=255)
                                                                                                                                        gAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCGTTTCCCAACGTGCAATgg          <  2:2393947/150‑1 (MQ=255)
                                                                                                                                           tggCCGGCATGACGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCa       >  1:2042214/2‑150 (MQ=255)
                                                                                                                                           gggCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCa       <  2:2296623/150‑1 (MQ=255)
                                                                                                                                             gCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATg           >  2:1090806/1‑144 (MQ=255)
                                                                                                                                             gCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATg           <  1:1090806/144‑1 (MQ=255)
                                                                                                                                                ggCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCACGGTg  <  1:2232266/150‑1 (MQ=255)
                                                                                                                                                ggCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCACGGTg  <  2:1239198/150‑1 (MQ=255)
                                                                                                                                                 gCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAAt                                 <  2:1220305/118‑1 (MQ=255)
                                                                                                                                                 gCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAAt                                 >  1:1220305/1‑118 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCCCAGTTGTTCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTGACGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCATGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCACGGTG  >  NC_002947/2571002‑2571295

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