Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,628,429 A→C 5.4% A16A (GCA→GCC clpP → ATP‑dependent Clp protease proteolytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,628,4290AC 5.4% 372.8 / 8.3 129A16A (GCA→GCCclpPATP‑dependent Clp protease proteolytic subunit
Reads supporting (aligned to +/- strand):  ref base A (56/65);  new base C (5/2);  total (61/68)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 131 total aligned reads displayed.

CTGCTGAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAACTGCTG  >  NC_002947/2628280‑2628576
                                                                                                                                                     |                                                                                                                                                   
ctgctgAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCTCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCa                                                                                                                                                     >  2:930333/1‑150 (MQ=255)
ctgctgAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCa                                                                                                                                                     <  2:2365837/150‑1 (MQ=255)
 tgctgAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCGTCGAGCTGGCTTGTGCGTATTCAAGAAATGCGTATTTGGGACTGAGGGGAGTACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAg                                                                                                                                                    <  1:376114/150‑1 (MQ=255)
       gTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGcc                                                                                                                                              >  1:1767505/1‑150 (MQ=255)
        tAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCt                                                                                                                                             <  1:442478/150‑1 (MQ=255)
        tAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCt                                                                                                                                             >  1:2605983/1‑150 (MQ=255)
         aaCGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTg                                                                                                                                            <  1:2196041/150‑1 (MQ=255)
         aaCGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCCTAATTCTTATATTCAGCAGAGCCCTGACATCCAGGCCGCAGgccgccgc                                                                                                                                            <  2:1652265/150‑3 (MQ=255)
           cGGTCTTCTCTCGTAGGAAACCTCCACAAGCTAGCCGTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGATGTGTCAAGTAATTCTTCTATTCACCATAGCTCTGACATCCACGTCGCAGCCGGCTTGGt                                                                                                                                          >  2:1391082/1‑150 (MQ=255)
             gTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTcc                                                                                                                                        <  2:2057113/150‑1 (MQ=255)
              tcttctCTCGTAGGAAACCTCCACAAGCCAGCCTTTGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTccc                                                                                                                                       <  1:796464/150‑1 (MQ=255)
               cttctCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAg                                                                                                                                                    >  2:393604/1‑136 (MQ=255)
               cttctCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAg                                                                                                                                                    <  1:393604/136‑1 (MQ=255)
                 tctctcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCCCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGAt                                                                                                                                    >  2:3058444/1‑150 (MQ=255)
                  ctctcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATg                                                                                                                                   >  1:1289762/1‑150 (MQ=255)
                  ctctcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATg                                                                                                                                   <  1:1888666/150‑1 (MQ=255)
                   tctcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATgg                                                                                                                                  >  2:2850733/1‑150 (MQ=255)
                     tcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATg                                                                                                                                   <  2:1105594/147‑1 (MQ=255)
                     tcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATg                                                                                                                                   >  1:1105594/1‑147 (MQ=255)
                     tcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGtt                                                                                                                                >  1:2377069/1‑150 (MQ=255)
                        tAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATc                                                                                                                             <  1:1573526/150‑1 (MQ=255)
                           gAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAg                                                                                                                          <  1:1018406/150‑1 (MQ=255)
                             aaCCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCa                                                                                                                        <  2:467954/150‑1 (MQ=255)
                               ccTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGt                                                                                                                      <  2:2925304/150‑1 (MQ=255)
                                    acaAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGcc                                                                                                                 <  2:1926929/150‑1 (MQ=255)
                                        gccagccTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTCGGGAGAGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCCGTCCGCCCGTg                                                                                                             >  2:2962855/1‑150 (MQ=255)
                                         ccagccTTCGAGCGGGCTTGTGCGTGTTCATGCCATGACTATTTGGGAGTGTGCAGAGGACATGTCCCGGAATTCGTATATTCAGCAGAACGGTGACATCCAGGCCGCAGGTGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTgg                                                                                                            <  2:918314/150‑1 (MQ=255)
                                         ccagccTGCGAACTCGCTTGTGCGGATTCAAGTCAAGACGATTGGGGAGTGAGCGCAGAACATGACCCGCAAGTCGTATATTCAGCACAGCTCTGTCATCCAGGCCACGGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCgggg                                                                                                            <  2:2921708/150‑1 (MQ=255)
                                               ttCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAAcg                                                                                                      >  2:2873680/1‑150 (MQ=255)
                                                  gAGCGGGCTTGTGCGTACTCAAGCCAGGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACgcgc                                                                                                   <  2:2184053/150‑1 (MQ=255)
                                                   aGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCt                                                                                                  >  2:2174500/1‑150 (MQ=255)
                                                        gCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGc                                                                                                                  >  2:489928/1‑129 (MQ=255)
                                                        gCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGc                                                                                                                  <  1:489928/129‑1 (MQ=255)
                                                          ttGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATTCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCAGTGGCGAACGCGCTTACGACa                                                                                           >  2:184831/1‑150 (MQ=255)
                                                          ttGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGccc                                                                                                                >  1:1217482/1‑129 (MQ=255)
                                                          ttGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACa                                                                                           >  2:1925812/1‑150 (MQ=255)
                                                            gtgCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATc                                                                                         <  1:187580/150‑1 (MQ=255)
                                                             tgCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAg                                                                                                                          <  2:852868/116‑1 (MQ=255)
                                                             tgCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAg                                                                                                                          >  1:852868/1‑116 (MQ=255)
                                                                gTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACAACCAGGCCGCCTGCGGCCTggat                                                                                                                                         >  2:2256310/1‑96 (MQ=255)
                                                                gTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACAACCAGGCCGCCTGCGGCCTggat                                                                                                                                         <  1:2256310/98‑3 (MQ=255)
                                                                 tATTCAAGCCATGACTATTTGGGAGTGAGCGCAGAACCGCGCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTc                                                                                    <  1:2463867/150‑1 (MQ=255)
                                                                   ttCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTcgc                                                                                  <  1:2174500/150‑1 (MQ=255)
                                                                      aaGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACCTCTACTTGCTcc                                                                               >  1:2351095/1‑150 (MQ=255)
                                                                           aTGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTg                                                                                                                                            <  1:2905795/84‑1 (MQ=255)
                                                                           aTGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTg                                                                          >  2:458050/1‑150 (MQ=255)
                                                                             gccTATTTTGGTGTGAGCTCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                        <  2:220283/148‑1 (MQ=255)
                                                                             gACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                        <  1:2873680/150‑1 (MQ=255)
                                                                                ggtgtGGGGGTGTGCGCAGGACATGTCCCGCAATTCATATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGa                                                                     >  2:899765/5‑150 (MQ=255)
                                                                                  tttGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGc                                                                   <  2:1289762/150‑1 (MQ=255)
                                                                                    tGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCgt                                                                 <  2:1579303/150‑1 (MQ=255)
                                                                                       gAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAg                                                                       >  2:107955/1‑141 (MQ=255)
                                                                                       gAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAg                                                                       <  1:107955/141‑1 (MQ=255)
                                                                                        aGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGa                                                             <  2:2377069/150‑1 (MQ=255)
                                                                                          tGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATc                                                           >  2:1516494/1‑150 (MQ=255)
                                                                                            aGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCtt                                                         >  2:2370673/1‑150 (MQ=255)
                                                                                               gcAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTCTCGAGCAGTCCGCCCGTCGAGAACGCGCTTACGACACCTACTCGCTTCTGCTGACGGAGCGAGTGCTCTTCCt                                                      >  2:82508/1‑150 (MQ=255)
                                                                                                cAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTg                                                     <  1:1924558/150‑1 (MQ=255)
                                                                                                cAGGACATGTCCCGAAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGTGCCTGTTGAAGGAGCGTGTGATCTTCCTg                                                     <  1:1391082/150‑1 (MQ=255)
                                                                                                   gACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTg                                                                                                                                            >  1:1794667/1‑60 (MQ=255)
                                                                                                   gACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTg                                                                                                                                            <  2:1794667/60‑1 (MQ=255)
                                                                                                   gACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGtt                                                  <  1:2186102/150‑1 (MQ=255)
                                                                                                    aCATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTg                                                 <  1:1185649/150‑1 (MQ=255)
                                                                                                        gTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGccc                                             <  1:1516494/150‑1 (MQ=255)
                                                                                                         tCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGt                                                   >  2:2998242/1‑143 (MQ=255)
                                                                                                         tCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGt                                                   <  1:2998242/143‑1 (MQ=255)
                                                                                                         tCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCg                                            <  1:1477816/150‑1 (MQ=255)
                                                                                                         tCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCCTGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCg                                            >  1:806465/1‑150 (MQ=255)
                                                                                                          cccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCgg                                           <  2:2350533/150‑1 (MQ=255)
                                                                                                           ccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACa                                                                                           <  1:2488027/101‑1 (MQ=255)
                                                                                                           ccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTTCGATGGGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGt                                          >  1:808019/1‑150 (MQ=255)
                                                                                                           ccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGt                                          >  2:312959/1‑150 (MQ=255)
                                                                                                            cGCAATTCTTATATTCAGCAGAGCTCTTACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGccc                                             >  2:2298692/1‑146 (MQ=255)
                                                                                                            cGCAATTCTTATATTCAGCAGAGCTCTTACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGccc                                             <  1:2298692/146‑1 (MQ=255)
                                                                                                            cGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTa                                         >  1:1476477/1‑150 (MQ=255)
                                                                                                            cGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTa                                         >  1:998955/1‑150 (MQ=255)
                                                                                                               aaTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                      <  1:458050/150‑1 (MQ=255)
                                                                                                                    ttATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTa                                 >  1:2511153/1‑150 (MQ=255)
                                                                                                                         ttCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATGGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATgg                            >  2:1806484/1‑150 (MQ=255)
                                                                                                                         ttCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACa                                                                                           <  1:879967/87‑1 (MQ=255)
                                                                                                                         ttCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCgtg                                                                <  2:1601816/114‑1 (MQ=255)
                                                                                                                         ttCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCgtg                                                                >  1:1601816/1‑114 (MQ=255)
                                                                                                                         ttCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATgg                            <  2:2549655/150‑1 (MQ=255)
                                                                                                                          tCAGCAGAGCCCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGc                           <  1:3058444/150‑1 (MQ=255)
                                                                                                                             gcagAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCaa                        <  2:1767505/150‑1 (MQ=255)
                                                                                                                              cagAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                      <  2:2611744/135‑1 (MQ=255)
                                                                                                                              cagAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                      >  1:2611744/1‑135 (MQ=255)
                                                                                                                                      tcgggaccAGGCCGCCTGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTcgtcgt               >  2:1755917/7‑150 (MQ=255)
                                                                                                                                      tGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTcgtcgt               >  1:652605/1‑150 (MQ=255)
                                                                                                                                          aTCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgcg           >  1:828434/1‑150 (MQ=255)
                                                                                                                                           accAGGCCGCCTGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCgtg                                                                >  1:1369091/2‑96 (MQ=255)
                                                                                                                                           accAGGCCGCCTGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCgtg                                                                <  2:1369091/95‑1 (MQ=255)
                                                                                                                                            ccAGGCCGCCTGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCa         >  1:2562304/1‑150 (MQ=255)
                                                                                                                                            ccAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTAGCCCGGTAGAGGACTACATTGCCAACCTCGTCGTTGCGCa         >  2:2464317/1‑150 (MQ=255)
                                                                                                                                             cAGGCCCCCTGCGGCCTGGTCCCTAGGGTTCCCGAGCAGTGGGCGCGGGGCTAACTGGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGCGGACTACATGGCCAACCTCGTCGTTGCGCaa        >  2:507219/1‑150 (MQ=255)
                                                                                                                                              aGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAAc       >  2:2681664/1‑150 (MQ=255)
                                                                                                                                              aGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAAc       >  2:1924143/1‑150 (MQ=255)
                                                                                                                                                gCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAActg     <  2:2605983/150‑1 (MQ=255)
                                                                                                                                                   gCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAActgctg  >  1:1757990/1‑150 (MQ=255)
                                                                                                                                                   gCAGGCGGACTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGAAActgctg  <  1:1518280/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTGCTGAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAACTGCTG  >  NC_002947/2628280‑2628576

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