Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,559,450 T→C 8.9% D570D (GAT→GAC) ‡ fepA → ferric enterobactin transport system outer membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,559,4500TC 8.9% 181.8 / 7.5 67D570D (GAT→GAC) ‡fepAferric enterobactin transport system outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base C (3/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

AATGGTGGCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGCTCG  >  NC_002947/2559304‑2559598
                                                                                                                                                   |                                                                                                                                                     
aaTGGTGGCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATAc                                                                                                                                                     <  1:436847/149‑1 (MQ=255)
    gtggCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGAt                                                                                                                                                 >  1:180231/1‑149 (MQ=255)
     tggCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTATCACCGAtg                                                                                                                                                <  2:437876/149‑1 (MQ=255)
       gCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCCGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTATCACCGAtga                                                                                                                                               <  7:137318/148‑2 (MQ=255)
        cTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGACACAAGATCATCGGTGATACCGAtgtg                                                                                                                                              <  1:3615/149‑1 (MQ=255)
        cTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTATCACCGAtgatc                                                                                                                                             <  6:169058/149‑4 (MQ=255)
             aCCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTa                                                                                                                                         <  4:379369/148‑1 (MQ=255)
              ccTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTac                                                                                                                                        <  2:230226/148‑1 (MQ=255)
               cTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTacac                                                                                                                                      <  5:432701/149‑1 (MQ=255)
                  cAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCa                                                                                                                                    <  2:316005/148‑1 (MQ=255)
                  cAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCAt                                                                                                                                   >  5:210853/1‑149 (MQ=255)
                      gCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTgg                                                                                                                                >  3:70393/1‑148 (MQ=255)
                      gCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTgg                                                                                                                                <  4:70393/148‑1 (MQ=255)
                      gCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGt                                                                                                                               <  1:362576/149‑1 (MQ=255)
                                 cTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCcg                                                                                                                   >  5:378407/1‑149 (MQ=255)
                                  tCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCcgc                                                                                                                  >  6:263810/1‑149 (MQ=255)
                                    aaGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGt                                                                                                                >  1:288731/1‑149 (MQ=255)
                                       cccGAGATCAGCGTCAACAAGGAAATCGTCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATTTGCTTTACAGCATTGGTACCGGCAGCCGCGTGAc                                                                                                             >  5:93127/1‑149 (MQ=255)
                                        ccGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCCCAACGACTACAATAACAAGATCATCGGTGATACCGATGTCCTTTACACCATTGGTACCGGCAGCCCCGTGACg                                                                                                            >  7:333851/1‑149 (MQ=255)
                                             aTCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTg                                                                                                       <  5:263810/149‑1 (MQ=255)
                                             aTCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACAGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTg                                                                                                       <  7:291376/149‑1 (MQ=255)
                                              tCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTgg                                                                                                      >  1:422238/1‑149 (MQ=255)
                                               cAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTgg                                                                                                      >  7:156024/1‑148 (MQ=255)
                                                aGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTg                                                                                                                                 <  7:276407/121‑1 (MQ=255)
                                                aGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTg                                                                                                                                 >  8:276407/1‑121 (MQ=255)
                                                   gTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACaa                                                                                                 <  2:422238/149‑1 (MQ=255)
                                                         aaGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCgg                                                                                           >  6:75621/1‑149 (MQ=255)
                                                           ggAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGc                                                                                          >  6:264437/1‑148 (MQ=255)
                                                                    ccTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCggg                                                                                >  6:152266/1‑149 (MQ=255)
                                                                    ccTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGAGCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCggg                                                                                <  5:36539/149‑1 (MQ=255)
                                                                      tgctgTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTg                                                                              >  4:83274/1‑149 (MQ=255)
                                                                           tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAgg                                                                         >  3:425269/1‑149 (MQ=255)
                                                                           tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAgg                                                                         >  1:205553/1‑149 (MQ=255)
                                                                           tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAg                                                                          <  8:156024/148‑1 (MQ=255)
                                                                           tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAg                                                                          <  5:264437/148‑1 (MQ=255)
                                                                           tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAg                                                                          <  6:378407/148‑1 (MQ=255)
                                                                              gACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCAtt                                                                                                                                  >  6:8384/1‑90 (MQ=255)
                                                                              gACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCAtt                                                                                                                                  <  5:8384/90‑1 (MQ=255)
                                                                                   tGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCg                                                                    <  6:321688/146‑1 (MQ=255)
                                                                                   tGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCg                                                                    >  5:321688/1‑146 (MQ=255)
                                                                                      tACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCa                                                              <  5:87786/149‑1 (MQ=255)
                                                                                          tGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAActt                                                          >  4:8586/1‑149 (MQ=255)
                                                                                              gcACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAActtctt                                                       >  6:169516/1‑148 (MQ=255)
                                                                                              gcACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACttcttc                                                      <  7:161433/149‑1 (MQ=255)
                                                                                               cACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTTATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACTCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCa                                                     <  1:120497/149‑1 (MQ=255)
                                                                                                 ccAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATc                                                   <  3:83274/149‑1 (MQ=255)
                                                                                                  cAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAg                                                                 >  5:39357/1‑134 (MQ=255)
                                                                                                  cAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAg                                                                 <  6:39357/134‑1 (MQ=255)
                                                                                                    gTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAg                                                >  4:60474/1‑149 (MQ=255)
                                                                                                              cTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGTCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCg                                       >  1:411875/1‑148 (MQ=255)
                                                                                                                   gCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACctg                                 <  3:8586/149‑1 (MQ=255)
                                                                                                                   gCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACctg                                 >  3:354500/1‑149 (MQ=255)
                                                                                                                       cGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACctggact                             <  3:60474/149‑1 (MQ=255)
                                                                                                                           tACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGaa                         >  1:177599/1‑149 (MQ=255)
                                                                                                                             caagaacaagaTCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGaa                         <  5:169516/147‑1 (MQ=255)
                                                                                                                                   acaagaTCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAac                  <  4:425269/149‑1 (MQ=255)
                                                                                                                                    caagaTCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAAcc                                   <  7:355234/131‑1 (MQ=255)
                                                                                                                                    caagaTCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAAcc                                   >  8:355234/1‑131 (MQ=255)
                                                                                                                                     aagcacatcGGTATCACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGGGGGACAACGCCGGCAAGGCGCGTGTGGAAGGTGTCGAAGGCAACTTCTTCATAGAGTTGGACCCGACCCTGGACTGGAACACCAacct                >  6:310194/6‑149 (MQ=255)
                                                                                                                                       gcacatcGGTATCACCGATGTGCTTTACACCATTGGTGACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAacct                >  1:296570/4‑148 (MQ=255)
                                                                                                                                        cacatcGGTATCACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAacctgac             >  6:60035/3‑149 (MQ=255)
                                                                                                                                            tcGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAAAGCCGTCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGTACACCAccctgacctg          >  2:257764/1‑148 (MQ=255)
                                                                                                                                              ggTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGa        <  4:354500/148‑1 (MQ=255)
                                                                                                                                              ggTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGAt       <  8:282706/149‑1 (MQ=255)
                                                                                                                                                tGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGc     >  3:293548/1‑149 (MQ=255)
                                                                                                                                                  aTACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGCTc   <  5:152266/149‑1 (MQ=255)
                                                                                                                                                  aTACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGCTc   <  2:205553/149‑1 (MQ=255)
                                                                                                                                                   tACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGCTcg  <  4:293548/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATGGTGGCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGCTCG  >  NC_002947/2559304‑2559598

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