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,8370GA94.4% 165.4 / ‑3.9 54D158N (GAT→AAT) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (17/34);  total (19/35)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGC  >  NC_000913/968596‑969072
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
tCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGaa                                                                                                                                                                                                                                      <  1:310327/250‑1 (MQ=255)
tCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTAGTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACTGCATTACCGTCTGGCGCGTAATGTGGaa                                                                                                                                                                                                                                      <  1:442953/250‑1 (MQ=255)
   cgTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAAACTGTGTTGATTTTTCAAAGCACCGAGGCGCCTGGTGCGGTTTCTTCCGCTCGTTTTGATGGGGGAAAAAACCATCTGGCGCAACACCCTGATGGGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTg                                                                                                                                                                                                                                  >  1:208287/1‑251 (MQ=255)
      gggggTCGTATCGCGGGGATAGGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCATCGCACTGATGTTCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTG‑‑GCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGtt                                                                                                                                                                                                                             <  1:428358/251‑1 (MQ=255)
           tCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAAACTGTGTTGATTTATCAACGCACTGATGCGCCTGTTTCGGTTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCCTCGTAACCGCCGACGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGTCGttttt                                                                                                                                                                                                                          >  1:221491/1‑248 (MQ=255)
               aTCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGAt                                                                                                                                                                                                                       <  1:286915/250‑1 (MQ=255)
                   cgGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAACCCGTGGTGATTTATTAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTTATGCGGTAAAAACCATTCTGGCGCAACACCCTGTTGGTCAGATCATCGTAACCCGCGACGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGTCGTTAtttttggtgg                                                                                                                                                                                                                  >  1:123551/1‑243 (MQ=255)
                      ggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAAAAGCACAGGCGGGTGATGAACCTGGGTTGATTTATTAACGCACTGTTGCGCCTGGTGCGGTTTCTCCCGGTCGTTTTGATGCGGGAAAAGCCCTTATTGCGCCACACCCTGATGTGCAGCACATCGTAACCGACGACGGGTTACAGCATTACCGTCTGGCGCGAAATGTGGGAATTGTCGGTATTTATGGTGTGCg                                                                                                                                                                                                               >  1:461635/1‑251 (MQ=255)
                       gATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAAACTGTGTTGGTTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCCGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTACCCTCTGGCGCGTGATGTGGAAATTTTCGGTATTGTTGGGGTGCgg                                                                                                                                                                                                              >  1:642167/1‑250 (MQ=255)
                              gggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTg                                                                                                                                                                                                       <  1:211593/250‑1 (MQ=255)
                               ggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGc                                                                                                                                                                                                                <  1:662160/241‑1 (MQ=255)
                                   aGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAAt                                                                                                                                                                                                  >  1:677095/1‑251 (MQ=255)
                                              tATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAAACTGTGTTGATTTATCAACCCACTGATGCGCCTGTTGCGGGTTCTCCCGGTTGTTCTGATGCGGTAAAAACCCTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACCACGGTTTACAACATTACCGACCTGCGCGTAATGGGGAAATTGTCGGTATTGATGGGGTGCGTTGCTTTGGCAATGGCTGGTGGtt                                                                                                                                                                                       >  1:410945/1‑251 (MQ=255)
                                                 ccGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGcc                                                                                                                                                                                    <  1:343774/251‑1 (MQ=255)
                                                    cTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGGACCTGGGTTTGTTTATTAACGCACAGTCGCACCTGTTGCGGGTTCTCCCGTTTGGTTTTTATAGGGAAAAAACATTTTGGGGCCAAACCCTGATGAGCAAAGCATCATAACCGACGACAGGTTTAAGCATTACCGTCTGACGCGGTAATAGGGAAATGTCGTATTTTATTGTGTGCGTCTCTTGTGCAATGGCTGGTGGGTGCaggc                                                                                                                                                                                 >  2:493915/1‑251 (MQ=255)
                                                            gTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAACCTGTGTTGGTTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCCGATGTGCAGATCCTCGTAAACGACGACGGTTTACAGCCTTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGGTGGTGTGCGTCGCTTTGGCAATGGCTGGGGGTTGCCGGCgggggccc                                                                                                                                                                         >  1:700069/1‑247 (MQ=255)
                                                                 cAGATACCACAACAGCACAGGCGGGTGATGAACCGGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTGCTGATGCGGTAAAAGCCATTCTGGCGCAACCCCCTGATGTGCAGATCATCGTAACCGACGTCGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAAGTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGt                                                                                                                                                                    <  2:150870/251‑1 (MQ=255)
                                                                      aCCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTgagc                                                                                                                                                                <  1:596202/250‑1 (MQ=255)
                                                                         acaacaGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTgagcgagc                                                                                                                                                            <  1:83354/251‑1 (MQ=255)
                                                                              agcccaGGCGGGGGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGgc                                                                                                                                                       <  2:623112/247‑1 (MQ=255)
                                                                               gCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCGCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGgcg                                                                                                                                                      <  1:603683/251‑1 (MQ=255)
                                                                               gCACAGGCGGGTGATGAACCGGTGTTGATTTATCAACGCACGGAGGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAAGTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGgcg                                                                                                                                                      <  2:133327/251‑1 (MQ=255)
                                                                                    ggCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGGGCTTa                                                                                                                                                  >  1:596711/1‑250 (MQ=255)
                                                                                              gaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCgggtgg                                                                                                                                       >  1:508313/1‑247 (MQ=255)
                                                                                               aaCCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGAt                                                                                                                                      <  1:456038/251‑1 (MQ=255)
                                                                                                    gtgtTGATTTATCTACTCCCTGTTGCGCCTGTTGCTGTTTCTCCCTTTCGTTCTGATGCGGTAATAGCCATTCTTGCGCTACTCCCTGCTTTGCAGCTCTTCGTATCCGACGTCGGTTTACAGCATTTCCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGTCAATTGCTGGTGGTTGCCTGCTGGTCCAATGCGTGAGCGATCTGGGCTCTTATACTCGGTTGATGCg                                                                                                                                   <  2:336910/249‑1 (MQ=255)
                                                                                                    gtggtGATTTATCAACGCATGGATGGTGTTGTTGCGGTTTCCCCCGTTCGTTCGGATGCGGTAAAAGCCAGTATGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCAGTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTGGCTGCggg                                                                                                                                 <  2:124647/247‑2 (MQ=255)
                                                                                                       ttGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTAATGCGGTAAAAGCCATTCTAGCGCAACCCCCTGATGTGCAGATCATCGTAACCGACGCCGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGGAATTGGCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCAGGGCCCAAGCGGGGGAGCGCGGGGCGCTTTAAGTCGGTTgcttcgggaa                                                                                                                               >  1:83599/1‑241 (MQ=255)
                                                                                                         gATTTATCAACGCACTGATGCGCCTGTTGCGTTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATc                                                                                                                             <  2:545629/250‑1 (MQ=255)
                                                                                                           tttATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGt                                                                                                                           <  2:221491/250‑1 (MQ=255)
                                                                                                                 aaCGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTgag                                                                                                                                                                 <  1:199911/206‑1 (MQ=255)
                                                                                                                     cACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGgtg                                                                                                                <  1:704637/251‑1 (MQ=255)
                                                                                                                      aCTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTTCCGGCGGGGCACATGTGTGTGCGAGCGGGGCGCTTAAAGATGGTTGTTGCGGtgttagtgtacgcgggggg                                                                                                               >  2:325218/1‑233 (MQ=255)
                                                                                                                         gATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTccc                                                                                                           <  1:154040/251‑1 (MQ=255)
                                                                                                                         gATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTccc                                                                                                           <  1:595195/251‑1 (MQ=255)
                                                                                                                          aTGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGAGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGGGCGAGCGGGGCGCTTTTAGGCGGTTGATGCGGTAATCCTCAAcggcgggggtccc                                                                                                           >  1:566953/1‑244 (MQ=255)
                                                                                                                                cTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGGGCGAGCGGGGCGCTTAAAGTCGGTTGATGGGGTAATCGTCAAcggcgggggcccccgccgc                                                                                                     >  1:179455/1‑238 (MQ=255)
                                                                                                                                 tGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTaaa                                                                                                                                                <  1:233934/207‑1 (MQ=255)
                                                                                                                                             tcCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTc                                                                                                                                             <  1:647187/198‑1 (MQ=255)
                                                                                                                                                 gttcgttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCa                                                                                   <  1:231058/251‑1 (MQ=255)
                                                                                                                                                              cccccacacgacgctcttccgatcccacCCTGATGTGCAGATCATCGTAACCGACGACGGTTGACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGGCGCTGTGGCAATGGCTGGTGGTTGCCGGCGGTGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTTATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCgg                                                                        <  2:76808/226‑1 (MQ=255)
                                                                                                                                                                     aaaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGTGCCCCCCCAGCGGGGGAATACCCATGCATCTTCTGCCGGGGCAggg                                                                 >  1:318488/1‑250 (MQ=255)
                                                                                                                                                                            ttCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCCGCGGTGGAATACCCATTCATCTTCTGCCGGGtgcggcggggagt                                                          >  1:627601/1‑239 (MQ=255)
                                                                                                                                                                                ggCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCTTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGTTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCTCTTAAAGTCGGTTGATGCTGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAAttt                                                        <  1:646607/249‑1 (MQ=255)
                                                                                                                                                                                ggCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTTAAGTCGGTTGATGCGGGAATCGTCAACGGCGGTGTCCCTCGCAGCGGGGGAATCCCCCCGCCTCTGCCGCCggggccgggggggacatttc                                                      >  1:84571/1‑234 (MQ=255)
                                                                                                                                                                                             gATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTAcg                                           <  2:700069/249‑1 (MQ=255)
                                                                                                                                                                                                        catcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGGCCCCCGCAGCGGTGAAAACCCCCTGCATCTGCTGCCGGGTCCGGCGGGGGATTTTCATACCGGGAAGCGTTGTGACggtg                              >  1:522382/1‑248 (MQ=255)
                                                                                                                                                                                                                         ggTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGt              <  1:282314/251‑1 (MQ=255)
                                                                                                                                                                                                                         ggTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGt              <  1:239554/251‑1 (MQ=255)
                                                                                                                                                                                                                              aCAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTg          <  2:208287/250‑1 (MQ=255)
                                                                                                                                                                                                                              aCAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTg          <  1:125458/250‑1 (MQ=255)
                                                                                                                                                                                                                              aCAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTAGTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCg                                                                                                   <  2:567007/161‑1 (MQ=255)
                                                                                                                                                                                                                                aGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCg       <  1:697060/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aCCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGc  <  1:448994/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
TCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGC  >  NC_000913/968596‑969072

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