Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,964 C→T 5.4% L196L (CTG→CTA aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913120,9640CT 5.4% 243.2 / 2.7 74L196L (CTG→CTAaroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base C (35/35);  new base T (2/2);  total (37/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.62e-01

GATGCGGTAGATAACCTGGTTAGTTGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACT  >  NC_000913/120817‑121111
                                                                                                                                                   |                                                                                                                                                   
gATGCGGTAGATAACCTGGTTAGTTGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCa                                                                                                                                                    <  3:389667/149‑1 (MQ=255)
          aTAACCTGGTTAGTTGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAaccacc                                                                                                                                          <  2:438136/149‑1 (MQ=255)
              ccTGGTTAGTTGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGa                                                                                                                                      <  3:213028/149‑1 (MQ=255)
               cTGGTTAGTTGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGAt                                                                                                                                     >  1:558690/1‑149 (MQ=255)
               cTGGTTAGTTGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGAt                                                                                                                                     >  4:217334/1‑149 (MQ=255)
                    tAGTTGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCcac                                                                                                                                <  2:474674/149‑1 (MQ=255)
                      gTTGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAg                                                                                                                              >  4:4093/1‑149 (MQ=255)
                       ttGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAgg                                                                                                                             >  4:568481/1‑149 (MQ=255)
                       ttGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAgg                                                                                                                             >  3:340219/1‑149 (MQ=255)
                        tGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGt                                                                                                                            <  3:568481/149‑1 (MQ=255)
                           tttCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCcaca                                                                                                                               <  2:198162/143‑1 (MQ=255)
                           tttCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCcaca                                                                                                                               >  1:198162/1‑143 (MQ=255)
                           tttCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGc                                                                                                                         <  3:201376/149‑1 (MQ=255)
                               tgTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCAGTGATCCCCACCAGTTCCTGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGTAACCACCCTGATCCCACATGTTGCTTAc                                                                                                                     <  4:62718/148‑1 (MQ=255)
                                gTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCTTCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGAt                                                                                                                                     >  1:270375/1‑132 (MQ=255)
                                gTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCTTCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGAt                                                                                                                                     <  2:270375/132‑1 (MQ=255)
                                   taCTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTc                                                                                                                 >  3:76990/1‑149 (MQ=255)
                                      tttGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGcc                                                                                                              >  3:363909/1‑149 (MQ=255)
                                           tCCGGGTTATCAGCTTCTGCTGCGGTGATCGCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCgg                                                                                                         >  2:679843/1‑149 (MQ=255)
                                           tCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGTAGAAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCgg                                                                                                         >  2:710860/1‑149 (MQ=255)
                                           tCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGTAGAAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCgg                                                                                                         >  4:53592/1‑149 (MQ=255)
                                             cGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGgc                                                                                                       >  4:616284/1‑149 (MQ=255)
                                                gTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGgccgc                                                                                                    >  4:323049/1‑149 (MQ=255)
                                                    tCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGtt                                                                                                <  1:186572/149‑1 (MQ=255)
                                                    tCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGtt                                                                                                >  2:473732/1‑149 (MQ=255)
                                                    tCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGtt                                                                                                <  2:466621/149‑1 (MQ=255)
                                                    gcAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGACCGCCGtt                                                                                                >  2:656164/2‑149 (MQ=255)
                                                        cTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCa                                                                                            >  1:619667/1‑149 (MQ=255)
                                                        cTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCa                                                                                            <  3:217334/149‑1 (MQ=255)
                                                        cTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCa                                                                                            >  2:642958/1‑149 (MQ=255)
                                                             gctgcGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAGCATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGaa                                                                                       <  3:86281/149‑1 (MQ=255)
                                                              ctgcGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAAt                                                                                      <  3:51269/149‑1 (MQ=255)
                                                                  ggTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATagca                                                                                  <  2:415907/149‑1 (MQ=255)
                                                                  ggTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATagca                                                                                  <  1:484472/149‑1 (MQ=255)
                                                                      aTCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGcagcca                                                                              >  4:93450/1‑149 (MQ=255)
                                                                         cccACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCa                                                                                            >  4:3019/1‑132 (MQ=255)
                                                                         cccACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCa                                                                                            <  3:3019/132‑1 (MQ=255)
                                                                          ccaccaGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAgccg                                                                          >  1:197205/1‑149 (MQ=255)
                                                                          ccaccaGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAgccg                                                                          >  1:356774/1‑149 (MQ=255)
                                                                          ccaccaGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAgccg                                                                          >  2:401888/1‑149 (MQ=255)
                                                                          ccaccaGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAgccg                                                                          >  2:370059/1‑149 (MQ=255)
                                                                            accaGTTCCAGACCATCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAgccgcc                                                                        >  4:558819/1‑149 (MQ=255)
                                                                                   ccAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATcc                                                                                                                                   <  4:270582/83‑1 (MQ=255)
                                                                                   ccAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATcc                                                                                                                                   >  3:270582/1‑83 (MQ=255)
                                                                                     aGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAgatgat                                                               <  2:558690/149‑1 (MQ=255)
                                                                                       accaccGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGaa                                                                     <  2:448215/141‑1 (MQ=255)
                                                                                       accaccGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGaa                                                                     >  1:448215/1‑141 (MQ=255)
                                                                                          accGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATagcag                                                                                 >  2:356457/1‑126 (MQ=255)
                                                                                          accGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATagcag                                                                                 <  1:356457/126‑1 (MQ=255)
                                                                                          accGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCg                                                          >  2:175769/1‑149 (MQ=255)
                                                                                                 gagaACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTaccac                                                   >  4:413331/1‑149 (MQ=255)
                                                                                                    aaCATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACcgc                                                <  1:370059/149‑1 (MQ=255)
                                                                                                       aTGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAgatg                                                                 <  1:478000/129‑1 (MQ=255)
                                                                                                       aTGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAgatg                                                                 >  2:478000/1‑129 (MQ=255)
                                                                                                         gATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATcc                                                                                                                                   <  1:603788/61‑1 (MQ=255)
                                                                                                         gATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATcc                                                                                                                                   >  2:603788/1‑61 (MQ=255)
                                                                                                            aaTCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACtt                                        >  3:537093/1‑149 (MQ=255)
                                                                                                             aTCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGTAGAAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACttt                                       <  3:53592/149‑1 (MQ=255)
                                                                                                             aTCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGTAGAAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACttt                                       <  1:710860/149‑1 (MQ=255)
                                                                                                             aTCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACttt                                       <  1:280577/149‑1 (MQ=255)
                                                                                                             aTCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACttt                                       <  2:404970/149‑1 (MQ=255)
                                                                                                                gCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGAt                                    <  4:7843/149‑1 (MQ=255)
                                                                                                                  catcatcatCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATaa                                  >  1:382969/1‑149 (MQ=255)
                                                                                                                    tcatcatcaCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATg                                <  3:4093/149‑1 (MQ=255)
                                                                                                                    tcatcatcaCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATg                                <  3:323049/149‑1 (MQ=255)
                                                                                                                    tcatcatcaCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATg                                >  3:241444/1‑149 (MQ=255)
                                                                                                                    tcatcatcaCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATg                                <  1:175769/149‑1 (MQ=255)
                                                                                                                           caCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGa                                     >  3:8763/1‑137 (MQ=255)
                                                                                                                           caCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGa                                     <  4:8763/137‑1 (MQ=255)
                                                                                                                                       gAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAACAGCCCGCCTCTGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACACCAtct             <  2:252987/149‑1 (MQ=255)
                                                                                                                                        aaGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCAtctc            <  3:403992/149‑1 (MQ=255)
                                                                                                                                              tGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCaaa      <  4:340219/149‑1 (MQ=255)
                                                                                                                                              tGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCaaa      <  4:363909/149‑1 (MQ=255)
                                                                                                                                                  gCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACt  >  1:726539/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATGCGGTAGATAACCTGGTTAGTTGCTTTCGGTATACTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACT  >  NC_000913/120817‑121111

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: