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,8370GA98.0% 160.2 / NA 52D158N (GAT→AAT) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/28);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCG  >  NC_000913/968602‑969073
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
ggggggCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCGCAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTc                                                                                                                                                                                                                                 <  2:293214/250‑1 (MQ=255)
 ggggTCGGATCGTGGGTATATGGTGTTAAGGCTGAATCTTATCCGCTGTTATTGTCGTCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGTCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTTGCGCGTAATG‑TGGAAATTGTCGt                                                                                                                                                                                                                               <  1:409483/251‑1 (MQ=255)
    gTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGTTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTAt                                                                                                                                                                                                                            <  1:89937/251‑1 (MQ=255)
      cGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTg                                                                                                                                                                                                                          <  1:58834/251‑1 (MQ=255)
           cgctgggATATGGTGGTAAGGCTGAAGCTTATCCGCTGTTATTTTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATgg                                                                                                                                                                                                                      <  1:64372/246‑1 (MQ=255)
            gcgGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGgtg                                                                                                                                                                                                                    <  1:31734/251‑1 (MQ=255)
              ggggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCTTCGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGtgtg                                                                                                                                                                                                                  <  2:346169/251‑1 (MQ=255)
              ggggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGATTTATCAAAGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGtgtg                                                                                                                                                                                                                  >  1:260098/1‑251 (MQ=255)
                           tAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATTCACCTGTGTTGATTTATTCACCCCCTTCTGCGCCCGTTGTGGTTTTTCCCGTTCGTTCTTATGCGGGAAAAGCCCTTCCTGCCCCACCCCCTGATGTGCAGATCCTCGTAAACAACGACGGTTTTCAGCCTTACCCTCTGGCCGGTGATG‑TGGGAAATGTCGGTATTGTTGGGGTGCGTCGCTTTGGcc                                                                                                                                                                                                     >  2:178871/1‑250 (MQ=255)
                             aGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGATGATTTACCCTTTTTTATTTTTCAACCCCCTGCTTCGCCTGTTGCCGTTTCTTCCGCTCGTTCTGTTGCGGTTAAAGCCCTTCTCGCGCAACCCCCTCCTGTTCACATCCTCGTTACCGACGACCGTTTTCAACATTTCCCTCTGGCGCGTAATG‑TGTACATTCTCGTTTTTTATCGTGTGCGTCGCTTTTGCact                                                                                                                                                                                                   >  2:233312/1‑249 (MQ=255)
                                 tGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCCGTGTTTTTTTATCAAAGCACTACTGGGCCCGTTGCGGCTTCTCCCCTTCCTTCTGATGCGGTAAAAACCATTCTGGCGCCACACCCTTATGTGCAGATCATCGTAACCGACGCCGGTTTACAGCATTACCCCCTTGCGCGTTATG‑GGGGCATTGTCGTTATTGGTGGTGCGCGTCGCTTTGGCCATGGCt                                                                                                                                                                                               >  1:336789/1‑251 (MQ=255)
                                   aaTCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGGTTTATTAAAGCAATGATGCGCCTGGTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGAAACACCCTTATGTGCAGATCAACGTAACCGACGACGGTTTACCGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGGGTGTGGCGCTTTGGCAATGGCt‑gg                                                                                                                                                                                            >  1:253405/1‑251 (MQ=255)
                                    aTCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCCACACCCTGATGTGCAGATCATCGTAACCGACGACGGGTTACCGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATTGGGTGCGTCGCTTTGGCAATGGCt‑ggg                                                                                                                                                                                           >  1:246369/1‑250 (MQ=255)
                                             gCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAACCTGTGTTGATTTATTAACGCACTGATGCGCCTGGTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGGTGTGCAGAGTATCGTAAACGCCGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGTTGGTGTGCGTCGCTTTGGCAAAGGCT‑GGGGGTGGCcgg                                                                                                                                                                                  >  2:74204/1‑251 (MQ=255)
                                                          gCAGATACCACAACAGCACAGGCGGGTGATGAACATGTGTTGGTTTATACACCCACTAAAGCGGCCGTTGCGGGTTTTCCCGTGTGTTCTTATGCGGTAAAAACCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTTACCACCGACGGATTACAGCCTTACCGTCTGGGGCGTAATG‑TAGAAATTTGCGTTTTTGATGGGGTGCGTCGCTTTGGCAATGGCTTGGTG‑TTGCCGGCgggggcgaaggg                                                                                                                                                                     >  2:294254/1‑243 (MQ=255)
                                                                aCCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTgagcg                                                                                                                                                               <  1:226733/251‑1 (MQ=255)
                                                                        aGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGgc                                                                                                                                                       <  1:193839/251‑1 (MQ=255)
                                                                            caGGCGGGTTATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACGCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCtt                                                                                                                                                   <  2:278455/251‑1 (MQ=255)
                                                                                   gTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCg                                                                                                                                            <  1:336694/251‑1 (MQ=255)
                                                                                       tgaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGAGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGGGCGGGGCGCTTAAAGGCGGTTg                                                                                                                                        >  1:366750/1‑251 (MQ=255)
                                                                                                gtTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTa                                                                                                                                <  1:233575/250‑1 (MQ=255)
                                                                                                 ttGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAAt                                                                                                                              >  1:243830/1‑251 (MQ=255)
                                                                                                       tATCAACGCACGGATGCGCCTGTTGCGGTTTCTCCCGTTCGTGCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCa                                                                                                                         <  2:211960/250‑1 (MQ=255)
                                                                                                          cAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAcgg                                                                                                                     <  1:382654/251‑1 (MQ=255)
                                                                                                          cAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAcgg                                                                                                                     >  1:199424/1‑251 (MQ=255)
                                                                                                              gCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAcggcggg                                                                                                                 >  1:516454/1‑250 (MQ=255)
                                                                                                               cACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTACAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGTGCGTAATG‑TGGCAATTGTTGTTATTGATGGAGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGTGGGAGCGGGGGGGGGGCTTTACGGCGGTTGGTGCGGTAATCGTCAGCCGCGgtg                                                                                                                >  1:195088/1‑251 (MQ=255)
                                                                                                                   gATGCGCCTGTTGCGGTTTCTCCCTTTCGTTCTGATGCGGTAAAAGCCATTCTGGCTCAATACCCTGATGTGCAGATCATCGTAACCGACGACGTTTTACAGCATTACCGTCTGGCGCGTAATG‑TGTAAAGTGTCGTTATTGATGGTGGTCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGTCGGTGTc                                                                                                              <  1:321170/249‑1 (MQ=255)
                                                                                                                   gATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTcc                                                                                                             <  1:103654/250‑1 (MQ=255)
                                                                                                                   gATGCGCCTGGTGCGGTTTCTCCCGTTTGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACCCCCTGATGTTCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTGGATGCGGTAATCGTCAACGGGGGTGTc                                                                                                              <  2:220769/249‑1 (MQ=255)
                                                                                                                     tGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTc                                                                                                          >  1:364245/1‑251 (MQ=255)
                                                                                                                      gcgtcTGTTGCGGTTTCTCCCGTTCGTTCGGATGCGGTAAAAGCCATTCTGGCGCAACCCCCTGATGTGCAGATCATCGTAGCCGACGACGGTTTACAGCATTACCGGCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCGTAAAGTCGGTTGAGGCGGTAATCGTCAACGGCGGTGTCCCTCg                                                                                                         <  2:336789/247‑1 (MQ=255)
                                                                                                                                   tttCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGGAATCGTCAACGGCGGTGTCCCTCGccaggggggaatc                                                                                            >  1:539672/1‑239 (MQ=255)
                                                                                                                                     tctccctgtcgttcTGACGCGGTAAAAGCCATTCGGGCGCAACTCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGGCTGGCGCGT‑ATGGTGGAAATTGTCGTTATTGATGGTGGGCGTCGCTTTGGCAATGGCT‑GGTTGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCATCGGCGGGGTCCCTCGCAGCGGTGAAATcc                                                                                           <  2:77155/242‑1 (MQ=255)
                                                                                                                                         ccGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGTAATCGTCAACGGCGGTGGCCCTCGCCGCGGGGAAATccccc                                                                                        >  1:234486/1‑248 (MQ=255)
                                                                                                                                              cgttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGGTGCGGTAATCGTCAAcggc                                                                                                                    >  1:147511/1‑216 (MQ=255)
                                                                                                                                              cgttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAcggc                                                                                                                    <  2:147511/216‑1 (MQ=255)
                                                                                                                                               gttcGGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTACTCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATctg                                                                                <  2:440284/251‑1 (MQ=255)
                                                                                                                                                      tGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCCATGCGTGAGCGAGCGGGGGGCTTAAAGTCGGTTGGTGCGGTAATCGTCAACGGCGGGGTCCCTCCCCGCGGGGAAATTCCCCTGCATCTgccgcccg                                                                         >  1:468908/1‑249 (MQ=255)
                                                                                                                                                                   aTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCCATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGGAATCGTCAACGGCGGGGTCCCTCGCAGCGGTGGAATTCCCATGCATCTGCTGCCCGGGCAGGCggggg                                                             >  1:447884/1‑247 (MQ=255)
                                                                                                                                                                            cAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGGGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTAc                                                   <  2:446204/251‑1 (MQ=255)
                                                                                                                                                                              acacCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCg                                                 <  2:260098/251‑1 (MQ=255)
                                                                                                                                                                                  ccTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTAc                                             <  1:74204/251‑1 (MQ=255)
                                                                                                                                                                                     gATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGGAATCCCCATGCATCTGCTGCCGGGGCAGGCGGGGAATTTAAGTACCGGGAcgcg                                          >  1:221893/1‑251 (MQ=255)
                                                                                                                                                                                          gCAGATCATCTTAATCGTCTACGGTTTACAGCATTACCTTCTGGCGCTTAATG‑TGTAAATTTTCGTTATTGATGGTGTGCTTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCACTGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACTGCGGTGTCCCTCGCATCGGTGAATTCCCCATGCATCTGCTGCCGGGTCTGGCGGTGATTTTACGTGCCGGTACTCGTtgtg                                     <  2:301750/251‑1 (MQ=255)
                                                                                                                                                                                               tcatcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGTAATCGTCAACGGCGGGGTCCCCCGCCGCGGTGGAATCCCCCTGCCTCTGCTGCCGGGGCAGGCGGGGGATTTTCCTACCGGGACGCGGTGTGGCGtt                                >  1:439619/1‑251 (MQ=255)
                                                                                                                                                                                                   gttAACCTACTACGGTTGACAGCATTTCCGGCTGTCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCTGTAATCGTCAACGGCGGTTTCCCTCGCAGCGGTGAAATCCCCATTCATCTTCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTc                            <  2:253405/249‑1 (MQ=255)
                                                                                                                                                                                                     tAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCCCGCAGCGGGGAAATCCCCATGCATCTGCTGCCGGGGCAGGGGGGGAATTTACGTAcccggaagcggtgggacggtggtccg                          >  1:537676/1‑227 (MQ=255)
                                                                                                                                                                                                      aaCCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGGCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGtt                                <  1:367604/244‑1 (MQ=255)
                                                                                                                                                                                                                     aCAGCATTACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTgg          <  1:57272/251‑1 (MQ=255)
                                                                                                                                                                                                                            tACCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATgg    <  2:364245/250‑1 (MQ=255)
                                                                                                                                                                                                                             aCCGTCTGGCGCGTAATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTCGCAGCGGGGAAATCCCCATGCATCTGCTGCCGGGGCAGGCGGTGGATTTTCGTACCGGTACGCGTTTTGACGTTGCTCCGCTTGAACATGTAGTGGCGATgggg  >  1:85635/1‑249 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
GGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATG‑TGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCT‑GGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCG  >  NC_000913/968602‑969073

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