Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,789,039 2 bp→AC coding (939‑940/960 nt) yibQ → divergent polysaccharide deacetylase domain‑containing protein YibQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,789,0390GA95.6% 130.6 / NA 45K313K (AAG→AAAyibQdivergent polysaccharide deacetylase domain‑containing protein YibQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/23);  minor base C (2/0);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
*NC_0009133,789,0400TC97.8% 136.7 / NA 45F314L (TTT→CTT) yibQdivergent polysaccharide deacetylase domain‑containing protein YibQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/23);  minor base A (1/0);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCCTC  >  NC_000913/3788798‑3789274
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                          
aaGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCt                                                                                                                                                                                                                                    >  1:88205/1‑251 (MQ=255)
        agcTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGCAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCaa                                                                                                                                                                                                                             >  1:426992/1‑250 (MQ=255)
                gaatgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGgc                                                                                                                                                                                                                    >  2:252910/1‑251 (MQ=255)
                 aatgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGgcg                                                                                                                                                                                                                   <  1:132737/251‑1 (MQ=255)
                  attaaCCGCAAGTTGTTACTTCTCCACCGCCGAATCACGCTTTGCCTGATGAACCGCGTCATTCTTTCCGTGGAGTGAAGCTGTGAAAACCGAAGACACCGATCGAACCTGTTTATGCTAATCGCTTCTTTGAATTATTAAGCGAAAGCATCAGACAGAGCACTCTGTTCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCAGAAGCGGCTAACCTTCACGCTATCGCATATTaatgacgc                                                                                                                                                                                                                  >  2:71717/4‑245 (MQ=255)
                  atgaaCCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGgcgc                                                                                                                                                                                                                  <  1:129777/251‑1 (MQ=255)
                        cGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTCCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTAcgcg                                                                                                                                                                                                            <  1:277768/251‑1 (MQ=255)
                                gATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCg                                                                                                                                                                                                     >  2:312198/1‑250 (MQ=255)
                                        tACACCACCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGACAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGACAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGCAATGCGATGGTGCg                                                                                                                                                                                             >  2:311199/1‑250 (MQ=255)
                                               ccgAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGtt                                                                                                                                                                                     >  1:354746/1‑251 (MQ=255)
                                                 gAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGACGATCTGCAAACCGAAGAAACCGACAGAAGCTGCTTATGCCATTCGCTTCTTTGCAGTATTAAGCGAAAGCATCAGCCAGAGAACGATGATCGTTTACTTACAGAATCAGTGGCAAGCCTGGGGCAAACAGCACGAAGCGGCGAAACTTCACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGtttt                                                                                                                                                                                   >  1:423749/1‑251 (MQ=255)
                                                           gtgCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCAt                                                                                                                                                                         <  1:319765/251‑1 (MQ=255)
                                                             gCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCTCGTTTTGTCGCGCCAtt                                                                                                                                                                        <  2:225001/250‑1 (MQ=255)
                                                                               tccgttccgtGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAgc                                                                                                                                                                                                                                 >  2:135350/1‑175 (MQ=255)
                                                                               tccgttccgtGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAgc                                                                                                                                                                                                                                 <  1:135350/175‑1 (MQ=255)
                                                                                           cGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTg                                                                                                                                         >  1:152727/1‑251 (MQ=255)
                                                                                            gTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTgg                                                                                                                                        >  1:255399/1‑251 (MQ=255)
                                                                                                      gCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAATGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCt                                                                                                                              >  2:290308/1‑251 (MQ=255)
                                                                                                       cAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATTGTGCGACCGCGTTTTGTCTCGCGATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCtt                                                                                                                             <  1:197092/251‑1 (MQ=255)
                                                                                                               agaaACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGc                                                                                                                                                                      <  2:244175/202‑1 (MQ=255)
                                                                                                               agaaACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGc                                                                                                                                                                      >  1:244175/1‑202 (MQ=255)
                                                                                                                      gATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTGCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACGTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACAtt                                                                                                              <  2:334818/251‑1 (MQ=255)
                                                                                                                        tAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTgg                                                                                                            >  1:202015/1‑251 (MQ=255)
                                                                                                                        tAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGACAGAGCACGCTGATCGTTTGCTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGGCACTTAACGCTAGCGCAGATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGATACCCACCTTCACGCGGCGCACATTgg                                                                                                            >  1:250811/1‑251 (MQ=255)
                                                                                                                                gTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAgg                                                                                                    <  2:154754/251‑1 (MQ=255)
                                                                                                                                          aTCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTataca                                                                                          <  2:93561/251‑1 (MQ=255)
                                                                                                                                                              aGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCg                                                                      <  2:152727/251‑1 (MQ=255)
                                                                                                                                                                        tCAGCCAGAGAACGCTGATCGTTTACTTCCAGAATCAGTGGCAGGGATGGGGCAACCAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCAATGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGAACTCAGGAGCTTATACATAACGGAAGTAAATATCTCGGCAAtcctcc                                                            >  1:334069/1‑251 (MQ=255)
                                                                                                                                                                           gCCAGAGAACTCTGATCGTTTACTTCCAGCATCAGTGGCCAGGCTGGGCCAAACAACCCGAAGCGGCGAACATTGACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGGCGCGCCATTGCCATCAGCGGACAGACCACATTACCTCCGGGCAACCCACACTCACGATGCGCACATTGGTAATCAGTCGCTTCTACCGACCGCAAGTAAATCTCTAGGCAATCATCCGtt                                                         >  2:403628/1‑251 (MQ=255)
                                                                                                                                                                            ccAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTg                                                        >  1:85418/1‑251 (MQ=255)
                                                                                                                                                                                   aCGCTGATCGTTTACTGCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAGCCCACCTTCACGCTGCGCACATTGGTAATCAGTCGCGTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTa                                                 <  1:275752/251‑1 (MQ=255)
                                                                                                                                                                                           cGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATc                                         >  2:129395/1‑251 (MQ=255)
                                                                                                                                                                                                aCTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAAGGCGATGGTGCGACCGCGTTTTGTCGCGCCAGTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAAGATCCGGc                                     <  1:371278/250‑1 (MQ=255)
                                                                                                                                                                                                                gCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATg                    <  2:306780/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGc                  <  2:151302/251‑1 (MQ=255)
                                                                                                                                                                                                                         gggCAAACAGCCCGAAGCGGCAAAACTTAACGCGAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACg           <  1:307353/251‑1 (MQ=255)
                                                                                                                                                                                                                                cAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGct    <  2:88205/251‑2 (MQ=255)
                                                                                                                                                                                                                                 aGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGctt   <  1:178189/251‑3 (MQ=255)
                                                                                                                                                                                                                                  gCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTGCACGCTGCGCACATGGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGGATGGCAAACTCGCAATGcttc  <  2:426992/251‑4 (MQ=255)
                                                                                                                                                                                                                                     cGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATCCATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATcc                                        >  1:86196/1‑210 (MQ=255)
                                                                                                                                                                                                                                     cGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATcc                                        <  2:86196/210‑1 (MQ=255)
                                                                                                                                                                                                                                     cGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCaa        >  2:123912/1‑242 (MQ=255)
                                                                                                                                                                                                                                     cGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCaa        <  1:123912/242‑1 (MQ=255)
                                                                                                                                                                                                                                      gAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGc                  >  2:268357/1‑231 (MQ=255)
                                                                                                                                                                                                                                      gAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGc                  <  1:268357/231‑1 (MQ=255)
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                          
AAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCCTC  >  NC_000913/3788798‑3789274

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