Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,559,449 A→T 9.2% D570V (GAT→GTT) ‡ 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,4490AT 9.2% 185.3 / 7.9 65D570V (GAT→GTT) ‡fepAferric enterobactin transport system outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base T (3/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

ACAATGGTGGCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGCTC  >  NC_002947/2559302‑2559597
                                                                                                                                                    |                                                                                                                                                     
acaaTGGTGGCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGa                                                                                                                                                       >  6:87786/1‑148 (MQ=255)
  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                                                                                                                               <  4:70393/148‑1 (MQ=255)
                        gCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTgg                                                                                                                               >  3:70393/1‑148 (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                                                                                                                                >  8:276407/1‑121 (MQ=255)
                                                  aGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTg                                                                                                                                <  7:276407/121‑1 (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                                                                         <  6:378407/148‑1 (MQ=255)
                                                                             tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAg                                                                         <  5:264437/148‑1 (MQ=255)
                                                                             tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAg                                                                         <  8:156024/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)
                                                                                                                                                    |                                                                                                                                                     
ACAATGGTGGCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGAACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGTACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGCTC  >  NC_002947/2559302‑2559597

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: