Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 751,156 G→A 5.2% R338R (CGC→CGT ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913751,1560GA 5.2% 246.9 / 2.9 77R338R (CGC→CGTybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (43/30);  new base A (2/2);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

ATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCGGGACTTTTTTCTGGT  >  NC_000913/751009‑751301
                                                                                                                                                   |                                                                                                                                                 
atcaGAAATCAGCAGGCCGCCGTACAGCGAGGTGGTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGc                                                                                                                                                  <  4:316239/149‑1 (MQ=255)
 tcaGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGc                                                                                                                                                  >  2:301616/1‑148 (MQ=255)
 tcaGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGc                                                                                                                                                  >  1:147373/1‑148 (MQ=255)
     aaaTCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTc                                                                                                                                             <  5:56363/149‑1 (MQ=255)
          agcagGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTCTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGaa                                                                                                                                        <  8:287864/149‑1 (MQ=255)
           gcagGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAg                                                                                                                                       >  8:316726/1‑149 (MQ=255)
           gcagGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGAGAGTCAGGAAg                                                                                                                                       >  8:316721/1‑149 (MQ=255)
           gcagGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGATTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAg                                                                                                                                       >  1:184682/1‑149 (MQ=255)
            cagGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAgg                                                                                                                                      <  4:310307/149‑1 (MQ=255)
                 cgccgTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGaa                                                                                                                                        <  2:350693/142‑1 (MQ=255)
                 cgccgTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGaa                                                                                                                                        >  1:350693/1‑142 (MQ=255)
                     gTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAtg                                                                                                                             <  6:25061/149‑1 (MQ=255)
                      tACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGCCCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAtgc                                                                                                                            >  1:109060/1‑149 (MQ=255)
                        cAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGCCCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAtgctg                                                                                                                          >  5:272928/1‑149 (MQ=255)
                                  tttGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTg                                                                                                                >  8:84294/1‑149 (MQ=255)
                                      agagCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGaaa                                                                                                            <  1:301616/149‑1 (MQ=255)
                                           aTCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAGGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGtt                                                                                                       >  2:120273/1‑149 (MQ=255)
                                           aTCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGtt                                                                                                       <  7:316726/149‑1 (MQ=255)
                                           aTCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGAGACATGGAAGGCCTTGGCTGACCCGCGTCCAATTTATAGCGAACCTGTCACTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAattatt                                                                                                       >  6:351947/1‑149 (MQ=255)
                                            tCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGAGATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGGTGCCGAAACCTGGAAATTGtt                                                                                                       >  5:142942/1‑148 (MQ=255)
                                             cccccAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTGAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGtt                                                                                                       <  7:316721/147‑1 (MQ=255)
                                              ccccAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAGGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCAc                                                                                                    <  1:120273/149‑1 (MQ=255)
                                               cccAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAAt                                                                                                           <  2:390128/141‑1 (MQ=255)
                                               cccAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAAt                                                                                                           >  1:390128/1‑141 (MQ=255)
                                               cccAGGAAACTTCATTGCTAAAAAAGGTTTCATTTACAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACATGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTCGAAATTCTTcccc                                                                                                   >  7:157363/1‑149 (MQ=255)
                                                       aCTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCAtc                                                                                           >  1:203872/1‑149 (MQ=255)
                                                        cTTCATTGCTAAAAAAGGGTTCATTTTCAGCTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCAtct                                                                                          <  6:132142/149‑1 (MQ=255)
                                                         ttCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTc                                                                                                      >  1:128785/1‑136 (MQ=255)
                                                         ttCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTc                                                                                                      <  2:128785/136‑1 (MQ=255)
                                                                cTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGCCCCGCGGCCAATTTATAGCGAAACTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATGTTCTTCCgg                                                                                  >  3:82494/1‑148 (MQ=255)
                                                                cTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGt                                                                                  >  8:204722/1‑149 (MQ=255)
                                                                 tAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCTCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCAGGAAATTGTTCACCCGCCCATCTTCTTCCGt                                                                                  >  1:54576/1‑148 (MQ=255)
                                                                 tAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGt                                                                                  >  3:262585/1‑148 (MQ=255)
                                                                   aaaaaGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCAc                                                                               <  5:353491/149‑1 (MQ=255)
                                                                   aaaaaGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCAc                                                                               >  5:216161/1‑149 (MQ=255)
                                                                        ggTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGa                                                                          <  6:119671/149‑1 (MQ=255)
                                                                         gTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGa                                                                          >  3:172377/1‑148 (MQ=255)
                                                                            tCATTTTCAGTTTGATGTGACATGGAAGGCCTTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTg                                                                           <  2:322628/144‑1 (MQ=255)
                                                                            tCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCACTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGCCCATCTTCTTCCGTCACTTTg                                                                           >  1:322628/1‑144 (MQ=255)
                                                                               ttttCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAg                                                                   <  7:204722/149‑1 (MQ=255)
                                                                                tttCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGc                                                                  >  7:120887/1‑149 (MQ=255)
                                                                                         gATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTg                                                          <  2:109060/148‑1 (MQ=255)
                                                                                         gATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTg                                                          >  3:118920/1‑148 (MQ=255)
                                                                                          aTGTGACATGGAAGGCCGTGGCTGACCCTCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGaa                                                        >  5:93607/1‑149 (MQ=255)
                                                                                             tgACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACa                                                      <  4:172377/148‑1 (MQ=255)
                                                                                              gACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGAGAACCTGTCCCTGGAGAGTCAGGAAGGGCGGCGATGCTGCCGAAACCTGGAAATTGTTCACCGGACCATCTTCTTCCGGCACTGTGACTTCCAGCGTGCCCTGaagctc                                                    >  5:20085/1‑145 (MQ=255)
                                                                                                   ggAAGGCCGTGGCTGCCCCGCGGCCAATTTATAGCGCACCTGTCCCTGGCGCGGCAGGAAGGGCGTCGCTGCTGCCCAAACCTGGAAATTGTGCACCCGACCATCTTCTTCCGTCACGTTGCCCTCCAGCGTGCCCTGACCAGCCTGag                                               >  4:128691/1‑148 (MQ=255)
                                                                                                   ggAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCCGCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACGGAt                                               >  4:123265/1‑149 (MQ=255)
                                                                                                   ggAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGACGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGAt                                               >  8:248862/1‑149 (MQ=255)
                                                                                                   ggAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGACGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGAt                                               >  4:304148/1‑149 (MQ=255)
                                                                                                     aaGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTg                                             <  2:147373/149‑1 (MQ=255)
                                                                                                      aGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGa                                            <  6:216161/149‑1 (MQ=255)
                                                                                                          cGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGtc                                        <  7:84294/149‑1 (MQ=255)
                                                                                                           gTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGtc                                        >  4:251474/1‑148 (MQ=255)
                                                                                                               cTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTCGCGAGTCAGGAAGGGCGTAGATGCTGCCGAAACCTGGAAATTGTTGACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCaa                                   <  5:351947/149‑1 (MQ=255)
                                                                                                                    ccGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTCGAAATTGTTCACCCGACCATATTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCaa                                   >  7:115785/1‑144 (MQ=255)
                                                                                                                    ccGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCGGCCGCGACCTGGAAATGGGGCACCCGCGCAGCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCaataat                               >  5:86054/1‑148 (MQ=255)
                                                                                                                     cgcgGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGATGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTATGCCATCCAGCGTGCCCAGAACAGACTGATTGAGGTCCTCAataata                              >  4:335092/1‑148 (MQ=255)
                                                                                                                      gcgGCCAATTTATCGCGTCGCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAACTTGTTCTCCGGACCATCTACTTCCGTCTCTTTGACATCCAGCGTGCGCTGTACAGACTGATTGTGGTCATCAAtaataat                            <  8:157363/149‑2 (MQ=255)
                                                                                                                         gCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTg                                                               >  1:26850/1‑111 (MQ=255)
                                                                                                                         gCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTg                                                               <  2:26850/111‑1 (MQ=255)
                                                                                                                              tttATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGCCCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTg                    >  7:15137/1‑149 (MQ=255)
                                                                                                                              tttATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTg                    >  2:224752/1‑149 (MQ=255)
                                                                                                                                  taGCGCACCTGTCCCTGGCGAGTCAGGTAGGGCGTCGATGCTGCCGCAACCTGGAAATGGGTCACCCGCCCATCTTCTTCCGTCACTTTGACCTCCCGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCGGGCg                 <  3:128691/148‑1 (MQ=255)
                                                                                                                                  taGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCg                 <  3:123265/148‑1 (MQ=255)
                                                                                                                                     cGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCGGGAc             >  7:54140/1‑149 (MQ=255)
                                                                                                                                       aaCCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACgg                         >  2:159824/1‑135 (MQ=255)
                                                                                                                                       aaCCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACgg                         <  1:159824/135‑1 (MQ=255)
                                                                                                                                         ccTGTCCCTGGCGAGGCCGGAAGGGCGGCGAGGCTGCCGAAACCTGGAAATTGTTCCCCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGGCATCAATAATAAACGGGCCTGGCGGGACtttt         <  5:378610/149‑1 (MQ=255)
                                                                                                                                          cTGTCCCTGGCGAGTCAGGAAGGGCGTCGATTCTGCCAAAACCTTTAAAATCATCAAACGACAATCTTCTTCCGTCACTATTACAACCAACCTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCGCGACttttt        >  2:82888/1‑149 (MQ=255)
                                                                                                                                             tCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCGGGACTTTTTTCt     <  3:251474/149‑1 (MQ=255)
                                                                                                                                                cTGACGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCGGGACTTTTTTCTGGt  <  7:248862/149‑1 (MQ=255)
                                                                                                                                                cTGACGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCGGGACTTTTTTCTGGt  <  3:304148/149‑1 (MQ=255)
                                                                                                                                                 tGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCGGGACTTTTTTCTGGt  >  8:257294/1‑148 (MQ=255)
                                                                                                                                                 tGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCGGGACTTTTTTCTGGt  >  3:354785/1‑148 (MQ=255)
                                                                                                                                                 tGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACGTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCGGGACTTTTTTCTGGt  >  2:170942/1‑148 (MQ=255)
                                                                                                                                                 tGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCAGACAATCTTCTTCCGTCACTTTGACATCCAGCGTGACCTGAACAGACTGATTGAGGTCATCAATAATAACCGGGCCTGGCGGGACTTTTTTCTGGt  >  5:116729/1‑148 (MQ=255)
                                                                                                                                                   gCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACgg                         >  7:97617/1‑123 (MQ=255)
                                                                                                                                                   gCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACgg                         <  8:97617/123‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATAAACGGGCCTGGCGGGACTTTTTTCTGGT  >  NC_000913/751009‑751301

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: