Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,985,726 A→G intergenic (‑134/‑121) yqeF ← / → yqeG putative acyltransferase/putative transporter YqeG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,985,7260AG98.2% 194.4 / NA 55intergenic (‑134/‑121)yqeF/yqeGputative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (30/24);  total (31/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 = 9.97e-01

GCTTTCACGACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAAGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTG  >  NC_000913/2985486‑2985960
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           
gCTTTCACGACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAAc                                                                                                                                                                                                                                  >  2:262063/1‑251 (MQ=255)
     gaCGACCAGCATAACAAGTTCCACGGCGGAAAGACCCGCTAACGCACCACGAAAGCAGCCGATAGGGGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACt                                                                                                                                                                                                                              >  2:157670/2‑250 (MQ=255)
     cACGACCAGACTACCAAGTTCCACGTCCGAATGACCCTCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGTT‑TTGCAACCACt                                                                                                                                                                                                                              <  2:630196/250‑1 (MQ=255)
        gACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAACAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAAcc                                                                                                                                                                                                                          >  1:522992/1‑251 (MQ=255)
              aCTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACccacca                                                                                                                                                                                                                     >  2:479407/1‑250 (MQ=255)
                   cAAGTTCCACGGCTGAATGACCCGCTAACGCACCCCGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAAGTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATagag                                                                                                                                                                                                               <  2:386522/251‑1 (MQ=255)
                                ggAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGc                                                                                                                                                                                                  >  1:21853/1‑251 (MQ=255)
                                 gAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACACTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTAAGAGACATAACTGGCAAGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCg                                                                                                                                                                                                 >  2:551275/1‑251 (MQ=255)
                                             aaCGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCAt                                                                                                                                                                                     <  1:465603/251‑1 (MQ=255)
                                                                 gATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAATCTGGTTTATGGTa                                                                                                                                                                 <  1:569972/251‑1 (MQ=255)
                                                                 gATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTa                                                                                                                                                                 <  1:305273/251‑1 (MQ=255)
                                                                      gtgtCCGTAACGCCCCGACAATCACAACGTCTTTCATCGCAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTAGCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGTGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTg                                                                                                                                                            <  2:606249/251‑1 (MQ=255)
                                                                                 gCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCAttct                                                                                                                                                 <  1:572026/251‑3 (MQ=255)
                                                                                   cccGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGt                                                                                                                                               >  1:535782/1‑251 (MQ=255)
                                                                                         aaTCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAAtt                                                                                                                                         <  1:460161/251‑1 (MQ=255)
                                                                                               aaCGTATTTCATCACAACCTCGCACCAATTAGCAGTACAATAGTAAAGGATTGTTAATTTCAGTTATCTTATATGTTTAAAAGAAGTGCTTTTTATCTCAAAGGAAATATGCCTCAGACA‑CGGTCATAGACCTAGCTGGCACGTAGGGT‑TTGCAACCACAAACCCACCAATGGAGGGGTAGATATGGCGTTAATCTCCCATACTTAACCTGGTTTTGGGGAAATTTCCCGCCAttctgcgtattttgtagc                                                                                                                                   <  2:574728/251‑17 (MQ=255)
                                                                                                      ttCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAgatgat                                                                                                                               <  2:194048/248‑1 (MQ=255)
                                                                                                        catcaCAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGtt                                                                                                                          >  1:115114/1‑251 (MQ=255)
                                                                                                             caACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaa                                                                                                                     <  1:479407‑M1/251‑5 (MQ=255)
                                                                                                                       ccAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccga                                                                                                           <  2:115114‑M1/251‑15 (MQ=255)
                                                                                                                         aaTTAACAGAGCACTCCTAAATGATTGTTCATTACAGTTATCTTAACTGTTTAAAAAAAGCGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTCACCCACCAATAGAGGTGGAGATAGGGCGTTAATCTCCCATCCTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgac                                                                                                          >  2:286651‑M1/1‑235 (MQ=255)
                                                                                                                                 gTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccgg                                                                                                 <  1:498839‑M1/251‑25 (MQ=255)
                                                                                                                                   gCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTGTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccgggg                                                                                               <  2:21895‑M1/251‑27 (MQ=255)
                                                                                                                                        aGTAAATTATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAGATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtg                                                                                          >  2:587347‑M1/1‑220 (MQ=255)
                                                                                                                                        aGTAAATGATTATTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACACAGGAACTATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGCAGGGT‑TTGCAACCACAAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaacaaccccgacatttgccggggttgtg                                                                                          >  2:102130‑M1/1‑220 (MQ=255)
                                                                                                                                              tGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcat                                                                                     >  1:147803‑M1/1‑214 (MQ=255)
                                                                                                                                               gATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcata                                                                                    >  1:598168‑M1/1‑213 (MQ=255)
                                                                                                                                                aTTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataac                                                                                  >  1:500893‑M1/1‑212 (MQ=255)
                                                                                                                                                    ttAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGCG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgta                                                                               >  1:595798‑M1/1‑208 (MQ=255)
                                                                                                                                                            aGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttatt                                                                      >  1:17187‑M1/1‑200 (MQ=255)
                                                                                                                                                             gTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattt                                                                     <  1:466261‑M1/251‑53 (MQ=255)
                                                                                                                                                                   ttAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccgg                                                               >  1:253805‑M1/1‑193 (MQ=255)
                                                                                                                                                                      aTTGTTTAAAAAAAGCGATTTTTATCACAAAGGAAATATGACTGAG‑TAACAGTCAGAACCATAACTGCCACGTAGGGA‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGCTGATACGTGCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgattattttaccggacg                                                            >  2:85031‑M1/1‑190 (MQ=255)
                                                                                                                                                                               aaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGAGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtcaaaaaaccccgacatttgccggggttgtgatcataacgtaatgcttattttaccggtcgcatctccgg                                                   >  2:3199‑M1/1‑181 (MQ=255)
                                                                                                                                                                                  aaGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaa                                                <  1:262063‑M1/251‑74 (MQ=255)
                                                                                                                                                                                            tATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacataataac                                       <  1:608301‑M1/250‑83 (MQ=255)
                                                                                                                                                                                             aTCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacg                                      <  1:615842‑M1/250‑84 (MQ=255)
                                                                                                                                                                                              tCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtc                                    >  2:13524‑M1/1‑166 (MQ=255)
                                                                                                                                                                                               cacaAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcg                                   >  1:284457‑M1/1‑165 (MQ=255)
                                                                                                                                                                                                acaAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgc                                  <  2:21853‑M1/251‑88 (MQ=255)
                                                                                                                                                                                                      gcaataTGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatgg                            >  1:321691‑M1/3‑158 (MQ=255)
                                                                                                                                                                                                      gAAATCTGCCTGAG‑CAGACGTCAGAGACCTACCTGGCACGTAGGGT‑ATGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTACCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatgg                            >  2:511005‑M1/1‑158 (MQ=255)
                                                                                                                                                                                                        aaTATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggta                          <  1:103143‑M1/251‑96 (MQ=255)
                                                                                                                                                                                                        aaTATGCATGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGA‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGAGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggt                           >  1:553629‑M1/1‑156 (MQ=255)
                                                                                                                                                                                                             gCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgctattttaccggacccatccccgggaacagaataacgtcgcggatggtatggccg                     >  1:601289‑M1/1‑151 (MQ=255)
                                                                                                                                                                                                               cTGAG‑CAGCAGTCAGAAACATAACTGGCCCGT‑AGATGTTGCAACCACTAACCCAACGATAGAGGGGAAGAGAGGGCGTTAATCTCCCATACTTAAACAGGTTTATGGTAAATTGCCCTCCATTCTGATTAATTTGTAGCTGATACGTTCtaaaaaaacccctacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggaaggtatggctg                    >  2:239020‑M1/1‑149 (MQ=255)
                                                                                                                                                                                                                  ag‑cagcAGTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttgg                <  2:522992‑M1/251‑106 (MQ=255)
                                                                                                                                                                                                                         agTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaaca          <  1:60970‑M1/251‑112 (MQ=255)
                                                                                                                                                                                                                         agTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCT‑CCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaac           >  1:359123‑M1/1‑139 (MQ=255)
                                                                                                                                                                                                                          gTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaaca          <  2:117574‑M1/250‑112 (MQ=255)
                                                                                                                                                                                                                          gTCAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaaca          >  1:117574‑M1/1‑139 (MQ=255)
                                                                                                                                                                                                                            cAGAGACATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagc        >  2:11069‑M1/1‑137 (MQ=255)
                                                                                                                                                                                                                                gaCATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGcc                                                                                                                                                          >  2:11347/1‑100 (MQ=255)
                                                                                                                                                                                                                                gaCATAACTGGCACGTAGGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGcc                                                                                                                                                          <  1:11347/100‑1 (MQ=255)
                                                                                                                                                                                                                                 aCATAAATGGCCCGTAGGGT‑TTGCAACCACTAACACAACAATAGAGGGGAAGCTAGGGCGTTAATCTCCCATACTTAACCTGGTCTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGCTGATACGTTCtcaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagattaacgtcgcggaaggtatggctgttggggacaagcattacc  >  2:517360‑M1/1‑132 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           
GCTTTCACGACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAG‑CAGCAGTCAGAGACATAACTGGCACGTAAGGT‑TTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTG  >  NC_000913/2985486‑2985960

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