Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,651,319 G→C 8.4% intergenic (+41/+29) PP_2323 → / ← aroF‑I hypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,651,3190GC 8.4% 254.7 / 13.4 83intergenic (+41/+29)PP_2323/aroF‑Ihypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (39/37);  new base C (3/4);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

GGCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACTGCAAATC  >  NC_002947/2651176‑2651467
                                                                                                                                                |                                                                                                                                                    
ggCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCg                                                                                                                                                <  1:1942672/150‑1 (MQ=255)
 gCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCgg                                                                                                                                               >  1:1981007/1‑150 (MQ=255)
  ctctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGc                                                                                                                                              <  1:1709699/150‑1 (MQ=255)
   tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCg                                                                                                                                             <  2:1101760/150‑1 (MQ=255)
   tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCg                                                                                                                                             <  2:1100629/150‑1 (MQ=255)
   tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATAACGTCTTGCTCATATGTGGAACGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCAACAGACACAACAAAAAACGCCGGGCAGTGCACGGCg                                                                                                                                             >  1:362910/1‑150 (MQ=255)
       gAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGtttt                                                                                                                                         >  2:321085/1‑150 (MQ=255)
        aGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGttttt                                                                                                                                        <  2:2115697/150‑1 (MQ=255)
            cGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCAtt                                                                                                                                    >  1:2007117/1‑150 (MQ=255)
             gCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCGCCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTc                                                                                                                                   >  2:1550128/1‑150 (MQ=255)
             gCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTc                                                                                                                                   >  2:179812/1‑150 (MQ=255)
                       tGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGtt                                                                                                                         >  1:2447451/1‑150 (MQ=255)
                                   tGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTa                                                                                                             <  1:1524123/150‑1 (MQ=255)
                                   tGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTa                                                                                                             <  2:963132/150‑1 (MQ=255)
                                        ccGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTg                                                                                                        >  2:1067718/1‑150 (MQ=255)
                                        ccGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTg                                                                                                        >  2:1048684/1‑150 (MQ=255)
                                          gTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCg                                                                                                      <  2:2007117/150‑1 (MQ=255)
                                             ttGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgca                                                                                                   <  1:804426/150‑1 (MQ=255)
                                                   atatGTGGAGCGTTACATTGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACg                                                                                                           <  1:2060989/136‑1 (MQ=255)
                                                   atatGTGGAGCGTTACATTGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACg                                                                                                           >  2:2060989/1‑136 (MQ=255)
                                                    tatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGcaca                                                                                              <  1:1312868/148‑1 (MQ=255)
                                                    tatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGcaca                                                                                              >  2:1312868/1‑148 (MQ=255)
                                                       gtgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                         <  2:2071541/150‑1 (MQ=255)
                                                       gtgCAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                         <  1:2453243/150‑1 (MQ=255)
                                                             cGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGATCGTGTTGGTCCGGCGTTTTTCATTCCCGCAATGTTCAGACGCGCTTACGCTGCGGCAGCACATCCTTGAGCtt                                                                                   >  1:42721/1‑150 (MQ=255)
                                                               ttACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTgg                                                                                 >  2:579270/1‑150 (MQ=255)
                                                               ttACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTgg                                                                                 <  1:1583295/150‑1 (MQ=255)
                                                                  cATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcat                                                                              >  2:952158/1‑150 (MQ=255)
                                                                   aTGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGt                                                                                                                          >  1:1036175/1‑105 (MQ=255)
                                                                   aTGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGt                                                                                                                          <  2:1036175/105‑1 (MQ=255)
                                                                       aGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGcatgcatg                                                                         <  2:1981007/150‑1 (MQ=255)
                                                                        gcgcCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACCGGGACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATTCTTGAGCTTGGCatgcatgc                                                                        >  2:2169207/1‑150 (MQ=255)
                                                                        gcgcCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTg                                                                                  >  1:391908/1‑140 (MQ=255)
                                                                        gcgcCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTg                                                                                  <  2:391908/140‑1 (MQ=255)
                                                                           ccAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCAtgctgc                                                                     >  1:221504/1‑150 (MQ=255)
                                                                              gcagcGGCATTCGCGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCTGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCa                                                                  <  2:1323240/150‑1 (MQ=255)
                                                                              gcagcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCATCACCTCCTGGAGCTTGGCGTGCATGCTGCGCa                                                                  >  2:124359/1‑150 (MQ=255)
                                                                               cagcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAg                                                                 <  2:221504/150‑1 (MQ=255)
                                                                                agcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAgg                                                                <  1:321085/150‑1 (MQ=255)
                                                                                agcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAgg                                                                >  2:1548081/1‑150 (MQ=255)
                                                                                        tCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtctc                                                        <  2:207207/150‑1 (MQ=255)
                                                                                        tCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtctc                                                        <  2:2227788/150‑1 (MQ=255)
                                                                                          gAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCgg                                                      >  2:1589887/1‑150 (MQ=255)
                                                                                           aGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatgcaag                                                                         <  2:1174500/130‑3 (MQ=255)
                                                                                           aGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGcatgcatg                                                                         >  1:1174500/1‑130 (MQ=255)
                                                                                            gCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGATCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTa                                                    <  1:2017222/150‑1 (MQ=255)
                                                                                            gCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTa                                                    <  1:1589887/150‑1 (MQ=255)
                                                                                             cAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAg                                                   >  2:2220945/1‑150 (MQ=255)
                                                                                                 gCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGa                                               >  1:1432387/1‑150 (MQ=255)
                                                                                                  ccTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATGCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGAc                                              <  2:42721/150‑1 (MQ=255)
                                                                                                    tGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGTGCACTGCCCGGCGtt                                                                                                                                           <  1:482073/55‑1 (MQ=255)
                                                                                                    tGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCa                                            <  1:1067718/150‑1 (MQ=255)
                                                                                                    tGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGtt                                                                                                                                           >  2:482073/1‑55 (MQ=255)
                                                                                                         tGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGa                                       <  1:1548081/150‑1 (MQ=255)
                                                                                                           acgacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTc                                         >  1:654946/1‑146 (MQ=255)
                                                                                                           acgacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTc                                         >  1:654667/1‑146 (MQ=255)
                                                                                                           acgacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTc                                         <  2:654667/146‑1 (MQ=255)
                                                                                                           acgacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTc                                         <  2:654946/146‑1 (MQ=255)
                                                                                                           acgacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGCTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATg                                     >  2:49843/1‑150 (MQ=255)
                                                                                                              acAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                                  <  1:1431912/150‑1 (MQ=255)
                                                                                                                aGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAg                                                                                      >  1:313790/1‑96 (MQ=255)
                                                                                                                aGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAg                                                                                      <  2:313790/96‑1 (MQ=255)
                                                                                                                     ccAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCg                           >  1:132332/1‑150 (MQ=255)
                                                                                                                      cAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGt                          >  1:393213/1‑150 (MQ=255)
                                                                                                                        gACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGa                        <  1:179812/150‑1 (MQ=255)
                                                                                                                           acaacaAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCg                     >  2:1767873/1‑150 (MQ=255)
                                                                                                                             aacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                                  <  1:1054464/135‑1 (MQ=255)
                                                                                                                             aacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                                  >  2:1054464/1‑135 (MQ=255)
                                                                                                                              acaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatgcat                                                                          <  1:784664/94‑1 (MQ=255)
                                                                                                                              acaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatgcat                                                                          >  2:784664/1‑94 (MQ=255)
                                                                                                                               caaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                         >  1:993751/1‑78 (MQ=255)
                                                                                                                               caaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                         <  2:993751/78‑1 (MQ=255)
                                                                                                                                aaaaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACg                <  1:579270/150‑1 (MQ=255)
                                                                                                                                 cggaccGCCGGGCCGGGCCGGGCGTTGTGCAGCCCCGCAGCGTAAAGCCGCGCTTACGCTGCGGCAGGACACCCTGGAGCTGGGCCTGCCTGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGGCACGc               <  2:558671/145‑1 (MQ=255)
                                                                                                                                 aaaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCa                                   <  1:997393/130‑1 (MQ=255)
                                                                                                                                 aaaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCa                                   >  2:997393/1‑130 (MQ=255)
                                                                                                                                  aaaaCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGcc              >  2:1129257/1‑150 (MQ=255)
                                                                                                                                      cGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATAc          >  1:193314/1‑150 (MQ=255)
                                                                                                                                      cGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGt                              <  2:2185108/130‑1 (MQ=255)
                                                                                                                                      cGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGt                              >  1:2185108/1‑130 (MQ=255)
                                                                                                                                       gCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACt         >  1:2422020/1‑150 (MQ=255)
                                                                                                                                             gCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACTGCaaa    >  1:2172319/1‑150 (MQ=255)
                                                                                                                                               aGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACTGCAAATc  >  1:2361328/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATACTGCAAATC  >  NC_002947/2651176‑2651467

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