Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,283,203 G→T 5.2% L249M (CTG→ATG)  yejK ← spermidine nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,283,2030GT 5.2% 260.1 / 2.2 76L249M (CTG→ATG) yejKspermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base T (2/2);  total (39/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 = 9.11e-01

GCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCGCTTGCAACAGCC  >  NC_000913/2283055‑2283349
                                                                                                                                                    |                                                                                                                                                  
gcatcAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTTACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAg                                                                                                                                                    >  5:332624/1‑149 (MQ=255)
    cAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGctct                                                                                                                                                >  5:262053/1‑149 (MQ=255)
        gTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTACAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc                                                                                                                                             >  3:358913/1‑148 (MQ=255)
          tAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGATGCTCACTTCGCTTACGCCTGCCAGCTCTGTCGa                                                                                                                                           >  1:204095/1‑148 (MQ=255)
           aaTAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGa                                                                                                                                           >  3:83541/1‑147 (MQ=255)
               gTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGa                                                                                                                                     >  4:169823/1‑149 (MQ=255)
                  aaGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGAttt                                                                                                                                  >  7:45135/1‑149 (MQ=255)
                  aaGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTCACGCCTGCCAGCTCTTTCGACAGCGAttt                                                                                                                                  <  8:169150/149‑1 (MQ=255)
                       aCCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGtt                                                                                                                             >  6:158541/1‑149 (MQ=255)
                       aCCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGtt                                                                                                                             >  4:150595/1‑149 (MQ=255)
                       aCCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGtt                                                                                                                             <  2:300920/149‑1 (MQ=255)
                        ccgccgCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTc                                                                                                                            <  8:338730/149‑1 (MQ=255)
                        ccgccgCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTc                                                                                                                            >  6:338754/1‑149 (MQ=255)
                          gccgcTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGa                                                                                                                          <  5:381608/149‑1 (MQ=255)
                          gccgcTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGa                                                                                                                          <  6:262053/149‑1 (MQ=255)
                           ccgcTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGAt                                                                                                                         <  8:322834/149‑1 (MQ=255)
                           ccgcTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTATTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGATCTTTCGACAGCGATTTCAGTTCGa                                                                                                                          >  6:200125/1‑148 (MQ=255)
                               tACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCTGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATctct                                                                                                                     >  3:260294/1‑149 (MQ=255)
                                 ccGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTc                                                                                                                   <  5:222458/149‑1 (MQ=255)
                                              tCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAgc                                                                                                      <  5:338754/149‑1 (MQ=255)
                                                 aCTGCCGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAgctgc                                                                                                   <  4:358913/149‑1 (MQ=255)
                                                       gCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCATCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTa                                                                                             >  2:21767/1‑149 (MQ=255)
                                                       gCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCATCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTa                                                                                             >  2:15352/1‑149 (MQ=255)
                                                       gCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTa                                                                                             <  8:21071/149‑1 (MQ=255)
                                                           tgtgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGAtctc                                                                                                                      <  2:376722/120‑1 (MQ=255)
                                                           tgtgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGAtctc                                                                                                                      >  1:376722/1‑120 (MQ=255)
                                                               cTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGCCAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGc                                                                                      >  7:123522/1‑148 (MQ=255)
                                                               cTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGc                                                                                      <  3:169823/148‑1 (MQ=255)
                                                               cTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGc                                                                                      >  4:69701/1‑148 (MQ=255)
                                                                tACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGAGATTTCACTTCGATCTCTTCGACTGCCTGCAGCTGATCATTACAGTAGCta                                                                                    >  5:282271/1‑149 (MQ=255)
                                                                 aCGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCta                                                                                    >  3:247036/1‑148 (MQ=255)
                                                                    gTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAg                                                                                       <  7:343363/142‑1 (MQ=255)
                                                                    gTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAg                                                                                       >  8:343363/1‑142 (MQ=255)
                                                                             gAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGgc                                                                       <  2:190936/149‑1 (MQ=255)
                                                                                      ttCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCAc                                                                   <  1:308241/144‑1 (MQ=255)
                                                                                      ttCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCAc                                                                   >  2:308241/1‑144 (MQ=255)
                                                                                      ttCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGtt                                                                <  4:247036/147‑1 (MQ=255)
                                                                                      ttCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCt                                                              >  4:138580/1‑149 (MQ=255)
                                                                                        ccAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGt                                                            >  4:369809/1‑149 (MQ=255)
                                                                                         cAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGt                                                            <  3:150595/148‑1 (MQ=255)
                                                                                           gTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGc                                                         >  8:37235/1‑149 (MQ=255)
                                                                                           gTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTGCTGTCGc                                                         <  2:199002/149‑1 (MQ=255)
                                                                                               gTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTg                                                     <  8:154304/149‑1 (MQ=255)
                                                                                               gTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTg                                                     >  3:324976/1‑149 (MQ=255)
                                                                                               gTAGCCTTTTTCAGCCGCAAACTCCGTGAAGATCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCATGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGTCGCAAGTTCTGTCGCTCTg                                                     >  8:371420/1‑149 (MQ=255)
                                                                                                tAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCTCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGc                                                    >  5:345758/1‑149 (MQ=255)
                                                                                                tAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGc                                                    <  7:219986/149‑1 (MQ=255)
                                                                                                tAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGc                                                    >  4:42357/1‑149 (MQ=255)
                                                                                                    cTTTATCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTa                                                <  1:404453/149‑1 (MQ=255)
                                                                                                         tCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACTCCTTCCAGCTCTTTCGAAAGCGATTTCAGTTCGAGCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCGGGCGCACGTTCTGTCGATCTGCTTTATCCa                                            <  6:363930/148‑1 (MQ=255)
                                                                                                          cAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAg                                           <  3:42357/148‑1 (MQ=255)
                                                                                                             ccGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAg                                                                                                       >  5:196509/1‑85 (MQ=255)
                                                                                                             ccGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTATTCGCCTGCCTGCAg                                                                                                       <  6:196509/85‑1 (MQ=255)
                                                                                                              cGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTtgtg                                      >  7:26543/1‑149 (MQ=255)
                                                                                                                  aaCTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTc                                  <  7:31546/149‑1 (MQ=255)
                                                                                                                  aaCTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTc                                  >  2:17320/1‑149 (MQ=255)
                                                                                                                   aCTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAACCCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCg                                 >  2:332738/1‑149 (MQ=255)
                                                                                                                    cTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCgg                                >  6:318109/1‑149 (MQ=255)
                                                                                                                    cTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCgg                                <  5:200125/149‑1 (MQ=255)
                                                                                                                     tCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCtata                                                                                  <  8:112342/98‑1 (MQ=255)
                                                                                                                     tCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCtata                                                                                  >  7:112342/1‑98 (MQ=255)
                                                                                                                       cGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGc                               <  1:262745/147‑1 (MQ=255)
                                                                                                                                 aCTTCGCTTACGCCTGCCATCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAatcatc                   <  1:15352/149‑1 (MQ=255)
                                                                                                                                 aCTTCGCTTACGCCTGCCATCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAatcatc                   <  1:21767/149‑1 (MQ=255)
                                                                                                                                 aCTTCGCGTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAatcat                    <  5:10304/148‑1 (MQ=255)
                                                                                                                                   ttCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCaa                 >  1:299084/1‑149 (MQ=255)
                                                                                                                                   ttCGCTTACACCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCaa                 >  1:221709/1‑149 (MQ=255)
                                                                                                                                   atCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGGTCATTACAGAAGCTATAAACCTGCTGGCGCACGGTCTGTCGATCTGCTTTATCAAGTTGTGCTTCGTCAGTGAATTCATCaa                 >  8:301742/2‑149 (MQ=255)
                                                                                                                                     cGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCaa                 <  7:37235/147‑1 (MQ=255)
                                                                                                                                     cGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAAc                <  1:332738/148‑1 (MQ=255)
                                                                                                                                        ttACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAAcc               <  4:88813/146‑1 (MQ=255)
                                                                                                                                        ttACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAAcc               >  3:88813/1‑146 (MQ=255)
                                                                                                                                        ttACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCGCt            >  4:136519/1‑149 (MQ=255)
                                                                                                                                          aCGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCGCtt           <  1:294251/148‑1 (MQ=255)
                                                                                                                                          aCGCCTGCCAGCTCTTGCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCGCtt           <  1:288858/148‑1 (MQ=255)
                                                                                                                                              cTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCGCTTGCAAc      <  4:194439/149‑1 (MQ=255)
                                                                                                                                                  cAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCGCTTGCAACAGcc  >  8:68923/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCGCCTGCCTGCAGCTGCTCATTACAGTAGCTATAAACCTGCTGGCGCACGTTCTGTCGCTCTGCTTTATCCAGTTGTGCTTCGGCAGTGAAATCATCAACCGCTTGCAACAGCC  >  NC_000913/2283055‑2283349

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: