Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,596,415 A→T 6.4% L213Q (CTG→CAG)  dctD‑III ← C4‑dicarboxylate transport transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,596,4150AT 6.4% 351.4 / 10.8 109L213Q (CTG→CAG) dctD‑IIIC4‑dicarboxylate transport transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (53/49);  new base T (3/4);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 110 total aligned reads displayed.

GCTCCTGTATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGTTCACCGGCA  >  NC_002947/1596266‑1596562
                                                                                                                                                     |                                                                                                                                                   
gCTCCTGTATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCCTGCCGTCGATTGCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACa                                                                                                                                                     <  2:1270830/150‑1 (MQ=255)
 cTCCTGTATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAg                                                                                                                                                    >  2:899740/1‑150 (MQ=255)
  tCCTGTATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGt                                                                                                                                                   >  1:2085779/1‑150 (MQ=255)
  tCCTGTATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGt                                                                                                                                                   >  2:1804590/1‑150 (MQ=255)
   ccTGTATCACCCGCAACACCTTCGCCTGTTGTGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACCCGTATCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGtt                                                                                                                                                  <  1:762624/150‑1 (MQ=255)
       tATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACTGTTCGGc                                                                                                                                              >  2:1592620/1‑150 (MQ=255)
         tCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACTGTTCGGcca                                                                                                                                            <  2:2262198/150‑2 (MQ=255)
         tCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCt                                                                                                                                            <  2:478595/150‑1 (MQ=255)
            cccGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACTGTTCGGccatga                                                                                                                                         <  1:1592620/150‑5 (MQ=255)
                 aaCAGCTTGGCCTGTTGGGCCGGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGc                                                                                                                                    <  1:1382144/150‑1 (MQ=255)
                 aaCAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGc                                                                                                                                    <  1:2275431/150‑1 (MQ=255)
                 aaCACCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCTGTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGc                                                                                                                                    <  1:779812/150‑1 (MQ=255)
                   cAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGa                                                                                                                                  >  2:229288/1‑150 (MQ=255)
                      cTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTc                                                                                                                               >  1:1349908/1‑150 (MQ=255)
                       ttGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCg                                                                                                                              >  1:1538732/1‑150 (MQ=255)
                        tGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCgg                                                                                                                             <  1:1804590/150‑1 (MQ=255)
                        tGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCgg                                                                                                                             <  2:1538732/150‑1 (MQ=255)
                         ggCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCggg                                                                                                                            <  1:1042029/150‑1 (MQ=255)
                           tattgtGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGa                                                                                                                          >  1:84238/6‑150 (MQ=255)
                           ccTGTTGGGCCAAGCGCATGCTATCGATTTCGTCGAGGAACACGGTGTCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGGTCCCGTGAAGGCGCCGCTCTCATGTCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGa                                                                                                                          <  2:321367/150‑1 (MQ=255)
                            cTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGAt                                                                                                                         <  1:1025686/150‑1 (MQ=255)
                             tGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATg                                                                                                                        <  2:364511/150‑1 (MQ=255)
                                     cAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGa                                                                                                                                  <  1:306347/132‑1 (MQ=255)
                                      aGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAg                                                                                                               >  2:132695/1‑150 (MQ=255)
                                        gggCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGtt                                                                                                             >  1:89141/1‑150 (MQ=255)
                                          gCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTaa                                                                                                           <  2:1224360/150‑1 (MQ=255)
                                           cATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAg                                                                                                          <  2:763158/150‑1 (MQ=255)
                                           cATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAg                                                                                                          <  2:1349908/150‑1 (MQ=255)
                                                  tcGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACa                                                                                                   >  2:1873628/1‑150 (MQ=255)
                                                    gATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGTCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGGTAAGGGCTACaaa                                                                                                 >  2:2130213/1‑150 (MQ=255)
                                                     aTTTCGACGAGGAACCTGTTGCCACCGCTGGCATACTCCAGCTGTCTGATGCGTTGACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAg                                                                                                <  2:573147/150‑1 (MQ=255)
                                                      tttCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACa                                                                                                                                      <  2:512258/111‑1 (MQ=255)
                                                      tttCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACa                                                                                                                                      >  1:512258/1‑111 (MQ=255)
                                                       ttCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAggg                                                                                              <  1:784836/150‑1 (MQ=255)
                                                           tcgAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACa                                                                                                                                      >  2:722540/1‑106 (MQ=255)
                                                           tcgAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACa                                                                                                                                      <  1:722540/106‑1 (MQ=255)
                                                           tcgAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTg                                                                                          >  2:186082/1‑150 (MQ=255)
                                                                     cGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACggg                                                                                <  1:166901/150‑1 (MQ=255)
                                                                           cACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGa                                                                          >  1:1644614/1‑150 (MQ=255)
                                                                           cACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGa                                                                          >  1:1644864/1‑150 (MQ=255)
                                                                            aCCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACTGTTCGGcca                                                                                                                                            <  2:2002244/83‑2 (MQ=255)
                                                                            aCCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACTGTTCGGcca                                                                                                                                            >  1:2002244/1‑82 (MQ=255)
                                                                            aCCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAg                                                                                                          >  2:362951/1‑117 (MQ=255)
                                                                            aCCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAg                                                                                                          <  1:362951/117‑1 (MQ=255)
                                                                            aCCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGaa                                                                         >  1:861600/1‑150 (MQ=255)
                                                                             ccGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAg                                                                        <  2:1644614/150‑1 (MQ=255)
                                                                             ccGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAg                                                                        <  2:1644864/150‑1 (MQ=255)
                                                                              cGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGt                                                                       >  2:1149007/1‑150 (MQ=255)
                                                                               gTTGGCATACTCCAGCTTGCCGATGCTTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTc                                                                      >  2:2395979/1‑150 (MQ=255)
                                                                                ttGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCg                                                                     <  1:186082/150‑1 (MQ=255)
                                                                                    cATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTgca                                                                 <  1:132695/150‑1 (MQ=255)
                                                                                        cTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTgcaggca                                                             >  1:2134369/1‑150 (MQ=255)
                                                                                         tCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTgcaggcag                                                            >  2:1888856/1‑150 (MQ=255)
                                                                                           cAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCg                                                          >  1:1965164/1‑150 (MQ=255)
                                                                                           cAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAAGGAGCg                                                          >  2:255487/1‑150 (MQ=255)
                                                                                            aGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACGCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCgg                                                         <  1:899740/150‑1 (MQ=255)
                                                                                                     aTGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGc                                                >  2:847304/1‑150 (MQ=255)
                                                                                                     aTGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGc                                                <  2:2170118/150‑1 (MQ=255)
                                                                                                      tGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCt                                               >  1:963419/1‑150 (MQ=255)
                                                                                                        cGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTcc                                             <  1:1453302/150‑1 (MQ=255)
                                                                                                         gTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTggcgg                                                                                      >  1:369655/1‑108 (MQ=255)
                                                                                                         gTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTggcgg                                                                                      <  2:369655/108‑1 (MQ=255)
                                                                                                          tttACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCtt                                           <  2:600103/150‑1 (MQ=255)
                                                                                                           ttACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTg                                          <  2:89141/150‑1 (MQ=255)
                                                                                                             aCCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGcc                                        >  1:132301/1‑150 (MQ=255)
                                                                                                               gcaatgcgcCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCAc                                      >  1:84648/6‑150 (MQ=255)
                                                                                                               ccTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGaa                                                                         >  2:1485939/1‑115 (MQ=255)
                                                                                                               ccTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGaa                                                                         <  1:1485939/115‑1 (MQ=255)
                                                                                                               ccTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCAc                                      >  1:1727942/1‑150 (MQ=255)
                                                                                                                cTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACt                                     >  1:2380984/1‑150 (MQ=255)
                                                                                                                  gggCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTaa                                                                                                           >  1:405441/1‑78 (MQ=255)
                                                                                                                  gggCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTaa                                                                                                           <  2:405441/78‑1 (MQ=255)
                                                                                                                      gcCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGCAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGt                               <  2:84238/150‑1 (MQ=255)
                                                                                                                      gcCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGt                               <  1:2143495/150‑1 (MQ=255)
                                                                                                                       cccGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGtt                              <  1:1149007/150‑1 (MQ=255)
                                                                                                                        ccGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGttt                             >  2:903233/1‑150 (MQ=255)
                                                                                                                           tGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGTAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCAc                          >  2:1619015/1‑150 (MQ=255)
                                                                                                                           tGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCAc                          >  2:523679/1‑150 (MQ=255)
                                                                                                                            gAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCAcc                         <  1:1619015/150‑1 (MQ=255)
                                                                                                                               ggCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACcgcg                      <  1:229288/150‑1 (MQ=255)
                                                                                                                                gcgcCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCtt                                           <  1:2221617/128‑1 (MQ=255)
                                                                                                                                gcgcCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCtt                                           >  2:2221617/1‑128 (MQ=255)
                                                                                                                                  gccgcTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGAt                   >  1:732148/1‑150 (MQ=255)
                                                                                                                                   ccgcTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATc                  >  2:2357365/1‑150 (MQ=255)
                                                                                                                                    cgcTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCa                 >  1:1362654/1‑150 (MQ=255)
                                                                                                                                       tctcATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCGGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAaca              >  2:1426369/1‑150 (MQ=255)
                                                                                                                                          cATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGaa                                                                         <  1:699295/88‑1 (MQ=255)
                                                                                                                                          cATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGaa                                                                         >  2:699295/1‑88 (MQ=255)
                                                                                                                                           gaggCCGAACTGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGaa                                                                         <  2:1320718/85‑1 (MQ=255)
                                                                                                                                           gaggCCGAACTGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGaa                                                                         >  1:1320718/3‑87 (MQ=255)
                                                                                                                                           aTGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGtt          >  1:326815/1‑150 (MQ=255)
                                                                                                                                            tGGCCGAACAGTTCGACCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAg                                                                                                          <  2:758295/53‑1 (MQ=255)
                                                                                                                                               ccGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGTTCAcc      <  2:2380984/150‑1 (MQ=255)
                                                                                                                                                cGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAggg                                                                                              <  2:653052/61‑1 (MQ=255)
                                                                                                                                                cGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAggg                                                                                              >  1:653052/1‑61 (MQ=255)
                                                                                                                                                cGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGc                                         <  1:920528/114‑1 (MQ=255)
                                                                                                                                                cGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGTTCACCg     >  2:172742/1‑150 (MQ=255)
                                                                                                                                                 gAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGTTCACCgg    <  2:1965164/150‑1 (MQ=255)
                                                                                                                                                   aCAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGTTCACCGGCa  <  1:2279710/150‑1 (MQ=255)
                                                                                                                                                   aCAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGTTCACCGGCa  <  2:1573735/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTCCTGTATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGTTCACCGGCA  >  NC_002947/1596266‑1596562

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