Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,054,760 G→T 10.1% K46N (AAG→AAT) ‡ exbD → TonB‑gated outer membrane transporter gating inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,054,7600GT10.1% 224.2 / 14.8 79K46N (AAG→AAT) ‡exbDTonB‑gated outer membrane transporter gating inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (40/31);  new base T (4/4);  total (44/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 = 3.04e-02

TGGGGTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGAC  >  NC_002947/6054612‑6054906
                                                                                                                                                      |                                                                                                                                                  
tGGGGTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAATACCTCGAAATCAACGTTACGCCTGTCATCGACGTGTGGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCGTTGGCTACGGTCGACATCAAg                                                                                                                                                    <  2:143450/149‑1 (MQ=255)
  gggTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAgg                                                                                                                                                   <  3:299868/148‑1 (MQ=255)
    gTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTc                                                                                                                                                 >  8:305388/1‑148 (MQ=255)
    gTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCAGCTCCCTTGGCTACGGTCGACATCAAGGTc                                                                                                                                                 >  1:156097/1‑148 (MQ=255)
     tAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGa                                                                                                                                               <  5:174958/149‑1 (MQ=255)
      aaGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGAc                                                                                                                                              >  3:167790/1‑149 (MQ=255)
        gCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACCTTGATGTCGAccg                                                                                                                                            <  5:375512/149‑2 (MQ=255)
        gCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCt                                                                                                                                            >  8:167133/1‑149 (MQ=255)
          cATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGc                                                                                                                                          >  6:94682/1‑149 (MQ=255)
           aTGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACCTTGATGTCGAccgtag                                                                                                                                         <  4:146636/149‑5 (MQ=39)
           aTGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGcc                                                                                                                                         >  2:14814/1‑149 (MQ=255)
               gCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGcc                                                                                                                                     >  5:91342/1‑149 (MQ=255)
                ccTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACCGGTCGACATCAAGGTCGACCTGCCGGcc                                                                                                                                     <  1:210080/149‑1 (MQ=255)
                ccTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCt                                                                                                                                    <  7:305388/149‑1 (MQ=255)
                  tGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTc                                                                                                                                   >  2:84275/1‑148 (MQ=255)
                        tcAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACCTTGATGTCGAccgta                                                                                                                                          <  7:154238/135‑4 (MQ=39)
                        tcAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACCTTGATGTCGAccgta                                                                                                                                          >  8:154238/1‑132 (MQ=39)
                         cAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTCCCGGCCTCAACCGcc                                                                                                                            >  8:305413/1‑148 (MQ=255)
                          aaCGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCa                                                                                                                           >  8:49737/1‑148 (MQ=255)
                          aaCGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCa                                                                                                                           >  8:43915/1‑148 (MQ=255)
                            cGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAAc                                                                                                                        >  7:59647/1‑149 (MQ=255)
                            cGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAAc                                                                                                                        >  8:141425/1‑149 (MQ=255)
                                  tGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAAcc                                                                                                                               >  3:113006/1‑136 (MQ=255)
                                  tGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAAcc                                                                                                                               <  4:113006/136‑1 (MQ=255)
                                   ggCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGgcgc                                                                                                                  >  7:304814/1‑148 (MQ=255)
                                     cgacgaCCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCaa                                                                                                               <  1:14814/149‑1 (MQ=255)
                                       acgacCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACCTTGATGTCGAccgta                                                                                                                                          <  5:75178/120‑4 (MQ=38)
                                       acgacCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACCTTGATGTCGAccgta                                                                                                                                          >  6:75178/1‑117 (MQ=38)
                                            cTCGGCGAATACCACGAAATCAACGTTACGACGTTCATCGAGGTGATGCTGGTGCTTCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGga                                                                                                        <  2:156097/149‑1 (MQ=255)
                                                 cGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAAc                                                                                                   >  5:59937/1‑149 (MQ=255)
                                                   aaaaCCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCg                                                                                                 <  6:91342/149‑1 (MQ=255)
                                                       ccACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGgtgt                                                                                             <  1:84275/149‑1 (MQ=255)
                                                         aCGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTc                                                                                           >  8:330832/1‑149 (MQ=255)
                                                         aCGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAAACGGTGTTc                                                                                           >  2:16525/1‑149 (MQ=255)
                                                             aaTCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCa                                                                                       <  8:59647/149‑1 (MQ=255)
                                                                 aaCGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGt                                                                                   >  4:339475/1‑149 (MQ=255)
                                                                 aaCGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGt                                                                                   <  3:339475/149‑1 (MQ=255)
                                                                  aCGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTc                                                                                  <  7:167133/149‑1 (MQ=255)
                                                                    gTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCg                                                                                                                                                <  7:33117/85‑1 (MQ=255)
                                                                    gTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCg                                                                                                                                                >  8:33117/1‑85 (MQ=255)
                                                                     ttACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAg                                                                               <  7:43915/149‑1 (MQ=255)
                                                                     ttACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAg                                                                               <  7:49737/149‑1 (MQ=255)
                                                                     ttACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAg                                                                               <  4:167790/149‑1 (MQ=255)
                                                                               cATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAg                                                                      <  5:94682/148‑1 (MQ=255)
                                                                                aTCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAg                                                                      >  5:235327/1‑147 (MQ=255)
                                                                                aTCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGaa                                                                    <  8:304814/149‑1 (MQ=255)
                                                                                 tCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAg                                                                   <  1:217114/149‑1 (MQ=255)
                                                                                    aCGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCTGACCAGAAGCTg                                                                >  7:297666/1‑149 (MQ=255)
                                                                                    aCGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTg                                                                <  7:141425/149‑1 (MQ=255)
                                                                                       tGATGCTGGTTCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTCCGTCAGCGTCAAGGCCGACCAGAAGCTGTAt                                                             <  3:194603/149‑1 (MQ=255)
                                                                                           gCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCg                                                         >  7:304417/1‑149 (MQ=255)
                                                                                            cTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCg                                                         <  7:305413/148‑1 (MQ=255)
                                                                                                   tgctgATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAg                                                                                                                                                    >  1:215098/1‑50 (MQ=255)
                                                                                                   tgctgATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAg                                                                                                                                                    <  2:215098/50‑1 (MQ=255)
                                                                                                     ctgATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGAt                                                 >  5:67122/1‑147 (MQ=255)
                                                                                                     ctgATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGAAATCAAGTTCTACATGCATTCATAAACATACAAAAAGGCTCAAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGATGTATGTCGGCGACGAt                                                 >  5:49809/1‑147 (MQ=255)
                                                                                                      tgATCATCTTCTTGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGt                                                                                         <  2:3662/106‑1 (MQ=255)
                                                                                                      tgATCATCTTCATGGTCGCCGCTCCCGTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGt                                                                                         >  1:3662/1‑107 (MQ=255)
                                                                                                        atcatcTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGt                                            >  1:200663/1‑149 (MQ=255)
                                                                                                        atcatcTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTAGACCTGCCTGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGt                                            >  4:67097/1‑149 (MQ=255)
                                                                                                            tcTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGATCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTtgct                                        >  8:152770/1‑149 (MQ=255)
                                                                                                            tcTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTtgct                                        >  3:179322/1‑149 (MQ=255)
                                                                                                                cATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTtgct                                        <  6:59937/145‑1 (MQ=255)
                                                                                                                 aTGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGCATGTCGGCGACGATCAGGTTGCTGCAc                                    <  1:16525/148‑1 (MQ=255)
                                                                                                                  tGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCAccc                                  >  4:90002/1‑149 (MQ=255)
                                                                                                                    gTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGa                                <  5:49314/149‑1 (MQ=255)
                                                                                                                    gTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGa                                >  5:92448/1‑149 (MQ=255)
                                                                                                                    gTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGa                                <  8:304417/149‑1 (MQ=255)
                                                                                                                       gccgcTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGAc                               >  8:189590/1‑147 (MQ=255)
                                                                                                                              ccTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCTGAGAAACCGGTGTTCGTCAGCGTCAAGGCTGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTgg                       <  8:297666/148‑1 (MQ=255)
                                                                                                                                  tGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCtt                         >  1:337332/1‑143 (MQ=255)
                                                                                                                                  tGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCtt                         <  2:337332/143‑1 (MQ=255)
                                                                                                                                    gCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGAt                 <  3:67097/149‑1 (MQ=255)
                                                                                                                                      tACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGc               >  6:52410/1‑149 (MQ=255)
                                                                                                                                       gcaggTCGACCTTGATGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGc               >  7:275225/4‑149 (MQ=255)
                                                                                                                                        caggTCGACCTTGATGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCt              >  5:160173/3‑149 (MQ=255)
                                                                                                                                            tCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCGGATGCTCGAc          >  3:55767/1‑149 (MQ=255)
                                                                                                                                               aCATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGcc       <  3:128029/149‑1 (MQ=255)
                                                                                                                                                 aTCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGcc       >  2:311486/1‑147 (MQ=255)
                                                                                                                                                 aTCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCaa     >  3:277416/1‑149 (MQ=255)
                                                                                                                                                   cAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTc             >  3:314780/1‑139 (MQ=255)
                                                                                                                                                   cAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTc             <  4:314780/139‑1 (MQ=255)
                                                                                                                                                     aGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCTGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGAc  >  5:222619/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TGGGGTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGAC  >  NC_002947/6054612‑6054906

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: