Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 968,837 G→A D158N (GAT→AAT)  lpxK → tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913968,8370GA90.5% 105.8 / 3.7 42D158N (GAT→AAT) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (15/23);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

TATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGC  >  NC_000913/968610‑969066
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                     
tATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACAACCTGATTTGCATATCATCGTAACC‑GACAAAGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGGAATTTTCGTTATTGAt                                                                                                                                                                                                                >  1:98901/1‑251 (MQ=255)
      ggggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAAACTGTGTTGAATTTTCAACGCAATGATGCGCCTGTTGCGGTTTCTTCCGTTCGTTCTGATGCGGGAAAAACCCTTCTGGCGCCACCCCCCGATGGGCCGATCATCGTAAAC‑GACGACGGGTTTCAGCATTACCCGCTGGCGCGGAATGTGGAAATTGTCGTTTTTGTTgggggg                                                                                                                                                                                                          >  1:92658/1‑247 (MQ=255)
       gggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAAACTGTGTTGATTTATCAACGCCCTGGTGCACCTGGTTCGGGTTCTCCCCCTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCCACACCCTGGTGTGCAGATCATCGTAAAC‑GACGACGGGTTAACGCATTACCGTCTGGCGCGGAATGTGGGAAATGTCGTTATTTATGGGGTGc                                                                                                                                                                                                         >  1:276478/1‑251 (MQ=255)
       gggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGc                                                                                                                                                                                                         <  1:151762/251‑1 (MQ=255)
            atGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTCGGTTGTTTTATCAAAAAACCGATGCGCCCGTTGCGGTTTTTTCCGTTAGTTTTGAAGCGGTAAAAACCATTCTGGAGCAAACCCCTGATGGGCAGACTATCGAAACCCGA‑GAAGGTTTACAGAATTAACGGCTGGAGCGTAATGTGGGAATTGTCGTTATTGATGGTTTGCGTCGc                                                                                                                                                                                                    >  1:188314/1‑251 (MQ=255)
                tggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAAACTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGTCGTTATTTGTGGTGTTCGGCGCTTTg                                                                                                                                                                                                >  1:352379/1‑251 (MQ=255)
                 ggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTgg                                                                                                                                                                                               <  1:66499/251‑1 (MQ=255)
                 ggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCGGTGTTGATTTATCAACGCACTGATGCGCCTGTTGGGGTTTCGCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTATCC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTgg                                                                                                                                                                                               <  2:285451/251‑1 (MQ=255)
                  gtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGTTTTATTAACGCACTGATGCGCCTGTTGCGGGTTCTCCCCTTTGTTCTGATGCGGGAAAAGCCATTCTGGCGCCACACCCTGATGTGCAGATCATCGTAACC‑CACGACGGGTTACAGCATTACCGTCTGGCGCGGAATGTGGGAATTGTCGTTTTTTATGGGGTGCGTGGCTTTggg                                                                                                                                                                                              >  1:257103/1‑250 (MQ=255)
                   tAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGTAACTGTGTTGATTTTTCAACCCCCTGGTTAGCCTGTTGCGGTTTCTACCGTTCGTTATGATGCGGGAAAAGCCCTTCTGGCGCAAAACCCCGATGTGCAGATCATCGTAAAC‑GAAGACGGGTTACAGCATTACAGTCCGGATGGTAATGTGGAAATTGTGGTTTTTGATGGTGGTCGTCGGTTTGGc                                                                                                                                                                                              >  1:274522/1‑250 (MQ=255)
                     aGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAAt                                                                                                                                                                                           <  2:332976/251‑1 (MQ=255)
                                tATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAACCTGTGTTGATTTTTTAAAGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTTTGGTGGGGTAAAAGCCCTTCGTGCGCAAAACCCTGATGTGCCAATCATCGGAAAC‑GACGACGGTTTACAGCCTTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTTATGGTGTGGGTCGCTTTTGGAATGGCtggggggt                                                                                                                                                                                >  1:22999/1‑246 (MQ=255)
                                       tGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAAACTGTGTTGGTTTATCAAAGCACTGATGCGCCTGTTGCGGTTTCTCCCGGTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGTTGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGTCGTTATTTATGGTGTGCGTCGCTTTGGGAATGGCTGGGGGTTGgcggcg                                                                                                                                                                         >  1:411327/1‑251 (MQ=255)
                                        gTTATTTTCGGCAGATACCACAACAGCACAGGCGTGTGATGAACCTGTGTTGATTTATCAACGCACTGTTGCGCCTGTTGCGGTTTCTCCCGTGCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGAGCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCcggcg                                                                                                                                                                         <  2:309041/250‑1 (MQ=255)
                                         ttATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCCGTGTTGATTTATCAACGCACTGGAGCGCCTGTTGCGGTTTCTCCCCTTCGTTCTGATGCGGTAAAAACCCTTCTGGCGCAACACCCCGGTGGGCAGATCATCGTCACC‑GACGACGGGTTACCGCATTACCGTCTGGCGCGTAATGTGGAGATTGTCGTTATTTATGGTGTGCGGCGCTTTGGCCATGGCAGGGTGTTGACGGCggg                                                                                                                                                                       >  1:323905/1‑251 (MQ=255)
                                                                 gCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGgc                                                                                                                                                <  1:439067/250‑1 (MQ=255)
                                                                            tgatgaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGGGCGGGGCGCTTAAAGTCgg                                                                                                                                    >  1:340762/1‑251 (MQ=255)
                                                                               tgaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGAGCGAGCGGGGGGCTTAAAGTCGGTTg                                                                                                                                 >  1:97548/1‑251 (MQ=255)
                                                                                 aaCCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTgcg                                                                                                                               >  1:355248/1‑249 (MQ=255)
                                                                                  aCCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCAGTTGATg                                                                                                                              <  1:455793/251‑1 (MQ=255)
                                                                                     tgagtTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATAGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGGCCATGCGTGGGGGGGGGGGGCGCGTAAAGTCgggtgttgggg                                                                                                                           >  2:292541/4‑242 (MQ=255)
                                                                                      gtgtTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCgg                                                                                                                           <  1:208644/250‑1 (MQ=255)
                                                                                       ttttGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGc                                                                                                                                                               <  2:2807/211‑1 (MQ=255)
                                                                                       tgtTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCgggggccatgc                                                                                                                                                               >  1:2807/1‑206 (MQ=255)
                                                                                                   aaCGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAcggc                                                                                                             <  1:169056/251‑1 (MQ=255)
                                                                                                     cacctggATGCGCCTGTTGCGGTTTCTCCCGTTCGTGCGGATGCGGTAAAAGACATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAcggcgg                                                                                                           <  2:239294/245‑1 (MQ=255)
                                                                                                                 ccTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGAt                                                                                                                                                                                                                <  2:162844/138‑1 (MQ=255)
                                                                                                                 ccTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGAt                                                                                                                                                                                                                >  1:162844/1‑138 (MQ=255)
                                                                                                                 ccTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAg                                                                                               <  2:97548/251‑1 (MQ=255)
                                                                                                                         ggTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGGCAACGGCGGTGTTCCTccccggggggga                                                                                        >  1:193588/1‑239 (MQ=255)
                                                                                                                         ggTTTCTACCGTTCGGTCTGATGCGGTAAAAGTCATTCTGGCGCAACCACCGGATGGGCAGTTCAGCGGAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGGCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGaaa                                                                                       <  2:143234/251‑1 (MQ=255)
                                                                                                                          gTTTCTCCCGTTGTTTCGGAGGTGGTAAAAGCCATTCTGGCGCAACCCCCTGATGGGTAGATCATCGGAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCTCTTGGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGGCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGaaa                                                                                       <  2:258891/250‑1 (MQ=255)
                                                                                                                                  cgttcgttCTGATGCGGTAAAATGCATTCTGGCGCAACTCCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTTGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGGTGCGGTAATCTGCATCGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGc                                                                              <  1:130022/251‑1 (MQ=255)
                                                                                                                                       gttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATctg                                                                         <  1:206987/251‑1 (MQ=255)
                                                                                                                                       gtgctgttGCGGTAAACGCCATTCTGGCGCAACCCCCTGGTGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATctg                                                                         <  2:88948/244‑1 (MQ=255)
                                                                                                                                               gCGGTAAAAGCTATTCTGGCTCAACCCCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCgg                                                                  <  2:257103/250‑1 (MQ=255)
                                                                                                                                                         ccATTCTGGCGCAACACCCTGATGTGCATATCATCGTAACC‑GACAACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGTTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTg                                                       <  2:98901/251‑1 (MQ=255)
                                                                                                                                                           aTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGaa                                                     <  2:76917/251‑1 (MQ=255)
                                                                                                                                                                               ggtgCAGAGCATCGTAACC‑GACGTCGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGtt                                 <  2:355248/250‑1 (MQ=255)
                                                                                                                                                                                                cc‑GACGTCGGTTGACAGGATTACCGGCTGGCGCGTAAGGTGGAAATTGTCGGTATTGATGTTGTGCGTCGCTTTGGCAATGGCTGGTGGTGGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACTCGTTGTGACGTTGCTCAGCtt                <  2:53392/251‑1 (MQ=255)
                                                                                                                                                                                                            ttACAGCATTTCCGTCTGGCGCGTAATGTGGAAATTGTCGGTATTGATGGTGTGCGTCGCTTTGGCAAGGGCGGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCTGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGGTGCTCAGCTTGAACATGTAg      <  2:391719/250‑1 (MQ=255)
                                                                                                                                                                                                               cAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGGTGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGc  <  1:134279/251‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                     
TATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACC‑GACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGC  >  NC_000913/968610‑969066

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