Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,980,726 G→A 15.4% intergenic (+30/‑161) wecG → / → yifK UDP‑N‑acetyl‑D‑mannosaminuronic acid transferase/putative APC family amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,7260GA15.4% 139.3 / 17.6 59intergenic (+30/‑161)wecG/yifKUDP‑N‑acetyl‑D‑mannosaminuronic acid transferase/putative APC family amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base A (5/4);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-03

GGATCGCCAAAGCAGGAGATCATCATGCGCGACTGCCGTCTGGTACATCCAGATGCGCTGTATATGGGCGTTGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGC  >  NC_000913/3980460‑3980987
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                     
ggATCGCCAAAGCAGGAGATCATCATGCGCGACTGCCGTCTGGTACATCCAGATGCGCTGTATATGGGCGTTGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCg                                                                                                                                                                                                                                                               >  4:491874/1‑275 (MQ=255)
                   tcatccTGCGCGCCGTCCCTCTGGGACCACCTGATGCGCTGTATATGGGCGTTGGCCGGACTTACCATGTTTTCACCGGTCACGTACAACGCGCACCGAAAATCTGGCAAACGCTGTGGCTGGAGTTGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTtta                                                                                                                                                                                                                                            <  4:279201/275‑1 (MQ=255)
                       catGCGGGACTGCCGTCTGGTACATCCAGATGCGCTGTATAGGGGCGTTGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTac                                                                                                                                                                                                                                        <  2:300954/275‑1 (MQ=255)
                                                    aTGCGCTGTATATGGGCGTTGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTTTCTTTGTAAAGCGACGGAGTGGGCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAt                                                                                                                                                                                                           >  1:78619/1‑275 (MQ=255)
                                                                      ttGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAGCGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGaa                                                                                                                                                                                          >  2:41594/1‑274 (MQ=255)
                                                                      ttGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTTTTTACGCTGGCACTACACCGGCCACCTATGATTTTCCTCTCATTGTAAAGCGACGGAGTGGTCACCCCGTCGCCTGCCTGCTTTTTTATTACACAAAGCATTAAAATTTTTAATGCTTTATTTGCCATTTCTCCTGat                                                                                                                                                                                          >  2:535923/1‑273 (MQ=255)
                                                                      ttGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGTACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGaa                                                                                                                                                                                          >  2:197366/1‑274 (MQ=255)
                                                                      ttGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCCCTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGACGCCTGTCTGCTTTTTTTTTCCACAAAGCATTCAAATTTTTAATGCTTTTTTTCCCATTTCTCCTaaa                                                                                                                                                                                          >  4:371390/1‑274 (MQ=255)
                                                                      ttGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCCACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACCCAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCCCCTGaa                                                                                                                                                                                          >  2:428353/1‑274 (MQ=255)
                                                                      ttGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGaa                                                                                                                                                                                          >  2:447829/1‑274 (MQ=255)
                                                                      ttGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTNACGCTGGCACTACACCGGCACCCTATGATTTTCCTCTCTTTGTAAAGCGCCGGAGTGGTCACTCCGTCGCCTGTNTCCTTTTTTCTTACACAAAGCATTCAAATTNNTAATTCTTTTTTTTACATTTTCCCTTAAt                                                                                                                                                                                         >  3:61621/1‑275 (MQ=255)
                                                                          cGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTc                                                                                                                                                                                                                                                                     <  4:558969/195‑1 (MQ=255)
                                                                                         tttCACGGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAaca                                                                                                                                                                      <  1:61621/275‑1 (MQ=255)
                                                                                         tttCACCGGTTACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAaca                                                                                                                                                                      <  1:104299/275‑1 (MQ=255)
                                                                                           tctccGGTCACTCAACAGGCGCACCGANNTTCGTGCGACCTCTGGGACTTGAGTTGCTTNCACTCCTGTCTTCGCAGCCGAGCCGCATTAAGCGTCAGTTTCGTGTGCTGGGTTGTTGACGATAGGACTACACCNGGAACGTATAATTTTACACTCTTTGTAAGACGACGGAGTGCTCACTCCGTAGCGTGTCTGCTTGATTATTAATCAAAGCAATCAAATATTTTATGCTTTATTTGCCATGTCTCATGAATTACGAAAACAttttataca                                                                                                                                                                      <  3:68552/270‑7 (MQ=255)
                                                                                               cgacagcGTAAAACGCGCACCGAAAATCGGGGAAACGCTGGGGCTGGGGTGGCCCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTAGTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGAt                                                                                                                                                                <  4:484647/269‑1 (MQ=255)
                                                                                                ggTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCGCTACACNGGCACGCTATGATTTTCCTATCTTTGTACAGCGCCGGGGTGATGACTCCGTCGCCTGTCTGCTTTTTTTNTACACCATGCCTTCAAATTTATAAGGCTTAATTTGCCAGTTTTCCNTNATTTATACAACATTTGCACCAATcaagt                                                                                                                                                               >  3:53556/1‑271 (MQ=255)
                                                                                                ggTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTTCCATTTCTCCTGAATTACGAAAACTTTTTCAACACTCGATg                                                                                                                                                               >  4:494187/1‑275 (MQ=255)
                                                                                                ggTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATg                                                                                                                                                               <  2:491875/275‑1 (MQ=255)
                                                                                                  tCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACAt                                                                                                                                                                              >  2:5812/1‑258 (MQ=255)
                                                                                                      gTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGc                                                                                                                                                                          <  4:5812/258‑1 (MQ=255)
                                                                                                       tAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGATTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGTACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCAt                                                                                                                                                        <  4:197365/275‑1 (MQ=255)
                                                                                                                cACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACCTTTGCAACACt                                                                                                                                                                    >  2:484648/1‑254 (MQ=255)
                                                                                                                           tGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCg                                                                                                                                                                                                                                                               >  4:175060/1‑152 (MQ=255)
                                                                                                                             gCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGa                                                                                                                                  >  2:111134/1‑275 (MQ=255)
                                                                                                                               aaaCGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGc                                                                                                                                                                                                                                                           <  2:175061/152‑1 (MQ=255)
                                                                                                                                 aCGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTgg                                                                                                                                                                                                                                                                       <  6:239540/138‑1 (MQ=255)
                                                                                                                                 aCGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTgg                                                                                                                                                                                                                                                                       >  5:239540/1‑138 (MQ=255)
                                                                                                                                  cGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTAAGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCa                                                                                                                             >  2:106/1‑275 (MQ=255)
                                                                                                                                      ggggCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCa                                                                                                                         <  2:37120/275‑1 (MQ=255)
                                                                                                                                             gatttgtcctACCGCCTGCTTTCGCGGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacacaa                                                                                                                  <  4:359603/267‑1 (MQ=255)
                                                                                                                                                       aCCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATCACCGGTAACAACGGAAAGCATGCAAACACAAAACGAGGAtt                                                                                                        >  4:147712/1‑275 (MQ=255)
                                                                                                                                                             tGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGc                                                                                                  >  2:537211/1‑275 (MQ=255)
                                                                                                                                                                 ttCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCa                                                                                                                                                                                                                                                                    >  5:134260/1‑109 (MQ=255)
                                                                                                                                                                 ttCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCa                                                                                                                                                                                                                                                                    <  6:134260/109‑1 (MQ=255)
                                                                                                                                                                        ccgagccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCGt                                                                                                                                                                                                                                                              >  6:968328/1‑108 (MQ=255)
                                                                                                                                                                        ccgagccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCGt                                                                                                                                                                                                                                                              <  5:968328/108‑1 (MQ=255)
                                                                                                                                                                          gagccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCa                                                                                     <  2:147713/275‑1 (MQ=255)
                                                                                                                                                                                     cgtgctcccTTCGTTTGGTGCGTTTTTTACTCTGGCACTACACCGGCGACCTATGTTTTTCCTCTCTTTGTAAAGCGACGGAGTGGGCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAGTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCg                                                                          <  4:537210/267‑1 (MQ=255)
                                                                                                                                                                                           aGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACCCGATGTACCCATAACGATAACCGGTACCACCGGAAAGCATGCAAACACAACACGAGGATTTATGGAAGATAACAAACCAGAGCAACAGCGTGGACt                                                                    >  2:494226/1‑275 (MQ=255)
                                                                                                                                                                                            gCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCt                                                                    <  2:494188/274‑1 (MQ=255)
                                                                                                                                                                                               tCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCGt                                                                                                                                                                                                                                                              >  5:253022/1‑85 (MQ=255)
                                                                                                                                                                                               tCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCGt                                                                                                                                                                                                                                                              <  6:253022/85‑1 (MQ=255)
                                                                                                                                                                                                cGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAg                                                                                                >  1:32073/1‑242 (MQ=255)
                                                                                                                                                                                                   ttGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCg                                                                                                                                                                                                                                                               <  7:466425/80‑1 (MQ=255)
                                                                                                                                                                                                   ttGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCg                                                                                                                                                                                                                                                               >  8:466425/1‑80 (MQ=255)
                                                                                                                                                                                                    ttttgcGTTATTTACGCTGGCACTACACCTGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTTGNCACTCCGTCGCCTGTCTGCTTTTTTATTACANAAAGCATTNNAATTTTTAATGCTTTATTTGCCANTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATaa                                                                                            <  3:32073/239‑1 (MQ=255)
                                                                                                                                                                                                       tgcGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGAc                                                         <  4:41593/274‑1 (MQ=255)
                                                                                                                                                                                                           ttATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTgg                                                                                                                                                                                                                                                                       <  7:565511/64‑1 (MQ=255)
                                                                                                                                                                                                           ttATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTgg                                                                                                                                                                                                                                                                       >  8:565511/1‑64 (MQ=255)
                                                                                                                                                                                                                tACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTttatta                                                                                                                                                                                                                                         >  3:106803/1‑89 (MQ=255)
                                                                                                                                                                                                                  ggctggaattaacccGGCAACCTATGATTTTGCTCTCCTTGTAAACCGACGGGGTGGTCATCCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTGGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCa                                             <  4:428352/262‑1 (MQ=255)
                                                                                                                                                                                                                  cGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCa                                             <  1:96306/275‑1 (MQ=255)
                                                                                                                                                                                                                  cGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCg                                                                                                                                                                                                                                                               <  6:830365/65‑1 (MQ=39)
                                                                                                                                                                                                                  cGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCg                                                                                                                                                                                                                                                               >  5:830365/1‑65 (MQ=39)
                                                                                                                                                                                                                    cTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTAcaca                                                                                                                                                                                                                                     <  1:106803/89‑1 (MQ=255)
                                                                                                                                                                                                                               ccGGCAACCGTTGGTTTTCCTTTCTTTGTAAAGCGACGGATTGGTCCCTCCGTCGCCTGTCTGCCTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                <  4:447828/275‑1 (MQ=255)
                                                                                                                                                                                                                                               tttctctctTTTTGAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACCCAAAGCATTCAAATTTTTCATGCTGTATTTGCCATTTATCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGc                <  4:106/272‑1 (MQ=255)
                                                                                                                                                                                                                                                tCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGcc               >  2:241041/1‑275 (MQ=255)
                                                                                                                                                                                                                                                         tGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATg       >  4:390216/1‑274 (MQ=255)
                                                                                                                                                                                                                                                           gaaatccaCGTAGTGGGTACTCCGTCGCCTCTCTGCTTTGTTATTACACAAAGCATTCAAATGTTTAATGCTTTATTTGCCAGTTCTCCTGAAGTACGAAAACATTTGCAACACTCGATGTACCCATAACGATACCCGGTAACACCGGAAATCATGCAAACACAACTCGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCCTTGGCGTCGGCCTGTTTATgggg    <  1:53556/268‑1 (MQ=255)
                                                                                                                                                                                                                                                             aaGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCAGTCAAATTTTTAATGCTTTATTTGCCATTTCTACTGAATTACGAAAACATTTGCAACAATCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc  <  2:371391/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATg                                                                                                                           >  5:174642/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                     
GGATCGCCAAAGCAGGAGATCATCATGCGCGACTGCCGTCTGGTACATCCAGATGCGCTGTATATGGGCGTTGGCGGGACTTACGATGTTTTCACCGGTCACGTAAAACGCGCACCGAAAATCTGGCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGC  >  NC_000913/3980460‑3980987

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