Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,235,070 A→G 5.5% H135R (CAC→CGC)  fecA → outer membrane ferric citrate porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,235,0700AG 5.5% 216.4 / 2.2 72H135R (CAC→CGC) fecAouter membrane ferric citrate porin
Reads supporting (aligned to +/- strand):  ref base A (32/36);  new base G (2/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCA  >  NC_002947/5234925‑5235216
                                                                                                                                                 |                                                                                                                                                  
ggCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCAccc                                                                                                                                                 <  5:294385/149‑1 (MQ=255)
    tGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAACAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGgc                                                                                                                                             >  3:49478/1‑149 (MQ=255)
      cAAGCACTCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGgcg                                                                                                                                            <  5:259619/148‑1 (MQ=255)
      cAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGgc                                                                                                                                             <  5:199138/147‑1 (MQ=255)
             gCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGAc                                                                                                                                    >  1:276305/1‑149 (MQ=255)
               cAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGt                                                                                                                                  <  2:269270/149‑1 (MQ=255)
                     aCCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTgc                                                                                                                             <  8:382853/148‑1 (MQ=255)
                     aCCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTgc                                                                                                                             >  7:382853/1‑148 (MQ=255)
                       cGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGcc                                                                                                                          >  2:397815/1‑149 (MQ=255)
                           aaCGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGa                                                                                                                      <  2:309191/149‑1 (MQ=255)
                                gTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAAt                                                                                                                 >  7:216399/1‑149 (MQ=255)
                                      ccTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAAc                                                                                                           >  4:385919/1‑149 (MQ=255)
                                       cTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACg                                                                                                          <  5:276661/149‑1 (MQ=255)
                                        tGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGaa                                                                                                            >  4:22082/1‑146 (MQ=255)
                                        tGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGaa                                                                                                            <  3:22082/146‑1 (MQ=255)
                                        tGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTACACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGAAAACGTATTCGAGCACCCCGGCGCCAGCGACGTAGTGCGCCGGGAGGAATTCGAACg                                                                                                          >  2:310777/1‑148 (MQ=255)
                                         gCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGc                                                                                                         >  3:145305/1‑148 (MQ=255)
                                          cAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGc                                                                                                         <  2:276305/147‑1 (MQ=255)
                                          cAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCa                                                                                                        <  2:72904/148‑1 (MQ=255)
                                             cccATGACAGCCGCTGCTGCTGGCGCATCAGTCGTGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAg                                                                                                                     <  4:37056/132‑1 (MQ=255)
                                             cccATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGTGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACgg                                                                                                    <  2:119419/149‑1 (MQ=255)
                                             cccATCACAGCCCCTGCTGCTGGCGAAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGCCGCCCGCGACGTAGTGCGCCGTGAg                                                                                                                     >  3:37056/1‑132 (MQ=255)
                                                aTCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGc                                                                                                 <  3:385919/149‑1 (MQ=255)
                                                 tCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCt                                                                                                <  1:397815/149‑1 (MQ=255)
                                                                ctgGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTacca                                                                                            >  5:367385/1‑138 (MQ=255)
                                                                ctgGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTacca                                                                                            <  6:367385/138‑1 (MQ=255)
                                                                   gcgcAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTgcg                                                                                                                            >  6:197645/1‑103 (MQ=255)
                                                                   gcgcAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTgcg                                                                                                                            <  5:197645/103‑1 (MQ=255)
                                                                   gcgcAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAg                                                                                <  3:163074/147‑1 (MQ=255)
                                                                    cgcAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGc                                                                             <  2:402158/149‑1 (MQ=255)
                                                                     gcagcaGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCt                                                                            >  5:172272/1‑149 (MQ=255)
                                                                     gcagcaGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCATGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCt                                                                            <  3:58611/149‑1 (MQ=255)
                                                                       agcagTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTc                                                                           >  6:129882/1‑148 (MQ=255)
                                                                          agTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAAcc                                                                       <  5:271773/149‑1 (MQ=255)
                                                                            tCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCGCCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCg                                                                      >  2:265991/1‑148 (MQ=255)
                                                                            tCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCGCCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCg                                                                      >  2:271672/1‑148 (MQ=255)
                                                                              gAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGGAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCAt                                                                   >  1:284647/1‑149 (MQ=255)
                                                                                   aGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCg                                                              >  7:124077/1‑149 (MQ=255)
                                                                                     gcgcCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCggg                                                            <  1:321420/149‑1 (MQ=255)
                                                                                     gcgcCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCggg                                                            <  5:129882/149‑1 (MQ=255)
                                                                                        ccTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTc                                                         <  4:49478/149‑1 (MQ=255)
                                                                                            cACCGTGGTCGGCGACTGGCTGGCCGAGGCGCCGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTCCCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACg                                                     >  5:986/1‑149 (MQ=255)
                                                                                            cACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCGCCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAAc                                                      <  1:271672/148‑1 (MQ=255)
                                                                                            cACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCGCCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAAc                                                      <  1:265991/148‑1 (MQ=255)
                                                                                            cACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAAc                                                      <  3:189554/148‑1 (MQ=255)
                                                                                            cACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACg                                                     >  5:98934/1‑149 (MQ=255)
                                                                                                 tGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGccccc                                                <  8:293394/149‑1 (MQ=255)
                                                                                                   gTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCg                                               <  5:360841/148‑1 (MQ=255)
                                                                                                   gTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGa                                              <  8:216399/149‑1 (MQ=255)
                                                                                                      ggCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGaa                                                                                 <  7:365353/111‑1 (MQ=255)
                                                                                                      ggCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGaa                                                                                 >  8:365353/1‑111 (MQ=255)
                                                                                                      ggCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCgggg                                                           >  8:223404/1‑133 (MQ=255)
                                                                                                      ggCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCgggg                                                           <  7:223404/133‑1 (MQ=255)
                                                                                                           ctggctggCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACg                                      <  4:145305/149‑1 (MQ=255)
                                                                                                            tggctggcCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACgg                                     >  3:281354/1‑149 (MQ=255)
                                                                                                              gctggcCGAGGCGCAGCAGGACAACGTATTCGACCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTa                                   >  6:8068/1‑149 (MQ=255)
                                                                                                                 ggcCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCg                                >  5:115829/1‑149 (MQ=255)
                                                                                                                  gcCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCg                                <  6:147182/148‑1 (MQ=255)
                                                                                                                    cGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCTACGTAGTGCGCCGTGAGGAATTCTAACGCAACGGTGCTACCACCGCCCGCTATGTGCTCAACCGCATCCCCGTGGTCAACGCCCCCGACAACAACGGTACCTGCa                             <  1:310777/149‑1 (MQ=255)
                                                                                                                     gAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAg                            >  6:353016/1‑149 (MQ=255)
                                                                                                                          gcagcaGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGt                                    >  7:358519/1‑136 (MQ=255)
                                                                                                                          gcagcaGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGt                                    <  8:358519/136‑1 (MQ=255)
                                                                                                                            agcagGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAg                                                                                                                     >  4:81667/1‑53 (MQ=255)
                                                                                                                                ggACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACACCCACGACAACAACGGTACCGGCAGCCATGACCTgg                 >  5:386445/1‑149 (MQ=255)
                                                                                                                                 gACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCACCCGACAACAACGGTACCGGCAGCCATGACCTGgc                >  5:186217/1‑149 (MQ=255)
                                                                                                                                        tATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAAc         <  6:186217/149‑1 (MQ=255)
                                                                                                                                         aTTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTg                                                                              >  7:121049/1‑79 (MQ=255)
                                                                                                                                         aTTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTg                                                                              <  8:121049/79‑1 (MQ=255)
                                                                                                                                         aTTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGCGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACt        >  4:57654/1‑149 (MQ=255)
                                                                                                                                            cGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTc      <  8:124077/148‑1 (MQ=255)
                                                                                                                                             gAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCgg    >  3:343159/1‑149 (MQ=255)
                                                                                                                                               gCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCa  >  5:2062/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCCTGCAAGCACGCCAGGAAACCGACAACGCGTTCAGCCTGCAACCCATCACAGCCCCTGCTGCTGGCGCAGCAGTCGAGCTAGGCGCCTCCACCGTGGTCGGCGACTGGCTGGCCGAGGCGCAGCAGGACAACGTATTCGAGCACCCCGGCGCCCGCGACGTAGTGCGCCGTGAGGAATTCGAACGCAACGGTGCTACCACCGCCCGCGAAGTGCTCAACCGCATCCCCGGGGTCAACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCA  >  NC_002947/5234925‑5235216

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