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

ATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGAT  >  NC_000913/968594‑969078
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
aTTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAAtggggg                                                                                                                                                                                                                                               >  1:129300/1‑247 (MQ=255)
   cgcgTGGGGGTCGTATCGCGGTGATATGGTGGTAAGGCTTAATCTGATCCGCTGTTATTGTCGGCTGATACCACATCAGCTCAGGCGGGTGATGAACCGTTGTTGATTTAACAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGGGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGTCGCGTAATTT‑GGatat                                                                                                                                                                                                                                           <  2:377748/251‑4 (MQ=255)
   cgcgTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTTGCGCAACGCCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGaaa                                                                                                                                                                                                                                            <  1:239718/250‑1 (MQ=255)
    gcgTGGGGGTCGTATCGCGGGGGTATGGTGGTAAGGCTGAATCGTATCCGCTGTTATTGTCGGCAGATACGACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGGGGCGGTTTCTCCCGTGCGTTCTGAGGCGGTATAATCCATTCTGGGGCAGCACGC‑TGATGTGCAGATCATGGTAACGGACGAGGGTTTGCAGCATTACCGTCTGTGGCGTAAGGT‑GGAAAtt                                                                                                                                                                                                                                          <  2:35175/251‑1 (MQ=255)
         ggggTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGt                                                                                                                                                                                                                                     <  1:16358/251‑1 (MQ=255)
                 aTCGCGGGGATATGGTGTGAAGGCTGAATCTTAGCCGCTGTTATTGTCGTCAGATACCACAACAGCACAGGCGGGTGATGAACCTGGGTTGATGTAGCAATGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACTCCA‑TGATGTGCAGATCATCGTAACCGATGCCGGTTTACAGCATTACCGGCTGGCGCGTAATGT‑GGAAATTGGCGTTATTGATg                                                                                                                                                                                                                             <  2:100968/251‑1 (MQ=255)
                 aTCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGAt                                                                                                                                                                                                                              <  2:253274/250‑1 (MQ=255)
                 aTCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGAAGAAACTGTGTTGATTTATCAACGCCCTGATGCGCCTGTTGCGGGTTCTCCCGTTTGTTCTGTTGCGGTAAAACCCATTCCGGCGCAACACCC‑TGGTGTGGAGATCATCGCAACCCACCACGGTTTACAGCCTTACCGTCATGCGCGTGAAGT‑GGGAATTGTCGTTTTTggtg                                                                                                                                                                                                                             >  1:127794/1‑248 (MQ=255)
                            atGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGAAACTGTGTTGAATTATCAACACACTGATGCGCATGGTGCGGTTTCTCCCGTTCGTTTTGGTGCGGGAAAAGCCATTCTGGCGCA‑CACCCCTGATGTGCAGATCCTCGGAACCGACGACGGTTTTAAGCATTACCGGCTGGAGCGGTATGT‑GGAAATTTGCGTTATTGTTGGGGTGGGTCGc                                                                                                                                                                                                                  >  1:343565/1‑251 (MQ=255)
                             tggtggTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTTTCAAAGCAATGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAAACGACGACGGTTTACAGCAATACCGTCTGGGGCGTAATGT‑GGAAATTGTCGTTTTTGATGGGGTGCGGCGCt                                                                                                                                                                                                                 >  1:24637/1‑251 (MQ=255)
                                tgatAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGCTGATCCTCTGTTGCGTTATCAACACACCGATGACCCTGTTGCGGTTTCTCCCGTTCGCTCTGATGCGGTAAAAACCCATCTGGCGCAAAACCC‑TGCTGTTCCGATCATCGTAACCGACGACGGGTTAACGCCTTACCCTCTGGCGCGTAATGTTGGAATTTG‑CGTTATTGTTGTTGTGCGGGGCTTTg                                                                                                                                                                                                              >  2:189181/4‑251 (MQ=255)
                                  gtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGAGAACCTGTGTTGATTTATCAACCCACTGATGCGCCTGTTGCGGTTTTTCCCGTTCGTTCTGATGCGGGACAAGCCCTTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGCCGAAGGTTTACAGCATTACCGTCTGGCGGGAAATGG‑GGTAATTGTCGTTATTGATTGTGTGCGTCGCTTTggg                                                                                                                                                                                                            >  2:417469/1‑250 (MQ=255)
                                    aaGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCa                                                                                                                                                                                                           <  1:494980/250‑1 (MQ=255)
                                                                   cAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATT‑ATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTg                                                                                                                                                                         <  1:10423/251‑1 (MQ=255)
                                                                        aCCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGGAATTGTCGGTATTGTTGGTGTGCGTCGCTTTGGCAATGGCTGGGGGTTGCCGGCGGGGGCCAA‑GCGTgaggg                                                                                                                                                                     >  1:259604/1‑249 (MQ=255)
                                                                           acaacaGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTGCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTgagcgagc                                                                                                                                                                  <  2:296024/251‑1 (MQ=255)
                                                                                    caGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGGGAGAGCGGGGGGCtt                                                                                                                                                         >  1:392906/1‑251 (MQ=255)
                                                                                        cGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCCATGCGTGAGCGGGCGGGGCGCTTaaa                                                                                                                                                      >  1:140997/1‑250 (MQ=255)
                                                                                              atgaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATTTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGtt                                                                                                                                               <  1:161998/251‑1 (MQ=255)
                                                                                                 aaCCTGTGTTGATTTATCAACGCACTGCTGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAACTTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCTGGGCGCTTAAAGTCGGTTGAt                                                                                                                                            >  1:62573/1‑251 (MQ=255)
                                                                                                      gtgtGGATTTATCAACGCACTGATGCGCCTGGTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCGGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGTATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGt                                                                                                                                       <  2:201300/251‑1 (MQ=255)
                                                                                                           gATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCg                                                                                                                                  <  1:136059/251‑1 (MQ=255)
                                                                                                                aTCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCCATGCGTGAGGGAGCGGGGCGCTTAAAGTCGGTGGATGCGGTAATCGTCAAc                                                                                                                             >  2:107978/1‑251 (MQ=255)
                                                                                                                     cGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCACCACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGCAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑GCAATGCGTGAGCGGGCGGGGCGCTTCAAGACGGTTGATGCGGTGATCGTCTAccgcgg                                                                                                                        >  2:222926/1‑251 (MQ=255)
                                                                                                                          tGATGTGCCCGTTGCGGTTTCTCCCGTTCGTTCTGATACGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACAGCCGACGGTTTACCGCATTACCGTCTGCCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGCGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGCGCGCGAGCGGGGCGCTTTAAGCCGGTTGTTGCggggactgtccaccgcggtgtcc                                                                                                                   >  2:178036/1‑230 (MQ=255)
                                                                                                                              gcgcCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTc                                                                                                               <  1:365342/250‑1 (MQ=255)
                                                                                                                                gcCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATTCGGGAGCGGGCGGGGCGCTTAAAGTCGGTTGATGGGGGAATCGTCAACGGCGGTGG‑ACCTCGc                                                                                                             >  1:119127/1‑250 (MQ=255)
                                                                                                                                     ttGCGGATTCTCCCGTTCGTTCTGATGAGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCCTCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCTATGCGTGCGCTAGCGGGGCGCTTAAAGTGGGTTGATGTGGTAATCTTGAAcgggggggtgccttgcacgggg                                                                                                        >  2:194725/1‑232 (MQ=255)
                                                                                                                                        cGGTTTCTCCCGTTCGTTCCAATGCGGCAAAAGCCATTCTGGCGAAACATAC‑TGATGTGCAGCTCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGA‑GGAACTTGTCGTTATTGCTGGTGTGCGTCGCTTTGGCACTCGCTGGTGGTTGCCGGCGGGG‑GCGATGCTGGAGCGAGCGGGGCGCTTTAAGACGGGTGATGAGGAAATAGTCCACGGAGgtgtggttccgaggcggtt                                                                                                      >  2:199366/1‑235 (MQ=255)
                                                                                                                                              ctcCCGTTCGTTCTGATGCGGAAAACGCCATTCTAGCGCAACCCCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGCCATTGTCGTTATTGATGGTGCGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGGGAGCGAGCGGGGGGCTTTAAGCCGGGTGAGGCGGGAATCGGCAACGGCGGGGTTCCCTCGcgcgcgggaacaccc                                                                                               >  2:230704/1‑237 (MQ=255)
                                                                                                                                                 ccGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGCCGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGGGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGAGGTAATCGTCAACGGCGGTGT‑CCCTCGGAGCGGTGAAATTCCCATg                                                                                           >  1:264531/1‑251 (MQ=255)
                                                                                                                                                      cgttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGTCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATct                                                                                      <  1:463952/251‑1 (MQ=255)
                                                                                                                                                      cgttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATc                                                                                       <  1:373089/250‑1 (MQ=255)
                                                                                                                                                       gttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATctg                                                                                     <  1:503016/251‑1 (MQ=255)
                                                                                                                                                            gATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATctgctg                                                                                  <  2:129300/249‑1 (MQ=255)
                                                                                                                                                                          cATTCTGGCTCAACACCC‑GGATGTGCAGATCAGCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GTAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGa                                                                  <  1:224041/251‑1 (MQ=255)
                                                                                                                                                                                   gcAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGTTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCTGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGTCGGTGAATTTTCGTa                                                         <  2:308227/251‑1 (MQ=255)
                                                                                                                                                                                       cacCC‑TGATGTGCAGGTCAGCGTAACCGACGACGGTTTGCAGGATTTCCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGATCGGGGCGCTTAAAGTCGGTGGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCTTGCATCTGCTGCCGGGTCAGGCGGTGAACTTACGTACCgg                                                     <  2:426462/251‑1 (MQ=255)
                                                                                                                                                                                             tGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCCATGCGTGAGCGAGCGGGGCGCTTTAAGTCGGTTGATGCGGTAATAGTCAACGGCGGTGT‑TCCCCCCAGAGGTGCAATCCCCCTACCTCTGCCTCCGGGTCCGGCGGGGAAttttaataccgcgacgc                                                >  2:170900/1‑237 (MQ=255)
                                                                                                                                                                                                         catcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑AGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCTGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTg                                    <  2:24637/251‑1 (MQ=255)
                                                                                                                                                                                                            cGTAACCGACGACTGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGAAGGTGTGCGGCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGTCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTGACGTACCGGTACGCGTTGTGACGTTGCTc                                 <  1:175243/251‑1 (MQ=255)
                                                                                                                                                                                                               aaCCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCCAAGCGTGAGCGAGCGGGGGGCTTAAAGTCGGTTGATGCGGGAATCGTCAACGGCGGGGT‑CCCTCGCAGCGGGGAAAACCCCCTGCGACTGCTGCCGGGGCAGGCGGGGGATTTTCGTTGCGGGACGCGTTGGGACggtgttccgc                              >  2:165494/1‑242 (MQ=255)
                                                                                                                                                                                                                     cgacgGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCCATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCTCCGCGGTGAAATCCCCATGCATCTGCCGCCCGGTCAGGCGGTGAATTTTCGTACCGGCACCCGCTTTGACGTTGCttcgctttcac                        >  1:188949/1‑241 (MQ=255)
                                                                                                                                                                                                                        cgGTTTACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGTAATCGTCAACGGCGGTGG‑CCCTCGGAGCGGGGGAAACACCCTGGATCTGCTGCCCGGGCAGGGGGGGGATTTACCACCCCGTAAGCGTTTTGACGTTGCTCCGCttttaca                       >  1:875/1‑244 (MQ=255)
                                                                                                                                                                                                                           tttACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGGGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGGCAGGCGGGGAATTTACGTACCGGGACGCGGTGTGACGTTTCCCAGCTTTGACATGTa                   >  1:321482/1‑250 (MQ=255)
                                                                                                                                                                                                                             tACAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGGGAAATCCCCATGCATCTGCTGCCGGGTCAGGGGGTGGATTTACGTAACGGGACGCGTTGTGACGTTGCTCCGCTTGTACATGTAGTg                >  1:128968/1‑251 (MQ=255)
                                                                                                                                                                                                                               cAGCATTACCGTCTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGc              <  2:259604/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTGCTTCCGGGTCAGGCGGTg                                                                   <  2:289535/186‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGGCGCGTAATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGT‑CCCTCGCAGCGGGGGAATCCCCCTGCCTCTGCTGCCGGGTCAGGCGGGGAATTTTCGTACCGGGACGCGTTGTGACGTTGCTCAGCTTGAACCTGTAGTGGCGGTGGCgggggg  >  1:33505/1‑249 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
ATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCC‑TGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGT‑GGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGG‑CCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGT‑CCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGAT  >  NC_000913/968594‑969078

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