Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 522,634 G→A 5.0% G740D (GGC→GAC)  ybbP → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913522,6340GA 5.0% 254.4 / 2.3 79G740D (GGC→GAC) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (40/35);  new base A (2/2);  total (42/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.67e-01

TGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGT  >  NC_000913/522486‑522782
                                                                                                                                                    |                                                                                                                                                    
tGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCgg                                                                                                                                                      >  5:202771/1‑149 (MQ=255)
   ggTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCt                                                                                                                                                   >  7:208471/1‑149 (MQ=255)
        gttgctgttgctAGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCGGAGCTGGTGGTGTGGCGCACACTCGGGGCGGGGAAAAAAAGGCTGCGTAACACGTTGTGGTGTGACTTCGCCAAGGGTGGGTTTGTTTCCGGCCTGGTGGc                                                                                                                                             >  8:72178/1‑149 (MQ=255)
              gttgctgGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCaa                                                                                                                                        >  8:108944/1‑149 (MQ=255)
               ttgctgGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAAt                                                                                                                                       >  2:416519/1‑149 (MQ=255)
               ttgctgGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAAt                                                                                                                                       >  4:318182/1‑149 (MQ=255)
                 gctgGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTg                                                                                                                                     >  4:406124/1‑149 (MQ=255)
                  ctgGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTg                                                                                                                                     <  6:335193/148‑1 (MQ=255)
                     gctcaggtgcaggtgGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGc                                                                                                                                 <  6:263532/146‑1 (MQ=255)
                       acaGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGc                                                                                                                                 >  1:135953/1‑147 (MQ=255)
                       acaGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGc                                                                                                                                 <  2:135953/147‑1 (MQ=255)
                       acaGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCg                                                                                                                                >  1:101712/1‑148 (MQ=255)
                        caggtgcaggtgGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCgg                                                                                                                               >  7:451582/1‑148 (MQ=255)
                                          cgtcagcgtcaTCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAgtgt                                                                                                             <  2:49488/148‑1 (MQ=255)
                                             cagcgtcaTCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCa                                                                                                         <  4:338566/149‑1 (MQ=255)
                                                 gtcaTCAGGAGCTGGTGGTGTGGCGCACAATCGGTGCGGGGAAAAAACTGCTGCGTACCCCGTTGTTGTGTGCGTTCGCCATGCGTGGGTTTGTTTCCGGCCTGGTGGCCGCAATAGGAGCGGAAACGGCGCTGGCAGTGGTGAAGGc                                                                                                      >  4:23864/1‑148 (MQ=255)
                                                    atcaGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGaaa                                                                                                  <  3:406124/149‑1 (MQ=255)
                                                       aGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAgtg                                                                                               >  2:380231/1‑149 (MQ=255)
                                                        ggAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAgtgt                                                                                              <  3:172938/149‑1 (MQ=255)
                                                          aGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGttt                                                                                            <  1:111166/149‑1 (MQ=255)
                                                            cTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGACCTGGTGGCCGCAATTGGTGCGGAAACGgcgc                                                                                                                      >  5:430545/1‑121 (MQ=255)
                                                            cTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGACCTGGTGGCCGCAATTGGTGCGGAAACGgcgc                                                                                                                      >  5:437348/1‑121 (MQ=255)
                                                            cTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGACCTGGTGGCCGCAATTGGTGCGGAAACGgcgc                                                                                                                      <  6:430545/121‑1 (MQ=255)
                                                            cTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGACCTGGTGGCCGCAATTGGTGCGGAAACGgcgc                                                                                                                      <  6:437348/121‑1 (MQ=255)
                                                               gtggtgTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGAttt                                                                                       <  2:321957/149‑1 (MQ=255)
                                                               gtggtgTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGAttt                                                                                       >  2:131074/1‑149 (MQ=255)
                                                                  gtgtGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTcc                                                                                    >  5:190104/1‑149 (MQ=255)
                                                                    gtgGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGt                                                                                  >  1:316013/1‑149 (MQ=255)
                                                                     tgGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGt                                                                                  <  2:316013/148‑1 (MQ=255)
                                                                     tgGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTg                                                                                 <  8:451582/149‑1 (MQ=255)
                                                                     tgGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTg                                                                                 <  1:131074/149‑1 (MQ=255)
                                                                       gcgcACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTgg                                                                                <  1:416519/148‑1 (MQ=255)
                                                                         gcACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTgg                                                                                                                   >  7:152588/1‑111 (MQ=255)
                                                                         gcACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTgg                                                                                                                   <  8:152588/111‑1 (MQ=255)
                                                                           acacTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGc                                                                            >  7:347795/1‑148 (MQ=255)
                                                                                 ggTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGAtt                                                                                        >  2:450758/1‑130 (MQ=255)
                                                                                 ggTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGAtt                                                                                        <  1:450758/130‑1 (MQ=255)
                                                                                 ggTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGAc                                                                       <  6:351944/147‑1 (MQ=255)
                                                                                 ggTGCGGGGAAAAAACTGCTGCGGACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACt                                                                      <  1:288814/148‑1 (MQ=255)
                                                                                   tGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCTTGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGc                                                                   <  6:202771/149‑1 (MQ=255)
                                                                                      ggggAAAAAACTGCTGCGGACCACGTTGTGGTGTGAGTTCGCCATGCTGGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGAt                                                                <  5:312740/149‑1 (MQ=255)
                                                                                       gggAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTCGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGAtt                                                               >  1:59191/1‑149 (MQ=255)
                                                                                        ggAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATtg                                                              >  6:374750/1‑149 (MQ=255)
                                                                                        ggAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATtg                                                              >  5:214958/1‑149 (MQ=255)
                                                                                         gAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATtgt                                                             <  2:39963/149‑1 (MQ=255)
                                                                                             aaaCTTCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGGTTGATTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGAt                                                         <  1:304919/149‑1 (MQ=255)
                                                                                               aCTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTg                                                       >  2:120485/1‑149 (MQ=255)
                                                                                                   ctgcGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGtt                                                     <  3:318182/147‑1 (MQ=255)
                                                                                                       gTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCg                                               >  7:419393/1‑149 (MQ=255)
                                                                                                       gTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGAATTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCg                                               <  6:190104/149‑1 (MQ=255)
                                                                                                          ccACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTg                                             >  2:5761/1‑148 (MQ=255)
                                                                                                           cACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCa                                           <  8:208471/149‑1 (MQ=255)
                                                                                                            aCGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAg                                          <  6:214958/149‑1 (MQ=255)
                                                                                                               ttGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCg                                        >  8:361210/1‑148 (MQ=255)
                                                                                                                tgtgGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGc                                                                            <  8:139406/111‑1 (MQ=255)
                                                                                                                tgtgGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGc                                                                            >  7:139406/1‑111 (MQ=255)
                                                                                                                 gtggtgTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGa                                      <  2:101712/148‑1 (MQ=255)
                                                                                                                    gtgtGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATtgt                                                             >  4:96533/1‑122 (MQ=255)
                                                                                                                    gtgtGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATtgt                                                             <  3:96533/122‑1 (MQ=255)
                                                                                                                     tgtgAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCgct                                 >  4:38979/1‑149 (MQ=255)
                                                                                                                     tgtgAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATGGTGGATTGTTCTGCCGTGCAGCGGAGCgct                                 >  6:247243/1‑149 (MQ=255)
                                                                                                                      gtgAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCgctg                                >  6:416603/1‑149 (MQ=255)
                                                                                                                         aGTTCGCCATGCTTGGGTTTGTTTCCGGCCAGGTGGCCGCAAATGGTGCGGAAACGGAGCTGGCAGTGTTGCAGGCGAAAGTGTTTGGTTTTCCGTGGGAGCAAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCgctgctg                             >  5:403342/1‑149 (MQ=255)
                                                                                                                           ttCGGCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCgctgctgct                           >  5:40887/1‑149 (MQ=255)
                                                                                                                              gCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCgctgctg                             >  5:353803/1‑144 (MQ=255)
                                                                                                                              gCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCgctgctg                             <  6:353803/144‑1 (MQ=255)
                                                                                                                              gCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGt                         <  1:380231/148‑1 (MQ=255)
                                                                                                                              gCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGt                         <  8:419393/148‑1 (MQ=255)
                                                                                                                              gCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTc                        >  8:436933/1‑149 (MQ=255)
                                                                                                                                   gCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTc                                                                                     >  8:82060/1‑83 (MQ=255)
                                                                                                                                   gCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTc                                                                                     <  7:82060/83‑1 (MQ=255)
                                                                                                                                     ttGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTc                 >  2:165388/1‑149 (MQ=255)
                                                                                                                                       gggTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTcg                <  1:5761/148‑1 (MQ=255)
                                                                                                                                         gtttgtttCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTcggcg             >  6:407117/1‑149 (MQ=255)
                                                                                                                                         gtttgtttCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTGGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTcggcg             >  1:291869/1‑149 (MQ=255)
                                                                                                                                          tttgtttCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGGTTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTcggcgg            <  6:403342/149‑1 (MQ=255)
                                                                                                                                          tttgtttCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTcggcgg            >  4:291727/1‑149 (MQ=255)
                                                                                                                                             gtttCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCggctg         >  5:380029/1‑149 (MQ=255)
                                                                                                                                                 ccGGCCTGGTGGCCGCAATTGGTGCGGAACCGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCggctggct      <  7:108944/148‑1 (MQ=255)
                                                                                                                                                 ccGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGgctggctg     >  3:65742/1‑149 (MQ=255)
                                                                                                                                                    gCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGt  <  3:291727/149‑1 (MQ=255)
                                                                                                                                                    gCCTGGTGGCCGCAATTGGAGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGCGAGCCAGACTGGCGATTGTGGATTGGTCTGCCGTGCAGCGGAGCGCTGATGCTGTCGCGTTTCGGCGGCTGGCTGGGt  >  2:183206/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGTTGTGGTGTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGT  >  NC_000913/522486‑522782

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