Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,559,448 G→A 9.2% D570N (GAT→AAT) ‡ 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,4480GA 9.2% 182.8 / 6.6 66D570N (GAT→AAT) ‡fepAferric enterobactin transport system outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (3/3);  total (33/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 = 8.88e-01

ACAATGGTGGCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGC  >  NC_002947/2559302‑2559595
                                                                                                                                                   |                                                                                                                                                    
acaaTGGTGGCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGa                                                                                                                                                     >  6:87786/1‑148 (MQ=255)
  aaTGGTGGCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATAc                                                                                                                                                  <  1:436847/149‑1 (MQ=255)
      gtggCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGAt                                                                                                                                              >  1:180231/1‑149 (MQ=255)
       tggCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTATCACCGAtg                                                                                                                                             <  2:437876/149‑1 (MQ=255)
         gCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCCGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTATCACCGAtga                                                                                                                                            <  7:137318/148‑2 (MQ=255)
          cTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGACACAAGATCATCGGTGATACCGAtgtg                                                                                                                                           <  1:3615/149‑1 (MQ=255)
          cTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTATCACCGAtgatc                                                                                                                                          <  6:169058/149‑4 (MQ=255)
               aCCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTa                                                                                                                                      <  4:379369/148‑1 (MQ=255)
                ccTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTac                                                                                                                                     <  2:230226/148‑1 (MQ=255)
                 cTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTacac                                                                                                                                   <  5:432701/149‑1 (MQ=255)
                    cAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCa                                                                                                                                 <  2:316005/148‑1 (MQ=255)
                    cAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCAt                                                                                                                                >  5:210853/1‑149 (MQ=255)
                        gCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTgg                                                                                                                             <  4:70393/148‑1 (MQ=255)
                        gCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTgg                                                                                                                             >  3:70393/1‑148 (MQ=255)
                        gCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGt                                                                                                                            <  1:362576/149‑1 (MQ=255)
                                   cTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCcg                                                                                                                >  5:378407/1‑149 (MQ=255)
                                    tCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCcgc                                                                                                               >  6:263810/1‑149 (MQ=255)
                                      aaGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGt                                                                                                             >  1:288731/1‑149 (MQ=255)
                                         cccGAGATCAGCGTCAACAAGGAAATCGTCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATTTGCTTTACAGCATTGGT‑ACCGGCAGCCGCGTGAc                                                                                                          >  5:93127/1‑149 (MQ=255)
                                          ccGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCCCAACGACTACAATA‑ACAAGATCATCGGTGATACCGATGTCCTTTACACCATTGGT‑ACCGGCAGCCCCGTGACg                                                                                                         >  7:333851/1‑149 (MQ=255)
                                               aTCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTg                                                                                                    <  5:263810/149‑1 (MQ=255)
                                               aTCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACAGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTg                                                                                                    <  7:291376/149‑1 (MQ=255)
                                                tCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTgg                                                                                                   >  1:422238/1‑149 (MQ=255)
                                                 cAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTgg                                                                                                   >  7:156024/1‑148 (MQ=255)
                                                  aGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTg                                                                                                                              >  8:276407/1‑121 (MQ=255)
                                                  aGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTg                                                                                                                              <  7:276407/121‑1 (MQ=255)
                                                     gTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACaa                                                                                              <  2:422238/149‑1 (MQ=255)
                                                           aaGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCgg                                                                                        >  6:75621/1‑149 (MQ=255)
                                                             ggAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGc                                                                                       >  6:264437/1‑148 (MQ=255)
                                                                      ccTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCggg                                                                             >  6:152266/1‑149 (MQ=255)
                                                                      ccTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGAGCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCggg                                                                             <  5:36539/149‑1 (MQ=255)
                                                                        tgctgTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTg                                                                           >  4:83274/1‑149 (MQ=255)
                                                                             tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAgg                                                                      >  1:205553/1‑149 (MQ=255)
                                                                             tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAgg                                                                      >  3:425269/1‑149 (MQ=255)
                                                                             tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAg                                                                       <  6:378407/148‑1 (MQ=255)
                                                                             tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAg                                                                       <  5:264437/148‑1 (MQ=255)
                                                                             tACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAg                                                                       <  8:156024/148‑1 (MQ=255)
                                                                                gACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCAtt                                                                                                                               >  6:8384/1‑90 (MQ=255)
                                                                                gACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCAtt                                                                                                                               <  5:8384/90‑1 (MQ=255)
                                                                                     tGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCg                                                                 >  5:321688/1‑146 (MQ=255)
                                                                                     tGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCg                                                                 <  6:321688/146‑1 (MQ=255)
                                                                                        tACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCa                                                           <  5:87786/149‑1 (MQ=255)
                                                                                            tGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAActt                                                       >  4:8586/1‑149 (MQ=255)
                                                                                                gcACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAActtctt                                                    >  6:169516/1‑148 (MQ=255)
                                                                                                gcACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACttcttc                                                   <  7:161433/149‑1 (MQ=255)
                                                                                                 cACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTTATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACTCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCa                                                  <  1:120497/149‑1 (MQ=255)
                                                                                                   ccAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATc                                                <  3:83274/149‑1 (MQ=255)
                                                                                                    cAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAg                                                              >  5:39357/1‑134 (MQ=255)
                                                                                                    cAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAg                                                              <  6:39357/134‑1 (MQ=255)
                                                                                                      gTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAg                                             >  4:60474/1‑149 (MQ=255)
                                                                                                                cTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGTCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCg                                    >  1:411875/1‑148 (MQ=255)
                                                                                                                     gCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACctg                              <  3:8586/149‑1 (MQ=255)
                                                                                                                     gCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACctg                              >  3:354500/1‑149 (MQ=255)
                                                                                                                         cGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACctggact                          <  3:60474/149‑1 (MQ=255)
                                                                                                                             tACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGaa                      >  1:177599/1‑149 (MQ=255)
                                                                                                                               caaga‑acaagaTCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGaa                      <  5:169516/147‑1 (MQ=255)
                                                                                                                                     acaagaTCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAac               <  4:425269/149‑1 (MQ=255)
                                                                                                                                      caagaTCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAAcc                                <  7:355234/131‑1 (MQ=255)
                                                                                                                                      caagaTCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAAcc                                >  8:355234/1‑131 (MQ=255)
                                                                                                                                       aagcacatcGGTATCACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGGGGGACAACGCCGGCAAGGCGCGTGTGGAAGGTGTCGAAGGCAACTTCTTCATAGAGTTGGACCCGACCCTGGACTGGAACACCAacct             >  6:310194/6‑149 (MQ=255)
                                                                                                                                         gcacatcGGTATCACCGATGTGCTTTACACCATTGGTGACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAacct             >  1:296570/4‑148 (MQ=255)
                                                                                                                                          cacatcGGTATCACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAacctgac          >  6:60035/3‑149 (MQ=255)
                                                                                                                                              tcGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAAAGCCGTCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGTACACCAccctgacctg       >  2:257764/1‑148 (MQ=255)
                                                                                                                                                ggTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGa     <  4:354500/148‑1 (MQ=255)
                                                                                                                                                ggTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGAt    <  8:282706/149‑1 (MQ=255)
                                                                                                                                                  tGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGc  >  3:293548/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACAATGGTGGCTGTTACCTGCAGGGCAACGCCGACCTCAAGCCCGAGATCAGCGTCAACAAGGAAATCGGCCTGCTGTACGACCGTGGTACCTGGCGCACCAGTGCCACCTACTTCCGCAACGACTACAAGA‑ACAAGATCATCGGTGATACCGATGTGCTTTACACCATTGGT‑ACCGGCAGCCGCGTGACGCAGTGGGACAACGCCGGCAAGGCGCGGGTGGAAGGTGTCGAAGGCAACTTCTTCATCGAGTTGGCCCCGAACCTGGACTGGAACACCAACCTGACCTGGATGC  >  NC_002947/2559302‑2559595

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: ‑