Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,999,172 Δ1 bp 16.7% intergenic (‑163/+37) PP_5252 ← / ← PP_5253 hypothetical protein/arylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,999,1720C.16.7% 557.4 / 133.5 210intergenic (‑163/+37)PP_5252/PP_5253hypothetical protein/arylesterase
Reads supporting (aligned to +/- strand):  ref base C (76/92);  new base . (17/18);  total (98/112)
Fisher's exact test for biased strand distribution p-value = 8.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-03

Only 100 of 211 total aligned reads displayed.

TAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCC  >  NC_002947/5999023‑5999320
                                                                                                                                                      |                                                                                                                                                    
 aGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCt                                                                                                                                                     <  1:4135959/150‑1 (MQ=255)
 aGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCt                                                                                                                                                     <  1:7628840/150‑1 (MQ=255)
    ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTa                                                                                                                                                  <  1:130756/150‑1 (MQ=255)
    ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTa                                                                                                                                                  <  1:2687115/150‑1 (MQ=255)
      ccGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACAAAGTATGTAgg                                                                                                                                                >  1:1378196/1‑150 (MQ=255)
       cGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGa                                                                                                                                               <  2:6504396/150‑1 (MQ=255)
         gCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGc                                                                                                                                             <  1:4419925/150‑1 (MQ=255)
         gCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGc                                                                                                                                             <  1:3203353/150‑1 (MQ=255)
         gCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGc                                                                                                                                             <  1:8779168/150‑1 (MQ=255)
           aGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCgg                                                                                                                                           <  2:3894734/150‑1 (MQ=255)
           aGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCgg                                                                                                                                           <  1:5142602/150‑1 (MQ=255)
           aGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCgg                                                                                                                                           <  1:2446760/150‑1 (MQ=255)
           aGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCgg                                                                                                                                           <  1:6649862/150‑1 (MQ=255)
           aGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGTCTGTAGGAGCgg                                                                                                                                           <  2:2998488/150‑1 (MQ=255)
           aGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCgg                                                                                                                                           <  2:8641628/150‑1 (MQ=255)
           aGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCgg                                                                                                                                           <  1:8425612/150‑1 (MQ=255)
           aGGGGGACAGCGAAGCATGCCGGTTGATGCTTGGCCTGGGCTTGAATGGATGTGCTTTGTGTTGAAGAAGGTTAGCCCTGGTCATGCTTTTTTCTGACTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGTCTTTGTTGGAGCgg                                                                                                                                           <  2:1443423/150‑1 (MQ=37)
            gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCggg                                                                                                                                          <  1:7080411/150‑1 (MQ=255)
            gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTATTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCggg                                                                                                                                          <  2:2699023/150‑1 (MQ=255)
                   tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCGTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTAccc                                                                                                                                   >  2:1042004/1‑150 (MQ=255)
                   tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTAccc                                                                                                                                   >  2:8590944/1‑150 (MQ=255)
                   tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTATTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCgg                                                                                                                                           <  2:2205865/142‑1 (MQ=255)
                   tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTATTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCgg                                                                                                                                           >  1:2205865/1‑142 (MQ=255)
                           aTGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCACGTCTGTAGGATCGGGTTTACCCGCGAATGt                                                                                                                           >  2:2931655/1‑150 (MQ=255)
                               cGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTcc                                                                                                                       <  1:8181773/150‑1 (MQ=255)
                                 gTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCCCGAAAGTTTcccc                                                                                                                     >  1:752537/1‑149 (MQ=255)
                                         aaGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTAcc                                                                                                             <  2:504810/150‑1 (MQ=255)
                                          aGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCg                                                                                                            <  2:4055498/150‑1 (MQ=255)
                                            gCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAAATAACGCt                                                                                                          >  1:6596818/1‑150 (MQ=255)
                                            gCCGGGGCTTGAATGGAGGGGGGGGGGGTTGAAAAAGGTTCGCCCGGGTCAAGCTTTTTGCTGCCGGGGCCGGCCTCTTCGGGGGTAAACCCGCCCCTACACGGTCTGTCGGAGGGGGTTTACCCGCGAATGTTTCCCACAATTACCGCt                                                                                                          <  1:4948039/150‑1 (MQ=255)
                                              cTGGGCTTTAATGTATGTGCGGTGGGTTGAAGAAGGTGCGCCCTGGCCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGTCTGTCGGAGGGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGa                                                                                                        <  1:3154647/150‑1 (MQ=38)
                                              cTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGa                                                                                                        <  2:5917524/150‑1 (MQ=255)
                                              cTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGa                                                                                                        <  1:1500842/150‑1 (MQ=255)
                                               tGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTAACGCTGaa                                                                                                       >  2:1611392/1‑150 (MQ=255)
                                                      gattGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGc                                                                                                <  2:8466164/147‑1 (MQ=255)
                                                      gAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGc                                                                                                <  2:5470328/150‑1 (MQ=255)
                                                        aTGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTTAAACCCCGCCAAAAACGTCTGTTTGAGCGGGTTTTAACGCCAATGTTTACCACAAATAACGCGGAAGGAAAGCCa                                                                                              >  1:1683954/1‑150 (MQ=255)
                                                            aTGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa                                                                                          <  2:6097373/150‑1 (MQ=255)
                                                            aTGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa                                                                                          <  2:4031767/150‑1 (MQ=255)
                                                            aTGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa                                                                                          <  2:2532178/150‑1 (MQ=255)
                                                            aTGTGCGGGGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTGTTCGGGGTTAAACGCGCTCCTACAACGTGTGTAGGTGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa                                                                                          <  1:1111289/150‑1 (MQ=255)
                                                               tgCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTc                                                                                       <  2:7006678/150‑1 (MQ=255)
                                                                 gggtggtTTGAATAAGGTTAGCCCTGGTCAAGATTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCAACGACTGTAGGAGCGTGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCt                                                                                     <  2:7651560/149‑1 (MQ=38)
                                                                  ggTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTc                                                                                    <  2:5110989/150‑1 (MQ=255)
                                                                       gTTGAAGAAGTTTAGCCCTGGTCAAGCTTTTTGCGGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGTGTTTACCCGCGAATGTTTCCCACTATTACCGCTGAAGGAACGCCAGCAAGTCCTAATTCa                                                                               <  2:6628039/150‑1 (MQ=38)
                                                                        ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAg                                                                              <  2:3606183/150‑1 (MQ=255)
                                                                        ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCGGAAGGAACGCCAGCAAGTCCTCATTCAg                                                                              >  2:7464750/1‑150 (MQ=255)
                                                                        ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAAAACCTCGGTAGGAGCGGGTTTACCCGCGAATGTTTACCACAAATACCGCTGAAGGAACGCCAGCAAGTACTCATTCAg                                                                              >  2:865273/1‑150 (MQ=38)
                                                                             gaagGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTTAAACCCCACCTAACACGGCGGGAGGAGCGGGTTTTCCCGCGAATGTTTCCCACCAGTTCCGCTGAAGGGACGCCCGCAAGTACTCATTCAgctgct                                                                         >  1:3489611/1‑150 (MQ=255)
                                                                                  ttAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCa                                                                    <  1:8590944/150‑1 (MQ=255)
                                                                                   tAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCAt                                                                   <  2:7569695/150‑1 (MQ=255)
                                                                                     gCCCTGGTCATGTTTTTTGCTGCCTGGGCCGGCCTCTTGGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGTGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGa                                                                 <  1:141690/150‑1 (MQ=255)
                                                                                     gCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGa                                                                 <  1:6453829/150‑1 (MQ=255)
                                                                                      cccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAg                                                                >  2:9060501/1‑150 (MQ=255)
                                                                                             cAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCCCGCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACCCCAGCAAGTCCTCATTCAGCAGCTGCGCATGAGTCACAGc                                                         >  2:3389622/1‑150 (MQ=38)
                                                                                              aaGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCCCTCCTACAACGTATGGAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTAACCCGGAATGCACGCCAGCAAGTACGCAATCAGCTGCTGCGCCAGAGTCAATGCa                                                        >  2:7698483/1‑150 (MQ=25)
                                                                                               aGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCAATCAGCTGCTGCGCATGAGTCACTGCaa                                                       >  2:6815962/1‑150 (MQ=255)
                                                                                                       gctgcCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                               <  2:1242452/150‑1 (MQ=255)
                                                                                                       gctgcCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAAAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                               >  2:537535/1‑150 (MQ=255)
                                                                                                          gcCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgcggc                                            <  1:4724872/150‑1 (MQ=255)
                                                                                                          gcCTGGGCCGGCCTCTTCGCGGGAAAACCCGCTCCTACAACGTCTGTCGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgcggc                                            <  1:3843407/150‑1 (MQ=255)
                                                                                                              gggCCGGCCTCTTCGCGGGTAACCCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGccc                                        <  1:1205140/150‑1 (MQ=255)
                                                                                                                   ggccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCa                                   <  1:2119179/150‑1 (MQ=255)
                                                                                                                    gccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCAt                                  >  2:3754665/1‑150 (MQ=255)
                                                                                                                     ccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATa                                 <  1:6058038/150‑1 (MQ=255)
                                                                                                                     ccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATa                                 <  2:1787144/150‑1 (MQ=255)
                                                                                                                         ttCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCCTACAcc                             >  2:5759169/1‑150 (MQ=255)
                                                                                                                           cgcgGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCATGCGGCGCCCCGGCATACACCtt                           <  2:2143771/150‑1 (MQ=255)
                                                                                                                              gggTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAg                        <  2:539048/150‑1 (MQ=255)
                                                                                                                                gTAAACCCGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCt                      >  1:4457774/1‑150 (MQ=255)
                                                                                                                                gTAAACCCGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCt                      >  1:6247461/1‑150 (MQ=255)
                                                                                                                                gTAAACCCGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCt                      >  1:2328771/1‑150 (MQ=255)
                                                                                                                                gTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCt                      >  2:6726967/1‑150 (MQ=255)
                                                                                                                                    aCCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc                  >  1:2974176/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg                 >  1:4780240/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg                 >  1:2680476/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg                 >  1:7708464/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg                 >  1:631609/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcgc                >  1:7966929/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGTGGCGCCCCGGCATACACCTTCAGCTCGgcgc                >  1:2842074/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcgc                <  2:714588/150‑1 (MQ=255)
                                                                                                                                       cGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGcc               >  2:1939736/1‑150 (MQ=255)
                                                                                                                                        gCTCCTACCACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCa              <  2:4789323/150‑1 (MQ=255)
                                                                                                                                        gCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCa              >  1:940205/1‑150 (MQ=255)
                                                                                                                                        gCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCa              >  1:6863179/1‑150 (MQ=255)
                                                                                                                                        gCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCa              >  1:487514/1‑150 (MQ=255)
                                                                                                                                         ctcctcCAACGTCTGTAGGAGCGGGATTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCAc             <  2:461249/150‑1 (MQ=255)
                                                                                                                                             tACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAt         <  1:9060501/150‑1 (MQ=255)
                                                                                                                                              aCAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATc        <  2:8905788/150‑1 (MQ=255)
                                                                                                                                              aCAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATc        <  1:4015097/150‑1 (MQ=255)
                                                                                                                                              aCAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATc        >  1:5036961/1‑150 (MQ=255)
                                                                                                                                              aCAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATc        >  2:5899971/1‑150 (MQ=255)
                                                                                                                                              aCAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATc        >  1:1570884/1‑150 (MQ=255)
                                                                                                                                              aCAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATc        >  2:1597874/1‑150 (MQ=255)
                                                                                                                                               cAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCa       >  2:8739046/1‑150 (MQ=255)
                                                                                                                                                aaCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGTGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAg      <  2:2842074/150‑1 (MQ=255)
                                                                                                                                                aaCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAg      >  1:5776231/1‑150 (MQ=255)
                                                                                                                                                    gTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTc   >  2:256950/1‑150 (MQ=255)
                                                                                                                                                     tCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTcc  >  1:8149351/1‑150 (MQ=255)
                                                                                                                                                     tCTGTAGGAGCGGGATTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTcc  <  2:9113537/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCC  >  NC_002947/5999023‑5999320

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