Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,347,206 C→T 6.5% intergenic (‑35/+20) PP_1166 ← / ← PP_1167 DMT superfamily permease/TRAP dicarboxylate transporter subunit DctM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,347,2060CT 6.5% 759.8 / 26.9 244intergenic (‑35/+20)PP_1166/PP_1167DMT superfamily permease/TRAP dicarboxylate transporter subunit DctM
Reads supporting (aligned to +/- strand):  ref base C (101/127);  new base T (9/7);  total (110/134)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 257 total aligned reads displayed.

AAGGCGATGCGCGCGGCGATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGATACTCACCTTGCCAATG  >  NC_002947/1347059‑1347355
                                                                                                                                                   |                                                                                                                                                     
aaGGCGATGCGCGCGGCGATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCaa                                                                                                                                                     <  2:5187002/150‑1 (MQ=255)
   gCGATGCGCGCGGCGATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCg                                                                                                                                                  >  2:477111/1‑150 (MQ=255)
     gATGCGCGCGGCGATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGc                                                                                                                                                <  2:2628264/150‑1 (MQ=255)
     gATGCGCGCGGCGATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGc                                                                                                                                                >  1:992444/1‑150 (MQ=255)
       tGCGCGCGGCGATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCCGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGccc                                                                                                                                              >  2:112116/1‑150 (MQ=255)
           cgcgGCGATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTATAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCgg                                                                                                                                          <  2:2487746/150‑1 (MQ=255)
           cgcgGCGATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCGGGGGCCGCTTGGCGGCCCCTgg                                                                                                                                          <  1:74990/150‑1 (MQ=255)
                 gttgtgGTTGCCAGCTAATGAACAGGCCTGGTGGGCGAGCAGGGGGATGGCTTGTTTGCGTGGGAAAAGCGGTACAGCCGAGGTGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATaa                                                                                                                                    <  1:5220653/148‑1 (MQ=255)
                 gATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATaa                                                                                                                                    <  1:5462929/150‑1 (MQ=255)
                 gATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATaa                                                                                                                                    >  2:4180661/1‑150 (MQ=255)
                    gtgGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTa                                                                                                                                 <  1:3200680/150‑1 (MQ=255)
                       gTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAg                                                                                                                              >  1:4060363/1‑150 (MQ=255)
                        ttGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGc                                                                                                                             >  1:2365277/1‑150 (MQ=255)
                          gCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGcac                                                                                                                           >  1:1512582/1‑150 (MQ=255)
                             aGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGcaccac                                                                                                                        >  2:1910609/1‑150 (MQ=255)
                             aGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGcaccac                                                                                                                        <  1:5399071/150‑1 (MQ=255)
                                gAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCt                                                                                                                     >  1:4864852/1‑150 (MQ=255)
                                   tGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAgg                                                                                                                  <  1:4639334/150‑1 (MQ=255)
                                    gAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGc                                                                                                                 <  2:5117905/150‑1 (MQ=255)
                                    gAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGc                                                                                                                 <  2:655935/150‑1 (MQ=255)
                                    gAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGc                                                                                                                 >  1:1766557/1‑150 (MQ=255)
                                     aaCAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCa                                                                                                                <  1:4446863/150‑1 (MQ=255)
                                       cAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGc                                                                                                              <  2:3605901/150‑1 (MQ=255)
                                          gCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCcac                                                                                                           >  2:4076096/1‑150 (MQ=255)
                                           ccAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCcaca                                                                                                          >  2:3427641/1‑150 (MQ=255)
                                            cAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGACAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCCGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAg                                                                                                         >  1:3108328/1‑150 (MQ=255)
                                            cAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAg                                                                                                         >  2:5585813/1‑150 (MQ=255)
                                            cAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAg                                                                                                         >  2:406002/1‑150 (MQ=255)
                                             aggaggGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCGCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGc                                                                                                        >  2:153404/1‑150 (MQ=255)
                                             aggaggGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGc                                                                                                        <  2:3953841/150‑1 (MQ=255)
                                             aggaggGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGc                                                                                                        >  2:745136/1‑150 (MQ=255)
                                               gaggGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCga                                                                                                      <  2:1261796/150‑1 (MQ=255)
                                                   gCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATc                                                                                                  <  2:161985/150‑1 (MQ=255)
                                                   gCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATc                                                                                                  <  1:429911/150‑1 (MQ=255)
                                                   gCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATc                                                                                                  >  1:2432767/1‑150 (MQ=255)
                                                    cGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCg                                                                                                 <  2:799344/150‑1 (MQ=255)
                                                     gAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGc                                                                                                <  1:4688052/150‑1 (MQ=255)
                                                      aGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCt                                                                                               <  2:4414374/150‑1 (MQ=255)
                                                      aGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCt                                                                                               <  2:5347935/150‑1 (MQ=255)
                                                           ggCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAAt                                                                                          <  2:3692540/150‑1 (MQ=255)
                                                           ggCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAAt                                                                                          <  2:1225157/150‑1 (MQ=255)
                                                             cGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGt                                                                                        <  2:3022558/150‑1 (MQ=255)
                                                             cGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGt                                                                                        <  2:1181630/150‑1 (MQ=255)
                                                               aTGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTaa                                                                                      <  1:4937909/150‑1 (MQ=255)
                                                               aTGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTaa                                                                                      <  2:2756235/150‑1 (MQ=255)
                                                               aTGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTaa                                                                                      <  1:3858439/150‑1 (MQ=255)
                                                               aTGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTaa                                                                                      <  2:5389574/150‑1 (MQ=255)
                                                                    tAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCAc                                                                                 >  2:3348618/1‑150 (MQ=255)
                                                                     aGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCACGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACt                                                                                <  2:4226080/150‑1 (MQ=255)
                                                                         tGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCa                                                                            <  2:96350/150‑1 (MQ=255)
                                                                           cGTGGGATAAGCGGTACAGCCGAGGAGGGATCTGCCATGTGCAGGATCCTTTTGCAAGAGCCCAGGGGCCGCCTAGCGGCCCCCGGTCATAATTACAGCACCACGCTTGGCAGCCTCAGCGAGATCGCTGGAATGTAAGTCACTGCCATc                                                                          <  2:4967638/150‑1 (MQ=255)
                                                                            gTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCa                                                                                                                <  1:1234457/111‑1 (MQ=255)
                                                                                gAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCAc                                                                     <  1:4390822/150‑1 (MQ=255)
                                                                                   aaGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAg                                                                  >  2:5074689/1‑150 (MQ=255)
                                                                                    agtggTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAgg                                                                 <  2:432476/147‑1 (MQ=255)
                                                                                      cGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGaa                                                               <  2:4548179/150‑1 (MQ=255)
                                                                                                 aggaggGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGt                                                    >  1:5217094/1‑150 (MQ=255)
                                                                                                  ggaggGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTa                                                   <  2:2222433/150‑1 (MQ=255)
                                                                                                     gggAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGaa                                                <  1:1417352/150‑1 (MQ=255)
                                                                                                        aaCAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACgg                                             >  1:351621/1‑150 (MQ=255)
                                                                                                        aaCAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCAGAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACgg                                             >  2:5026540/1‑150 (MQ=255)
                                                                                                          cAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCa                                           <  2:217212/150‑1 (MQ=255)
                                                                                                          cAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCa                                           <  1:467380/150‑1 (MQ=255)
                                                                                                           aGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCAt                                          >  1:2220141/1‑150 (MQ=255)
                                                                                                            gCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATc                                         <  1:5102229/150‑1 (MQ=255)
                                                                                                                 gtgCGGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATTGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAgcgc                                    >  1:158104/1‑150 (MQ=255)
                                                                                                                 gtgCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATTGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAgcgc                                    >  1:4128563/1‑150 (MQ=255)
                                                                                                                   gCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCt                                  >  2:3159670/1‑150 (MQ=255)
                                                                                                                    cAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCtt                                 >  2:4596206/1‑150 (MQ=255)
                                                                                                                      ggATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTACa                               >  2:1674986/1‑150 (MQ=255)
                                                                                                                       gATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCAc                              >  1:4964783/1‑150 (MQ=255)
                                                                                                                        aTCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCAcc                             <  1:1025147/150‑1 (MQ=255)
                                                                                                                        aTCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCAcc                             >  1:1449060/1‑150 (MQ=255)
                                                                                                                           cATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTc                          <  2:4704031/150‑1 (MQ=255)
                                                                                                                           cATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTc                          <  2:3167618/150‑1 (MQ=255)
                                                                                                                            aTTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCg                         <  1:484217/150‑1 (MQ=255)
                                                                                                                             tttGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGa                        >  1:3374169/1‑150 (MQ=255)
                                                                                                                             tttGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGa                        <  2:71208/150‑1 (MQ=255)
                                                                                                                                tcAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTc                     <  2:5561226/149‑1 (MQ=255)
                                                                                                                                gCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTc                     <  2:4651211/150‑1 (MQ=255)
                                                                                                                                gCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTc                     <  1:2172473/150‑1 (MQ=255)
                                                                                                                                  ttatgaccAGGGGCCGCCAAGCGGCCCCCGGTCATGATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGa                   >  1:1663321/7‑150 (MQ=255)
                                                                                                                                  aaGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGa                   >  2:884895/1‑150 (MQ=255)
                                                                                                                                  aaGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGa                   >  1:2733935/1‑150 (MQ=255)
                                                                                                                                    gagCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGATa                 <  2:2913211/150‑1 (MQ=255)
                                                                                                                                      gCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTaa                                                                                      >  1:2570944/1‑79 (MQ=255)
                                                                                                                                      gCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTaa                                                                                      <  2:2570944/79‑1 (MQ=255)
                                                                                                                                      gCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGATACt               <  1:4156011/150‑1 (MQ=255)
                                                                                                                                       cccAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGATACTc              >  2:5252674/1‑150 (MQ=255)
                                                                                                                                       cccAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGATACTc              >  1:627530/1‑150 (MQ=255)
                                                                                                                                        ccgggggCCGCTTGGCGGCCCCTGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTc                     <  1:3987790/139‑1 (MQ=255)
                                                                                                                                        ccgggggCCGCTTGGCGGCCCCTGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTc                     >  2:3987790/4‑142 (MQ=255)
                                                                                                                                        ccAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGATACTCa             >  2:5338177/1‑150 (MQ=255)
                                                                                                                                        ccAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGATACTCa             >  1:325569/1‑150 (MQ=255)
                                                                                                                                         cAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGATACTCAc            >  2:265389/1‑150 (MQ=255)
                                                                                                                                          gggggCCGCTTGGCGGCCCCTGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGATACTCAcc           >  2:2412275/2‑150 (MQ=255)
                                                                                                                                            gggCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTaa                                                                                      <  2:3755745/73‑1 (MQ=255)
                                                                                                                                               ccgccAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTaa                                                                                      <  1:2831833/70‑1 (MQ=255)
                                                                                                                                                 gccAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGAt                  >  1:4293475/1‑136 (MQ=255)
                                                                                                                                                   cAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGaa                                                               >  2:1760264/1‑89 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAGGCGATGCGCGCGGCGATGTGGTTGCCAGCGAATGAACAGGCCAGGAGGGCGAGCAGGGCGATGGCTAGTTTGCGTGGGAAAAGCGGTACAGCCGAGGAGGGAACAGCCATGTGCAGGATCCATTTGCAAGAGCCCAGGGGCCGCCAAGCGGCCCCCGGTCATAATTACAGCACCACGCTAGGCAGCCACAGCGAGATCGCTGGAATGTAAGTCACTGCCATCAGCACCAGGAACAGCGCCAGGTAGAACGGCATCAGCGCCTTCACCGTCGACTCGATACTCACCTTGCCAATG  >  NC_002947/1347059‑1347355

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