Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,111,727 2 bp→TT coding (865‑866/2544 nt) opgH → osmoregulated periplasmic glucans biosynthesis protein H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,111,7270CT98.5% 213.0 / NA 65P289M (CCG→TTG) opgHosmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/29);  minor base G (0/1);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009131,111,7280CT98.5% 196.8 / NA 65P289M (CCG→TTG) opgHosmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/29);  minor base G (0/1);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGA  >  NC_000913/1111501‑1111964
                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                            
aaTGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCaa                                                                                                                                                                                                                        >  1:221342/1‑250 (MQ=255)
   gggCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCaaaaa                                                                                                                                                                                                                     <  1:355635/250‑1 (MQ=255)
           tGCAACTGCTTCTTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGCAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCtggatg                                                                                                                                                                                                            >  1:172530/1‑251 (MQ=255)
             cAACTGCTTATTGGTCGCGATAAATACAGTATATTTGCGTCAACAGTTGGCGATGACCCAGTAAACCCGGAGCATCGCACGGCGCTGATCATGCCTATCTGTAACGAAGACGTGAGCCGTGGTATTGCTGGCCAGCCTGCAACGTGGGAATCAGTAAAAGCCACCAGGAAAGCCAAACACGATGATGTCAACATTCTAATTGACAGTGATAACTTTGATATCTGCGTCGCAGAACAAAAAGCCtcgatgga                                                                                                                                                                                                          <  2:456021/251‑1 (MQ=255)
                  gCTTATTGGTCGCGATAAATACAGTAGATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCAGCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCt                                                                                                                                                                                                       <  2:221342/249‑1 (MQ=255)
                   gttATTGGTCGCGAGAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGATCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGTGAATCATGAAAAGCCTCCGTGAATGCCAAACACGTTGATGTCTACATTCTTAGTGACGGTTATATCTTGGATAGCTGCGGTGCAGAGCGAAAAGCCttggtggttctgat                                                                                                                                                                                                    <  2:479993/250‑14 (MQ=255)
                     tATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCg                                                                                                                                                                                                  >  2:121998/1‑251 (MQ=255)
                        tGGTCGCGATAAATACCGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGCGCTTATCGCTg                                                                                                                                                                                               >  2:426938/1‑251 (MQ=255)
                         ggTCGCGATAAATACAGTATATCTGCGTCCACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTg                                                                                                                                                                                               >  1:177676/1‑250 (MQ=255)
                         ggTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTg                                                                                                                                                                                               <  1:487451/250‑1 (MQ=255)
                                   aaTACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACTTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGc                                                                                                                                                                                    <  2:120057/251‑1 (MQ=255)
                                    aTACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGCAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGACGTCGGTTGCg                                                                                                                                                                                   >  1:417510/1‑251 (MQ=255)
                                        aGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGAATCGCACGGCGTTGATCATGCCTATCTGTAAAGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCAACCGGGAATGCCAAACCCTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAgg                                                                                                                                                                               >  1:219510/1‑251 (MQ=255)
                                           atatCTGCGTCAACAGTTGGTGATGAACCAGTAAACCCGGAGCATCGCACGGCGTTGATCATTCCTATCTGTAACGAAGACGTGAACCGTGTTTGTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGGCAGTGATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCGGAAGTCGGTGGCGAAGGTCa                                                                                                                                                                            <  2:328766/251‑1 (MQ=255)
                                           atatCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACAGTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCa                                                                                                                                                                            <  2:428753/251‑1 (MQ=255)
                                              tcagcggcAACAGATGGAGATGAACCATTCCAACCGGAGCATCGCACGGCGTTGATCATGCCAATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGACTGCGTGCAACGTGGGAATCAGTAACAGCCACCGGGAATGCCCAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAGGCATGGAAGGAGCTTATCGCTGAAGTCGGAGGCGAAGGTc                                                                                                                                                                             >  2:475766/8‑247 (MQ=255)
                                                  cGTCAACAGTTGGCGATGCACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTc                                                                                                                                                                     >  2:313026/1‑251 (MQ=255)
                                                  cGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTc                                                                                                                                                                     <  1:18650/251‑1 (MQ=255)
                                                   gTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCATGGATGGAGCTTATCGCTGCAGTCGGTGGCGAAGGGCAGCTTTTCt                                                                                                                                                                    >  2:260476/1‑251 (MQ=255)
                                                        cAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCCTGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCCAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGCTTTTCTATcgc                                                                                                                                                               >  2:505795/1‑251 (MQ=255)
                                                         agctagCGATGAACCATTAAACACGGATCGTCGCTCGGTGTTGATCATGCCTATCTGTAACGAAGACGGGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGTAATCAGTAAGAGCCACCGGGAATGCCAAACACTTTGATGTCTACATGCTTAGTGACAGTCATAACTTGGATATCTGCGTCGCAGATCAAAAAGGCTGGATGGAGATTATCACTGAAGTCGGTGGCGAAGGTCACATTTTCTATcgca                                                                                                                                                              <  2:91876/246‑2 (MQ=255)
                                                             ggCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATcgccgcc                                                                                                                                                           <  1:320409/250‑1 (MQ=255)
                                                                gATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGccgccgccg                                                                                                                                                       >  1:120254/1‑251 (MQ=255)
                                                                  tGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTc                                                                                                                                                     >  1:389613/1‑251 (MQ=255)
                                                                     aCCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTTAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTcgcc                                                                                                                                                  <  1:426938/251‑1 (MQ=255)
                                                                        aTTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGTGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCcgcg                                                                                                                                               <  1:52766/251‑1 (MQ=255)
                                                                                        cGCACTGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCg                                                                                                                                <  1:357946/250‑1 (MQ=255)
                                                                                         gCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCGACCGGGAATGCCAAACACTTGGATGGCGACAGTCTTAGTGACAGGTATAACTTGGATATCTGCGGCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGGAGTCGGTgg                                                                                                                                                                                     <  2:238863/196‑1 (MQ=255)
                                                                                         gCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGAAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGTTGGAGCTTATCGCTGAAGTCGGTgg                                                                                                                                                                                     >  1:238863/1‑196 (MQ=255)
                                                                                         gCACGGCGTGGATCGTGCCTATCTGTAACGTAGACGTTAACCGTGTTTTTGCTGTCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTGATAACTTGGATATCTGCGGCGCAGAGCAAAAAGCCCGGATGGGGCGTATCTGTGAAGTCGGTGGCGAAGGTCAGAGTTTCTATCGTCGCCGCCGTAGCCGCGGGAAGCGTAAAAGCGGt                                                                                                                              <  2:178104/251‑1 (MQ=255)
                                                                                               cGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATc                                                                                                                        <  1:245196/251‑1 (MQ=255)
                                                                                                  tGATCATGCCTATCTGTAAAGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGACACCGGGAATGCCAAACACTTTGATGTATACATTCTTAGTGACAGTTATAACTTGGATATATGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGAt                                                                                                                     >  1:306177/1‑251 (MQ=255)
                                                                                                    aTCAGGCCTATCGGTAACGAAGACGTGAACCGTGATGTGGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGACACCGGGAAGGCAAAACACTTTGGTGTCTACATTCGTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTTGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGa                                                                                                                   <  1:121057/251‑1 (MQ=255)
                                                                                                         gCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGGAAAAGCGGTAATATCGATGACTTCt                                                                                                              <  2:432032/251‑1 (MQ=255)
                                                                                                               cTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTc                                                                                                        >  1:52314/1‑251 (MQ=255)
                                                                                                                 gTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTc                                                                                                        <  2:422589/249‑1 (MQ=255)
                                                                                                                        aGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGc                                                                                                >  1:358976/1‑250 (MQ=255)
                                                                                                                                aCCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTAc                                                                                       >  1:264945/1‑251 (MQ=255)
                                                                                                                                 ccGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTTGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATAGCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACa                                                                                      <  2:264945/251‑1 (MQ=255)
                                                                                                                                   gtgtTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGc                                                                                    >  1:128532/1‑251 (MQ=255)
                                                                                                                                     gtTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCt                                                                                   >  2:337853/1‑250 (MQ=255)
                                                                                                                                          tGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACAtgg                                                                             >  2:452605/1‑251 (MQ=255)
                                                                                                                                                   gCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGTCGCCGCCGTCGTCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTgg                                                                    <  1:219145/251‑1 (MQ=255)
                                                                                                                                                   gCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTg                                                                     <  1:55335/250‑1 (MQ=255)
                                                                                                                                                   gCGTGAAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACCTTGATGTCTACATTCTTAGTGACAGTTATACCTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGATCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTTCCGTCGCGGGGGCAGCCAGTACAGCTACATGGTGGTGCTgg                                                                    <  2:339687/251‑1 (MQ=255)
                                                                                                                                                    cGTGCAACGAGGGAATCAGTAAACGCCACCGGGAATGCCAAACACTTTGATGTCTACATTATTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTgg                                                                    >  2:318192/1‑250 (MQ=255)
                                                                                                                                                         aaCGTGGAAATCAGTAAAAGCCACCGGGAATGCAAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGCCGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTCATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTg                                                              >  2:331043/1‑251 (MQ=255)
                                                                                                                                                                   tCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGCGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAAt                                                    >  1:328272/1‑251 (MQ=255)
                                                                                                                                                                     aGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATg                                                   >  2:166660/1‑250 (MQ=255)
                                                                                                                                                                        aaaaGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCg                                               >  1:57856/1‑251 (MQ=255)
                                                                                                                                                                              cACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGAtt                                         >  1:405081/1‑251 (MQ=255)
                                                                                                                                                                              cACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTGGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTGAGCGCTGAAGTCGGGGGCGAAGGGCAGATTTTCTATCGCCGCCGCCGGCGCCTCGTGAAGCGTAAAAGCTGGAATATCGATGACTTCTGCCGGCGCTGGGTCAGCCAGTACAGCTACATGGTGGTGCTGGAGGTTGTCTCTGTGAGGACCGGTGAtt                                         <  2:498106/251‑1 (MQ=255)
                                                                                                                                                                               aCCGGCAATGCAAGACAATTTGCTGTCTACATTCTTAGTGACCGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGAGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTAAAGCTACCTGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGAttg                                        >  2:460590/1‑251 (MQ=255)
                                                                                                                                                                                  gggAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGAttgttt                                     >  1:451325/1‑251 (MQ=255)
                                                                                                                                                                                     aaTGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTtgtg                                  >  1:180799/1‑251 (MQ=255)
                                                                                                                                                                                        gCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCg                                <  2:358976/250‑1 (MQ=255)
                                                                                                                                                                                             actcTTTGAGGTATACATTCTTAGTGACAGGTATAACTTGGAAATCTGCGGCGCAGGGCAAAAAGCGTGGATGGGGCTTATCGCTGAAGTCGGTGGCGAAGGGCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGGTTGGGCGGGCTgg                          <  2:177676/248‑1 (MQ=255)
                                                                                                                                                                                                      tGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTg                  >  2:393825/1‑250 (MQ=255)
                                                                                                                                                                                                         cTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGGGATTGTTTGTGCGGGCTGGTGCGCCTGATgg              >  1:14724/1‑251 (MQ=255)
                                                                                                                                                                                                            cATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACCGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGa             >  2:140524/1‑249 (MQ=255)
                                                                                                                                                                                                                  tAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAAc      <  1:51167/250‑1 (MQ=255)
                                                                                                                                                                                                                  tAGTGAAAGTTATAACGTGGATATCTGCGTCTCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTAACGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTTCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGTTAATGACAGGTGATTGTTTGTGCGGGCTGGTGCGCTTGATGGAAGCCAAc      <  2:334251/250‑1 (MQ=255)
                                                                                                                                                                                                                    gTGACCGTTATAACTTGGAGATCTGCGTCGCAGAGCACAAAGCCTGGATGGAGCTTATCGCTGCCGACGGTGGCGGAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGACGCTGAGGCAGCCAGTAAAGCTACCTGGTAGTGCTGGATGCTGACTCGGTAATGTCAGGGGATTGTTTTTGTGGGCTGTTGCGCATGACAGAAGCCAACCCg   >  2:412699/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGGCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGACGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCg   >  1:499437/1‑250 (MQ=255)
                                                                                                                                                                                                                      gACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGa  <  1:384104/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                            
AATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGA  >  NC_000913/1111501‑1111964

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