Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,564,981 C→A 5.1% intergenic (+14/‑50) fepA → / → PP_2243 ferric enterobactin transport system outer membrane subunit/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,564,9810CA 5.1% 1237.7 / 37.1 364intergenic (+14/‑50)fepA/PP_2243ferric enterobactin transport system outer membrane subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (171/174);  new base A (8/11);  total (179/185)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 380 total aligned reads displayed.

CTACGACGTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCC  >  NC_002947/2564832‑2565129
                                                                                                                                                     |                                                                                                                                                    
cTACGACGTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGgccgcc                                                                                                                                                      <  2:2523697/150‑1 (MQ=255)
    gacgTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGc                                                                                                                                                  >  2:1497463/1‑150 (MQ=255)
       gTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGc                                                                                                                                               <  1:7580976/150‑1 (MQ=255)
           ggggAGGCGTGAACACCGTCTTCGACAAGCAGATCCTGCGCGGCGGCACCGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCAGGGCGGccctt                                                                                                                                           <  2:8456594/150‑3 (MQ=255)
            ggTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCACGGGCCGCAGGGCGGcccttg                                                                                                                                          <  1:1469095/150‑4 (MQ=255)
              tAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCAGGGCGGcccttggg                                                                                                                                        <  2:1263684/150‑6 (MQ=255)
                ggCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATcgcg                                                                                                                                      <  2:63572/150‑1 (MQ=255)
                 gCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCgg                                                                                                                                     >  2:5068728/1‑150 (MQ=255)
                 gCGTGAACAACGTCTTAGACAAGCAGATCCTGCGCGGCGGCAACGACAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATAGCgg                                                                                                                                     >  2:1880270/1‑150 (MQ=255)
                   gTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGtt                                                                                                                                   <  2:3880735/150‑1 (MQ=255)
                       acaacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAAt                                                                                                                               >  2:8376791/1‑150 (MQ=255)
                       acaacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAAt                                                                                                                               <  2:5181655/150‑1 (MQ=255)
                        caacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATc                                                                                                                              >  1:1052336/1‑150 (MQ=255)
                        caacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATc                                                                                                                              <  2:5703170/150‑1 (MQ=255)
                         aacGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCg                                                                                                                             >  1:1873586/1‑150 (MQ=255)
                           cGTCTTCGACAAGCAGATCCTGCGCGGCGGCTACGCAAGCAGCTCGGGGGCAAACACTTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCggg                                                                                                                           <  1:7641726/150‑1 (MQ=255)
                           cGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCggg                                                                                                                           <  1:5191441/150‑1 (MQ=255)
                           cGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATAACAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCggg                                                                                                                           >  1:3971270/1‑150 (MQ=255)
                                  gACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCACTCGCGGTTCAATCGGGTCGCGaa                                                                                                                    >  1:6289385/1‑150 (MQ=255)
                                   aCAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAAt                                                                                                                   <  1:4463152/150‑1 (MQ=255)
                                    cAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATg                                                                                                                  <  2:4904886/150‑1 (MQ=255)
                                    cAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATg                                                                                                                  <  1:5544886/150‑1 (MQ=255)
                                     aaGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATgg                                                                                                                 >  2:3414263/1‑150 (MQ=255)
                                      aGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGc                                                                                                                <  2:2608636/150‑1 (MQ=255)
                                       gCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGcc                                                                                                               >  1:2987303/1‑150 (MQ=255)
                                       gCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGcc                                                                                                               >  1:7962535/1‑150 (MQ=255)
                                          gATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCcgac                                                                                                            >  2:2490259/1‑150 (MQ=255)
                                            tCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCcgacga                                                                                                          <  1:1635771/150‑1 (MQ=255)
                                                   cggcggCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCtt                                                                                                   <  2:1114111/150‑1 (MQ=255)
                                                   cggcggCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCtt                                                                                                   <  2:5497512/150‑1 (MQ=255)
                                                     gcggcAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGc                                                                                                 <  2:4828383/150‑1 (MQ=255)
                                                      cggcAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCa                                                                                                <  2:2011253/150‑1 (MQ=255)
                                                      cggcAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCa                                                                                                <  2:3589288/150‑1 (MQ=255)
                                                       ggcAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAg                                                                                               <  1:6646678/150‑1 (MQ=255)
                                                           aCGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCaa                                                                                           >  2:726083/1‑150 (MQ=255)
                                                            cGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCaaa                                                                                          <  1:1211359/150‑1 (MQ=255)
                                                             gCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAAc                                                                                         >  1:6775814/1‑150 (MQ=255)
                                                               cagcagCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCcga                                                                                                             >  2:7613821/1‑128 (MQ=255)
                                                               cagcagCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGt                                                                                       >  1:4590639/1‑150 (MQ=255)
                                                                 gcagcTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGtct                                                                                     <  2:1654922/150‑1 (MQ=255)
                                                                 gcagcTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGtct                                                                                     >  1:5330183/1‑150 (MQ=255)
                                                                  cagcTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGtctc                                                                                    <  2:6464814/150‑1 (MQ=255)
                                                                   agcTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCg                                                                                   >  1:8309017/1‑150 (MQ=255)
                                                                      tCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAg                                                                                >  1:3115723/1‑150 (MQ=255)
                                                                       cGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGc                                                                               <  1:3668451/150‑1 (MQ=255)
                                                                       cGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGc                                                                               >  1:8189053/1‑150 (MQ=255)
                                                                       cGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGc                                                                               >  1:6098209/1‑150 (MQ=255)
                                                                         ggggCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGt                                                                             <  1:3813759/150‑1 (MQ=255)
                                                                          gggCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTg                                                                            <  1:2013386/150‑1 (MQ=255)
                                                                          gggCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTg                                                                            >  1:5268333/1‑150 (MQ=255)
                                                                              aaaCACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTaa                                                                        >  2:7983828/1‑150 (MQ=255)
                                                                              aaaCACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTaa                                                                        >  2:2825085/1‑150 (MQ=255)
                                                                               aaCACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTaaa                                                                       >  1:185937/1‑150 (MQ=255)
                                                                                acacGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAAc                                                                      >  1:2730492/1‑150 (MQ=255)
                                                                                 cacGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAAcc                                                                     <  1:494925/150‑1 (MQ=255)
                                                                                   cGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCAt                                                                   <  1:7148887/150‑1 (MQ=255)
                                                                                    gTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATg                                                                  >  2:379782/1‑150 (MQ=255)
                                                                                    gTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATg                                                                  <  2:3484759/150‑1 (MQ=255)
                                                                                      ataACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCa                                                                >  2:8736189/1‑150 (MQ=255)
                                                                                        aaCCAGCCGGGCAGTGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGTCCGCCCTGCGGCCCAATCGCGGTTCAATAGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCAcc                                                              >  2:476369/1‑150 (MQ=255)
                                                                                        aaCCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCAcc                                                              >  2:7465218/1‑150 (MQ=255)
                                                                                          ccAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCAcccc                                                            >  2:8092400/1‑150 (MQ=255)
                                                                                          ccAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCAcccc                                                            >  1:2336500/1‑150 (MQ=255)
                                                                                          ccAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCAcccc                                                            >  1:6987707/1‑150 (MQ=255)
                                                                                            aGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAg                                                          <  1:4631878/150‑1 (MQ=255)
                                                                                                  gCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCt                                                    >  1:4758768/1‑150 (MQ=255)
                                                                                                    aGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGc                                                  >  2:4712264/1‑150 (MQ=255)
                                                                                                      ggCGGTGTGTGCTGCGCTGAATATCAACGTCGGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCCGCga                                                <  1:2624234/150‑1 (MQ=255)
                                                                                                         ggTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGaa                                             <  1:58794/150‑1 (MQ=255)
                                                                                                             tttGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGc                                         <  1:6890747/150‑1 (MQ=255)
                                                                                                               tgctgcGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCt                                       <  2:8385577/150‑1 (MQ=255)
                                                                                                                gctgcGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTg                                      <  1:8699463/150‑1 (MQ=255)
                                                                                                                 ctgcGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTgg                                     <  1:2884648/150‑1 (MQ=255)
                                                                                                                  tgcGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGt                                    >  1:1693137/1‑150 (MQ=255)
                                                                                                                      cTGAATATCAACTTCTGATCCCACGGGCCGCCCTGCGGCCCAATCGCGGGTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGaaa                                <  1:1033644/150‑1 (MQ=255)
                                                                                                                      cTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGaaa                                <  2:6760091/150‑1 (MQ=255)
                                                                                                                       tGAATATCAACTTCTGATCCCAAGTGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAg                               >  2:7144689/1‑150 (MQ=255)
                                                                                                                            atCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGGTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATc                          >  1:994119/1‑150 (MQ=255)
                                                                                                                             tgtaCTTCTGTTCCGAAGGTCAGCCCAGCGGCCCAATCGCGGTTGAATCTGTTCTCTAAAGGCAGACGACAATCTTGAAGGCAAACTTCTCGCAGCGTGCTAATCCATGCAACCAAGCCCGCTGCGAGAAATGCCTGGTGAATGCCaaca                         >  2:3054414/4‑147 (MQ=255)
                                                                                                                             tCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATcc                         >  2:6701210/1‑150 (MQ=255)
                                                                                                                             tCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATcc                         >  1:5465901/1‑150 (MQ=255)
                                                                                                                             tCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATcc                         >  2:2720569/1‑150 (MQ=255)
                                                                                                                             tCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACACCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATcc                         >  2:437645/1‑150 (MQ=255)
                                                                                                                              cAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCg                        >  2:6136902/1‑150 (MQ=255)
                                                                                                                               aaCTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGc                       >  1:2137969/1‑150 (MQ=255)
                                                                                                                               aaCTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGc                       >  1:773932/1‑150 (MQ=255)
                                                                                                                                aCTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGcc                      <  2:2046523/150‑1 (MQ=255)
                                                                                                                                aCTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGcc                      >  2:326069/1‑150 (MQ=255)
                                                                                                                                 cTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCCTCCGCCt                     >  2:5264110/1‑150 (MQ=255)
                                                                                                                                      gATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGAccc                >  2:6905487/1‑150 (MQ=255)
                                                                                                                                       cgcgattgggCCGCAGGGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCt               >  2:8503137/8‑150 (MQ=255)
                                                                                                                                        tCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTg              >  2:7834764/1‑150 (MQ=255)
                                                                                                                                          gattgggCCGCAGGGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGAt            >  1:92670/5‑150 (MQ=255)
                                                                                                                                          ccAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGAt            >  1:2389651/1‑150 (MQ=255)
                                                                                                                                           cAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATc           >  2:218686/1‑150 (MQ=255)
                                                                                                                                           cAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATc           >  2:1443764/1‑150 (MQ=255)
                                                                                                                                            aaGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCt          >  2:6791385/1‑150 (MQ=255)
                                                                                                                                            aaGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCt          <  2:5957459/150‑1 (MQ=255)
                                                                                                                                                  cgccCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATg    >  1:7282922/1‑150 (MQ=255)
                                                                                                                                                    cccTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGcc  <  1:846185/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTACGACGTGCGGGTAGGCGTGAACAACGTCTTCGACAAGCAGATCCTGCGCGGCGGCAACGCCAGCAGCTCGGGGGCAAACACGTATAACCAGCCGGGCAGGGCGGTGTTTGCTGCGCTGAATATCAACTTCTGATCCCAAGGGCCGCCCTGCGGCCCAATCGCGGTTCAATCGGGTCGCGAATGGCCGACGACATTCTTGCAGCCAAACGTCTCGCAGCGTGCTAAACCATGCACCCCAGCCCGCTGCGAGAAATGCCTGGTGAAAGCCATCCGCCTGACCCTGATCTGCCATGCC  >  NC_002947/2564832‑2565129

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