Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,794,217 G→A A77A (GCG→GCA)  rfaD → ADP‑L‑glycero‑D‑mannoheptose 6‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,794,2170GA98.1% 172.4 / NA 52A77A (GCG→GCA) rfaDADP‑L‑glycero‑D‑mannoheptose 6‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/24);  minor base C (1/0);  total (28/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

ATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCACTACTGCCTGGAGCGTGAAATCCCGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTG  >  NC_000913/3793970‑3794455
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
aTTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTAAGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGc                                                                                                                                                                                                                                             <  2:73673/251‑1 (MQ=255)
        gTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATCAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGACGCGATTTT‑CCACGAAGGCGCATGCTCTTcccc                                                                                                                                                                                                                                     >  1:280302/1‑249 (MQ=255)
             aGTTAGGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACAGCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCg                                                                                                                                                                                                                                <  2:166820/251‑1 (MQ=255)
               ttATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGACCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATTGCGAAAGGATTTTTCCA‑GAAGGCGCATGCTCTTccaaccacgcg                                                                                                                                                                                                                              >  2:353819/1‑243 (MQ=255)
                          gTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGGGCATGCTCTTCCACCACCGAGTGGGATGGCa                                                                                                                                                                                                                    >  1:205347/1‑250 (MQ=255)
                             aCCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTTAATGATAACGGCCTCCCCGATATTCTGGGGGTGGAACACCTGAAAGAAGGCACCAAGTTTGTGAAACTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGGGATGTCGGAGAGCTTTT‑TCACCAAAGCGCCTGCCCTTCCACCCCCgcgggggatggccagt                                                                                                                                                                                                                 >  2:256821/1‑235 (MQ=255)
                             aCCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAATACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGAGGGCAAGta                                                                                                                                                                                                                <  2:248359/251‑1 (MQ=255)
                             aCCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt                                                                                                                                                                                                                 >  1:103130/1‑250 (MQ=255)
                                        gggCTTTATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGAt                                                                                                                                                                                                      >  1:14095/1‑250 (MQ=255)
                                            tttATCGGCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATaaca                                                                                                                                                                                                  >  1:313521/1‑250 (MQ=255)
                                                  ggCAGCAACATCGTTAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCGCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCa                                                                                                                                                                                           <  2:153269/251‑1 (MQ=255)
                                                        aaCATCGTTAAGGCACTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAATGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTc                                                                                                                                                                                     <  1:345626/251‑1 (MQ=255)
                                                               ttAAGGCCCTGAATGATAAAGGCATCACCGATATTCTGGCGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCaaag                                                                                                                                                                                >  1:14012/1‑247 (MQ=255)
                                                                           aTGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCactac                                                                                                                                                                  >  2:301444/1‑251 (MQ=255)
                                                                             gATAAAGGCATCTCCGATATTCTGGTGGTGAACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATAGGGATAAGGAAGACTGCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCaaatagcttccgcaatact                                                                                                                                                                 <  2:385317/250‑18 (MQ=255)
                                                                                    gCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAg                                                                                                                                                         <  1:206229/251‑1 (MQ=255)
                                                                                       tCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCACGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCgc                                                                                                                                                      >  1:234911/1‑250 (MQ=255)
                                                                                           cGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGaa                                                                                                                                                   >  2:209972/1‑250 (MQ=255)
                                                                                            gATATTCTGGTGGTGGACAACCTGAAAGACGGCAACAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAAt                                                                                                                                                 <  2:436916/251‑1 (MQ=255)
                                                                                                          ggACAACCACAAAGCCAGCACCACCTTTGTGACCCTGGTGGTTCTCGATCTCGCGGACTATATCGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATTTCGAAGCGATTTT‑CCACGAAGCCGCATGCTCTTCCACCACCGAGAGGGATGTAAAGTATATGATGGGAAACACCTATCAATACTCCAAACACCTGCTGCAATACTGTCTGGAGCACGAAATCCCATTCCTGTa                                                                                                                                     >  2:146211/1‑249 (MQ=255)
                                                                                                                       gACGGCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGcggc                                                                                                                      <  1:171138/251‑4 (MQ=255)
                                                                                                                          ggCACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCAc                                                                                                                   >  2:330023/1‑251 (MQ=255)
                                                                                                                            cACCAAGTTTGTGAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTG‑CCACGAAGGCGCATGCTCTTCCACCAGCGAGGGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCGGTTCCTGTATGCAGCGTCCGCGGCCACCt                                                                                                                 <  1:146211/251‑1 (MQ=255)
                                                                                                                                aaGTTTGTGAACCTGGTGGATATGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CAACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGAGCGAAATCCCGTTCATGTAGGCATCTTACGCGGCCACCTAcgg                                                                                                             >  1:165330/1‑251 (MQ=255)
                                                                                                                                     tgtgAACCTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGAc                                                                                                        >  1:256254/1‑251 (MQ=255)
                                                                                                                                           ccTGGTGGATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCAccc                                                                                                  >  1:71251/1‑250 (MQ=255)
                                                                                                                                              ggtggATCTGGATATCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCg                                                                                               >  1:108466/1‑251 (MQ=255)
                                                                                                                                                       ggatatCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTAt                                                                                        <  2:103130/247‑1 (MQ=255)
                                                                                                                                                       ggatatCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGGCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCGCCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTAGCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAACCCCGTTCCTGGATGCATCTTACGCGGCCACCTACGGCGGACGCATCTCCGCCtt                                                                                           <  2:385571/244‑1 (MQ=255)
                                                                                                                                                       ggatatCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAAGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATTGCAAGTATATGATGGATAACAACTATCATTACTCCAAAGAGCTATTGCACTACTGTCTGGAGCGCGAAATACCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCAGACGCACCTCCGACTTTATTg                                                                                      >  1:81556/3‑251 (MQ=255)
                                                                                                                                                          tatCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGCCGCACCCCCGGCttt                                                                                          >  2:207623/1‑244 (MQ=255)
                                                                                                                                                          tatCGCGGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACttt                                                                                          <  1:207623/244‑1 (MQ=255)
                                                                                                                                                                     aaatGGATAAGGAAGACGTCCTGATCCAGATTAGGGCTGGCGAAGAGTGCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTACACGACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGt                                                                          <  2:81556/247‑1 (MQ=255)
                                                                                                                                                                                gAAGATTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAACGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGACAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTACTGCACTACTGTCTGGAGCGCGAACTCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCCACTTTATTGAATTCCGCGAGTAAGAAAAACCGtt                                                             >  1:51677/1‑251 (MQ=255)
                                                                                                                                                                                          tGATCCAGATTATGGCTGGCGAAGAGTTCGGCGCTGTCGAAGCGCTTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTTTTTAATCCCGCGGTTACGAACAACCGTTGAACGTCTAc                                                   >  2:138722/1‑251 (MQ=255)
                                                                                                                                                                                           gATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACg                                                  <  2:384328/251‑1 (MQ=255)
                                                                                                                                                                                                 gATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGGGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGAATTTTTTTAATCCCGCGGGGACGCAAAAACGTTTAACGTCTACGGTTACt                                            >  2:239518/1‑251 (MQ=255)
                                                                                                                                                                                                    tATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGCCTTTATTGAATCCCGCGCGTACGAAAAACCGTTTAACGTCTACGGTTACTCa                                          >  2:251807/1‑250 (MQ=255)
                                                                                                                                                                                                     aTGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCaaa                                        <  1:301444/251‑1 (MQ=255)
                                                                                                                                                                                                        gctggcGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATACCAACTATAAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGTATCCCGCGAGGACGAAAAACCGATGAACGTATACGGGTACTCAAAAtt                                     >  1:282498/1‑251 (MQ=255)
                                                                                                                                                                                                          tggcGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTgg                                                                                                                                                           >  2:425097/1‑131 (MQ=255)
                                                                                                                                                                                                          tggcGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTgg                                                                                                                                                           <  1:425097/131‑1 (MQ=255)
                                                                                                                                                                                                            gcGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTg                                 <  2:238543/251‑1 (MQ=255)
                                                                                                                                                                                                            gcGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTg                                 <  1:435384/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGt               <  1:251807/250‑1 (MQ=255)
                                                                                                                                                                                                                                gAAGCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATCCTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGCCTTTATTTAATCCCGCGAGTACGAAAAACCGGTGAACGTCCCCGGTTACTCAAAATTCCTGTTTGCTGAATATGTTCGTc              >  2:385/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAAt          >  2:144672/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCGATTTT‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAAt          <  1:289460/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ttt‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGCTGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGCTGAATATGTTCGTCAAATcctgc     >  2:63012/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ttt‑CCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATcctg      <  2:14012/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ccACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTTGAGCGCGAAATCCCGTTCCTGGATGCATCTTCCGCGGTCACCTACGGCGGACGCCCCTCCGACTTTAGTGAATCCCGCGAGTACGAAAAACCGTTGAGCGTCTTCGGGTACGCAAAAGTCCTGTTTCATGAATATGTTCGTCAAATcctgcct   <  2:189585/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ccACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTATATGATGGATAACAACTATCAATACTCCAAAGAGCTGCTGCACTACTGTCTGGAGCGCGAAATCCCGTTCCTGTATGCATCTTCCGCGGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTGTGATGAATATGTTCGTCAAAGcctgcctg  <  2:306845/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
ATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAAGCCCTGAATGATAAAGGCATCACCGATATTCTGGTGGTGGACAACCTGAAAGACGGCACCAAGTTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTT‑CCACGAAGGCGCGTGCTCTTCCACCACCGAGTGGGACGGCAAGTATATGATGGATAACAACTATCAATACTCCAAGGAGCTGCTGCACTACTGCCTGGAGCGTGAAATCCCGTTCCTGTACGCTTCTTCCGCAGCCACCTACGGCGGACGCACCTCCGACTTTATTGAATCCCGCGAGTACGAAAAACCGTTGAACGTCTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTG  >  NC_000913/3793970‑3794455

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