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,8370GA97.1% 195.5 / NA 69D158N (GAT→AAT) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/49);  minor base C (2/0);  total (20/49)
Fisher's exact test for biased strand distribution p-value = 8.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTGGGC  >  NC_000913/968593‑969083
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
gaTTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGTCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGAGTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAAtgtg                                                                                                                                                                                                                                                    <  1:805083/249‑1 (MQ=255)
        tGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGGTGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAACCATTTTGGAGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGTc                                                                                                                                                                                                                                          >  1:68059/1‑251 (MQ=255)
             ggcGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGAGTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACTACGGTTTACAGCATTACCGGCTGGCGCGTAATGTGGAAATTGTCGTTAt                                                                                                                                                                                                                                     <  2:49996/249‑1 (MQ=255)
                gTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGAGATGAACCTGTGTTGATTTATCAAACCACTGATGCGCCTGTTGCGGTTTCTCCCCTTCGTTCTGATGCGGTAAAAGCCATTCTTGCGCAACACCCTGCTTTGTAGATCATCGTAACCGGAGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGTCGTTATTgt                                                                                                                                                                                                                                  >  1:842961/1‑250 (MQ=255)
                       ggggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGtgtg                                                                                                                                                                                                                           <  1:772604/251‑1 (MQ=255)
                             atGGTGGTAAGGCTGAATCTTTTCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCg                                                                                                                                                                                                                      <  1:809545/250‑1 (MQ=255)
                               ggtggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCtt                                                                                                                                                                                                                   <  1:589365/251‑1 (MQ=255)
                                          tGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAAACTTTGTTGATTTATCAACCCAATGATGGGCCTGTTGCGGGTTTTCCCGTTTGGTCTGACGCGGGGAAAACCATTCTGGGGCAACAACCTGGTGTGCAGATCATCGTTACCCGAGGCCGGTTACAGCAATACCGTCTGGCGCGTAATGTGGGAATTGTCGGTATTTGTGGTGTGGGTCGCTTTTGCAATgggt                                                                                                                                                                                                        >  1:419421/1‑249 (MQ=255)
                                           gAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGGTGATTTATCAAAGCCATGATGAGCCTGTTGCGGGTTCTCCCGTTCGTTCTGGTGCGGGAAAAGCCATTCTGGCGCCACACCCCAATTTGCAGATCATCGTAACCGAAGACGGGTTACAGCATTAACGTCTGGCGCGTAATTTGGGAAATTTAGTTATTGATGGGGTGCGGCGCTTTGGGAATGGCtg                                                                                                                                                                                                       >  1:798370/1‑251 (MQ=255)
                                            aaTCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAACCTGTGTTGATTTATCCACGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGGAATGTGGGAATTGTCGTTATTGATTGTGTGCGTCGCTTTGGCAATGGCtgg                                                                                                                                                                                                      >  1:266784/1‑251 (MQ=255)
                                                   tCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGGACCTGTGTTGATTTATCCACCCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAACCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTACCGGCTGGGGCGTAATTTGGAAATTGTTGTTATTGGTGGTGTGCGTCGCTTTGGCAATGGGTGGTGGTTgt                                                                                                                                                                                               >  1:620027/1‑250 (MQ=255)
                                                   tCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATAAAAGCACTGATGCGCCTGTTGCGGTTTCTCCCGGTCGTTCTGATGCGGGAAAAGCCATTCTTGCGCAACAACCTGATGTGCAGATCATCGTAACCGCCGACGGTTTACAGCATTAACGTCTGGCGGGTAATGTGGAACTTGTCGGTTTTTATGGTGTGCGTCGCTTTTGCAATAGCTGGGGGTTGc                                                                                                                                                                                               >  1:480712/1‑251 (MQ=255)
                                                    ccGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGc                                                                                                                                                                                               <  1:723271/250‑1 (MQ=255)
                                                    ccGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGc                                                                                                                                                                                               <  1:179016/250‑1 (MQ=255)
                                                    ccGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGAGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCAGTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGcc                                                                                                                                                                                              <  1:831882/251‑1 (MQ=255)
                                                            aTTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCgggg                                                                                                                                                                                       <  2:180102/250‑1 (MQ=255)
                                                               gTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCaa                                                                                                                                                                                   <  2:649964/251‑1 (MQ=255)
                                                                  ggCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAACCCTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTACCGTCTGGCGCGTCATGTGGAAATTGTCGTTATTTTTGGTGTGCGTCGCTTTGGCAATGGCTGGGGGTTGCCGGCgggggcgatgg                                                                                                                                                                                >  1:229350/1‑244 (MQ=255)
                                                                                    acaGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCt                                                                                                                                                              <  1:166461/251‑1 (MQ=255)
                                                                                        gCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACCCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTaaa                                                                                                                                                          <  2:266784/251‑1 (MQ=255)
                                                                                           ggTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTc                                                                                                                                                       <  1:523941/251‑1 (MQ=255)
                                                                                            gTGATGAACCTGTGTTGACTTATCAACGCACTGATGCGCCTGTTGAGGTTTCTCCCGTTCGTTCTGATGCGGTAAACGCCATTCTGGCGCAAAACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTCATGTGGAAATTGTCGTTATTGATGGAGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGACCGAAGCGTGAGCAAGAGGGGCGGTTagcatcg                                                                                                                                                      >  2:407074/1‑245 (MQ=255)
                                                                                                  aaCCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGAt                                                                                                                                                <  1:548376/251‑1 (MQ=255)
                                                                                                   atcgttgtggATTTAGCAACGTACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATTCGGTAAAAGCCATTCTGGCGCAACACCCGGATGTGCAGATCATCGTATCCGACGACGTTTTACAGCATGACCGTCTGGTGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTGTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAATTCGGTTGATg                                                                                                                                               <  2:803163/242‑1 (MQ=255)
                                                                                                         gtTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCGTAAAGTAGGGTGATGCggac                                                                                                                                          >  1:180059/1‑248 (MQ=255)
                                                                                                            gATTTATCAACGGACTGATGCGCCTGTGGCGGTTTCTCCCGTTCGTGGTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCg                                                                                                                                      <  2:344295/251‑1 (MQ=255)
                                                                                                            gATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCa                                                                                                                                                                                                              <  1:371212/179‑1 (MQ=255)
                                                                                                            gATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCa                                                                                                                                                                                                              >  2:371212/1‑179 (MQ=255)
                                                                                                             ttttATCAACGCATTGATGGGCCTGTTGCGGTTGCTCCCGTTCGGTCTGTTGCGGTAAAAGCCATGCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAGTCg                                                                                                                                      <  2:44506/249‑1 (MQ=255)
                                                                                                                tATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCACTGCGTGAGAGAGAGGGGGGCTTAAAGTCGGTTGTTGCGGTAAACGTCa                                                                                                                                   >  1:833986/1‑250 (MQ=255)
                                                                                                                tATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGGGCGGGGGGCTTTAAGTCGGTTGTTGCGGTAATCGTCa                                                                                                                                   >  1:672089/1‑250 (MQ=255)
                                                                                                                      agcaatGATGCGCCTGTGGTGGTTTATCACGATCGTGCTGATGCGGTAAAAGTCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTATAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAAGCGTCAAcggcgg                                                                                                                            <  1:605617/246‑1 (MQ=255)
                                                                                                                       ggaCTGATGCGCTTGTTGCGGTTTCTCCCGTTCGTTCTGATGCTGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGgt                                                                                                                           <  2:478831/249‑1 (MQ=255)
                                                                                                                        cACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGGCAATGCGTGAGCGCGCGGGGGGCTTAAAGGCGGGTGGTGCGGGAATCGGCAAcggcgggg                                                                                                                          >  1:36447/1‑249 (MQ=255)
                                                                                                                                   cTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATACGTGAGCGAGCGGGGGGCTTAAAGGCGGGTGATGCGGGAATCGTCAAaggcgggggccccccgcgc                                                                                                               >  1:605212/1‑238 (MQ=255)
                                                                                                                                     gTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGc                                                                                                               <  1:459467/249‑1 (MQ=255)
                                                                                                                                        gCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATGCGGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGGCGTTATTGATGGTGTGCGTCGCTGTGGCAATTGCTGGTGGTTGCCGTCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTg                                                                                                           <  1:51663/250‑1 (MQ=255)
                                                                                                                                          ggTTTCTCCCGTTCGTTCGGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGaaa                                                                                                        <  2:237611/251‑1 (MQ=255)
                                                                                                                                               cgcccgtttgatcggATGCGGTAAATGTCAGTCTGGCCCAACACCCTGATGTGCCGATCTTCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAAGGTGGAAATTGTCGTTATTTGTGGTGTGCTTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATcccc                                                                                                   <  2:625947/237‑1 (MQ=255)
                                                                                                                                                 cccGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTCGCCGCgggggaaacccccc                                                                                                 >  1:287754/1‑239 (MQ=255)
                                                                                                                                                 cccGTTCGTTCTGATGCGGTAAAAGCCATTCGTGCGCAACACTCTGATGGTCAGATAATCGGATCCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGGCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTTGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCGTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCAt                                                                                                 <  2:283633/251‑1 (MQ=255)
                                                                                                                                                  ccgtccgtTCTGATGCTGTAAATGCCATTCTGGCGCAATCCCCTGATGTGCAGATCACCGTAACCGACGACGGTTTACAGCATTACCGTCTGTCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGACCCTCGCAGCGTTGAAATCCCCATg                                                                                                <  2:458368/251‑1 (MQ=255)
                                                                                                                                                      gtgtgttggtGCGGTATAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCAGTACCGTCTGGCGCGTAATGTGGAAAGTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGTCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGGCGGTGGATGCGGTAATCGTCAACGGCGGTGTCCCTCTCAGCGGTGAAATCCCCATGCa                                                                                              <  2:709994/240‑1 (MQ=255)
                                                                                                                                                             gATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATctgct                                                                                        <  1:154047/248‑1 (MQ=255)
                                                                                                                                                                gCGGTAAAAGCCATTCTGGCGCTACTCCCTGATGTGCAGAGCATCGTAACCGACGTCGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGGCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTTGTGGTTGCCGGCGTGGCCAATGCGTGAGCGAGCGGGTCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCGCGCATCGGTGAAATCCCCATGCATCTGCTGCCggg                                                                                  <  2:217906/251‑1 (MQ=255)
                                                                                                                                                                      aaaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGc                                                                            <  1:310506/251‑1 (MQ=255)
                                                                                                                                                                       aaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCAt                                                                                                 <  2:824839/229‑1 (MQ=255)
                                                                                                                                                                       aaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCg                                                                           <  1:655108/251‑1 (MQ=255)
                                                                                                                                                                       aaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGGGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGTCGGGGCCCAAGACTGAACGAGAGGGGCGCTTAAAGTCCGTTGTTGCGGTAAGCATCAACGGCAGTGTCCCCCACAGCAGGGAAATCCCCactccacctcgcacgggtacagcg                                                                           >  1:48144/1‑228 (MQ=255)
                                                                                                                                                                       aaGCCATTCGGGCGCAACACCCTGATGTGCATATCATCGTAGCCGACGACGGTTTACAGCATTACCGGCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTTCCGGCGGGGCCTATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCg                                                                           <  1:496055/251‑1 (MQ=255)
                                                                                                                                                                            aTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGaa                                                                      <  1:498717/251‑1 (MQ=255)
                                                                                                                                                                                      aaCACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGGAATCGTCAACGGCGGGGTCCCCCGCAGCGGTGGAATTCCCCTGCATTTGCTGCCGGGGCCGGCGCGGAATTTatggaac                                                            >  1:160799/1‑245 (MQ=255)
                                                                                                                                                                                       acacCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCg                                                           <  1:446522/251‑1 (MQ=255)
                                                                                                                                                                                           ccttatttgttgctcatcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCGGGTTGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAAGCGTCATCGGCGGTGTCCCCCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTAc                                                       <  2:842961/238‑1 (MQ=255)
                                                                                                                                                                                              gATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTAcg                                                      <  2:672089/249‑1 (MQ=255)
                                                                                                                                                                                                     aGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGGGAGGGCGGGGGGCTTAAAGTCGGTTGATGTGGTAATCGTCCACGGCGGTGTCCCTCGCAGCGGGGAAAACCCCCTGCCTATGCTTCCGGGTCAGGCGGTGAAttttcatacccggacgcggtttg                                               >  1:812378/1‑229 (MQ=255)
                                                                                                                                                                                                       agtatcGTATCCGGCGTCGGTTTACAGCAGTACCGGCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTTGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCTATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGt                                           <  1:816747/248‑1 (MQ=255)
                                                                                                                                                                                                         catcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGGGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGtt                                          <  2:783723/250‑1 (MQ=255)
                                                                                                                                                                                                                  ccGACGACGGTTGACAGCATTACCGTCTGGCGCTTAATGTGGAAATTGTCGTTATTGATTGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGTCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGGAAGCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCTGTACGCGTTGTGACGTTGCTCAGCtt                                 <  2:798370/251‑1 (MQ=255)
                                                                                                                                                                                                                        ccgttgtACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCTGTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACTGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGTCGTTGAATTTACGTACCGGTACGCGTTGTGTCGTTGCTCAGCTTGAACAt                           <  2:375647/245‑1 (MQ=255)
                                                                                                                                                                                                                         cgGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATg                          <  1:839517/251‑1 (MQ=255)
                                                                                                                                                                                                                              tACAGCATTTCCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGt                      <  1:492268/250‑1 (MQ=255)
                                                                                                                                                                                                                                 aGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCg                  <  1:37429/251‑1 (MQ=255)
                                                                                                                                                                                                                                   catgaCCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGGCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGAt                <  2:68059/247‑1 (MQ=255)
                                                                                                                                                                                                                                     ttACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGGCCCTCCCCGCGGGGAAATCCCCATGCATCTGCTGGCGGGGCAGGGGGGGGAATTACGTACCCGGACGCGTTGTGGAGTTGCCCAGccttaaacttgagaggcggtgg              >  1:247948/1‑230 (MQ=255)
                                                                                                                                                                                                                                       aCCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGc             <  1:750410/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCgggg         <  1:808480/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ggCGCGTAATGTGGAAATTTGCTGTATTGATGGTGTTCGGCGCTTTGGCAATGGCTGGTGGTTGCCCGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGTTTGATGCTGTAATCGTCAACGTCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTTAATTTTCGTGCCGTTACGCGTTGTTACGTTGCTCAGCTTGATCATGTAGTGGCGATGGCGGGGAtt      <  2:847177/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 gcgTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTGGGc  <  1:824591/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
TATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTGGGC  >  NC_000913/968593‑969083

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