Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,054,757 C→T 9.9% I45I (ATC→ATT) ‡ 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,7570CT 9.9% 223.2 / 11.5 81I45I (ATC→ATT) ‡exbDTonB‑gated outer membrane transporter gating inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (40/33);  new base T (4/4);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

AAGTGGGGTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAA  >  NC_002947/6054609‑6054903
                                                                                                                                                      |                                                                                                                                                  
aaGTGGGGTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATc                                                                                                                                                    <  5:66231/149‑1 (MQ=255)
   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                                                                                                                     >  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                                                                                                                                       >  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                                                                                <  3:339475/149‑1 (MQ=255)
                                                                    aaCGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGt                                                                                >  4:339475/1‑149 (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                                                                            <  4:167790/149‑1 (MQ=255)
                                                                        ttACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAg                                                                            <  7:49737/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                                                                                                                                                 <  2:215098/50‑1 (MQ=255)
                                                                                                      tgctgATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAg                                                                                                                                                 >  1:215098/1‑50 (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                             <  8:304417/149‑1 (MQ=255)
                                                                                                                       gTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGa                             <  5:49314/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          >  3:314780/1‑139 (MQ=255)
                                                                                                                                                      cAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTc          <  4:314780/139‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
AAGTGGGGTAAGCCATGGGCCTGCATCTCAACGAAGGTGGCGACGACCTCGCCGAAAACCACGAAATCAACGTTACGCCGTTCATCGACGTGATGCTGGTGCTGCTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAA  >  NC_002947/6054609‑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: