Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,054,758 A→G 9.7% K46E (AAG→GAG) ‡ 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,7580AG 9.7% 224.2 / 12.8 82K46E (AAG→GAG) ‡exbDTonB‑gated outer membrane transporter gating inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (40/34);  new base G (4/4);  total (44/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 = 4.36e-01

TGGGGTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAA  >  NC_002947/6054612‑6054903
                                                                                                                                                    |                                                                                                                                                 
tGGGGTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAATACCTCGAAATCAACGTTACGCCTGTCATCGACGTGTGGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCGTTGGCTACGGTCGACATCAAg                                                                                                                                                 <  2:143450/149‑1 (MQ=255)
  gggTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAgg                                                                                                                                                <  3:299868/148‑1 (MQ=255)
  gggTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCACTTGGCTACGGTCGACATCacagt                                                                                                                                               >  4:194603/1‑145 (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                                                                                                                                       >  8:154238/1‑132 (MQ=39)
                        tcAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACCTTGATGTCGAccgta                                                                                                                                       <  7:154238/135‑4 (MQ=39)
                         cAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTCCCGGCCTCAACCGcc                                                                                                                         >  8:305413/1‑148 (MQ=255)
                          aaCGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCa                                                                                                                        >  8:43915/1‑148 (MQ=255)
                          aaCGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCa                                                                                                                        >  8:49737/1‑148 (MQ=255)
                            cGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAAc                                                                                                                     >  8:141425/1‑149 (MQ=255)
                            cGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAAc                                                                                                                     >  7:59647/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                                                                                                                                       >  6:75178/1‑117 (MQ=38)
                                       acgacCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACCTTGATGTCGAccgta                                                                                                                                       <  5:75178/120‑4 (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                                                                                                                                             >  8:33117/1‑85 (MQ=255)
                                                                    gTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCg                                                                                                                                             <  7:33117/85‑1 (MQ=255)
                                                                     ttACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAg                                                                            <  4:167790/149‑1 (MQ=255)
                                                                     ttACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAg                                                                            <  7:49737/149‑1 (MQ=255)
                                                                     ttACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAg                                                                            <  7:43915/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                             <  8:304417/149‑1 (MQ=255)
                                                                                                                    gTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGa                             >  5:92448/1‑149 (MQ=255)
                                                                                                                       gccgcTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGAc                            >  8:189590/1‑147 (MQ=255)
                                                                                                                              ccTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCTGAGAAACCGGTGTTCGTCAGCGTCAAGGCTGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTgg                    <  8:297666/148‑1 (MQ=255)
                                                                                                                                  tGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCtt                      <  2:337332/143‑1 (MQ=255)
                                                                                                                                  tGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCtt                      >  1:337332/1‑143 (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          <  4:314780/139‑1 (MQ=255)
                                                                                                                                                   cAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTc          >  3:314780/1‑139 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGGGGTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAA  >  NC_002947/6054612‑6054903

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: