Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,421,767 G→T 20.6% E149* (GAG→TAG)  yhdY → putative amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,421,7670GT20.6% 109.6 / 2.5 61E149* (GAG→TAG) yhdYputative amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (5/42);  new base T (12/2);  total (17/44)
Fisher's exact test for biased strand distribution p-value = 2.65e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCGTTGTTTATGGCAGAA  >  NC_000913/3421532‑3421994
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        
cacaAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCCCCAATGTTCTGGAAAATACTCACGCATCGCGGTCGCTATATTGCGGCCTGGGCGG‑GATTTACCCACTGATTGTCTGGTGGCGTGTGTTTGGC‑GGTTTTTTTGTGC‑TTgtgag                                                                                                                                                                                                                                      >  1:557520/1‑235 (MQ=255)
           tttgcgtgtgtggtttCATCCGCTGACGATTTGGTCAGTTTATGTTTGGGGTTTACCCGACCGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAAcc                                                                                                                                                                                                                           <  2:628313/225‑1 (MQ=255)
            ttccTTTTGGGTTTTTATCCACGAACGGTTTGGTAAGTTTATGGTTGGGCTTTACCCACCCGCCCAACGTTGGCGAATTAACCTCGCATTACTGAGTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAAcc                                                                                                                                                                                                                           <  2:717587/236‑1 (MQ=255)
                tGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGCGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCa                                                                                                                                                                                                                      <  1:48268/239‑1 (MQ=255)
                       tCTTCATCCACGCACGATTTGGGCAGTTTATGTATGGGCTTTACCCACAAGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGGTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATgggg                                                                                                                                                                                                                <  1:666997/238‑1 (MQ=255)
                       tCTTCACCCACGAATGATTTGGTCAGTTTATGTATGGGCTTTACCCAACCGACCAAGGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTGTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATgggg                                                                                                                                                                                                                <  1:376565/238‑1 (MQ=255)
                       gtctCAGCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGGGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATgggg                                                                                                                                                                                                                <  1:883286/235‑1 (MQ=255)
                        cTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAAGTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGACAATACTCCCGCATCGCGGTCGCTATATTTCGGCCGGGGCGGTGATTTACCCACTGATTCTCTGGTGGCTTATGTTTGGG‑GGTTTTTTTGGTG‑TTGAAGCGGGtgaaagacggggaagggggg                                                                                                                                                                                                               >  2:614786/1‑220 (MQ=255)
                            tagttgtttataagggactggTTTATGTTTGGGCTTTACCCACCCGGACAACGCTGGCGAATTAACCTCGCATTTCTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCt                                                                                                                                                                                                         <  2:220459/219‑1 (MQ=255)
                            ctCCCCGAACGGTTTGGTTATTTTTTGTATGGGCTTTCCCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGAGGGTTGAAACCCGGCAATGGGGC‑GGGCt                                                                                                                                                                                                         <  2:322683/238‑1 (MQ=255)
                                 cgaacgaTTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGTC‑GGTTTTTTTTCTC‑TTT‑AGCGGGTTGAAACCCGGCACCGgggg‑gggcgtcgc                                                                                                                                                                                                     >  1:771006/1‑233 (MQ=255)
                                  tgatggttttgttagttttattatgGGTCTTACCCCCCGGACCAGCCCTGGCCAATTAACCTCGCATTAATGATTGGGCTTGTGTCGATCGCCACAATGTTGTGGAAAATACGCCCGCATCGCGTGCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCTGCAATGGGGC‑GGGCTGACGCt                                                                                                                                                                                                   <  2:780782/218‑1 (MQ=255)
                                  gaacgtTTTGGCTAGTTTATGTTTGGGCTTTACCCACCCGACCACCGCCTGCCAATTAACCTCGGATTACCGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATTCTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCt                                                                                                                                                                                                   <  2:735700/239‑1 (MQ=255)
                                   acctaTTTTGTCCTTTTTTGTTTGGGCTTTGCCCCCCCTACCAACGCGGGCGAATTAACCTCGCATTACGGATTGGGCTTGTGGCGTTCGCACGAGTGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTa                                                                                                                                                                                                  <  2:161784/235‑1 (MQ=255)
                                    acgaTTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGGGATTTACCCACTGATTGTCTGGTGGCTGAGGTATTGC‑GGTTTTTTTTTTC‑TTT‑AGAGGGGTGAAACACCGGGATggggg‑ggggggggggcaa                                                                                                                                                                                                 >  2:591167/1‑229 (MQ=255)
                                        tttGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTTCGGCCTGGGCGGTGATTTACCCACTGCTTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTTCTCTTTT‑TGCGGGTTGAAACCCGGAAATGgggc‑gggcgaggctcacctt                                                                                                                                                                                              >  2:679181/1‑227 (MQ=255)
                                              cAGTTTTTTTTTGGGCTTTACCCCCACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCCCCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAattat                                                                                                                                                                                      <  2:221262/239‑1 (MQ=255)
                                                    ttgttataagagacagcaccCGGCCAACGGTGGCGAATTAACCTCGCTTTACTGATTGGGCTTGTGTCGGTCGCCCCAATGTTCTGGAAAATTCTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCa                                                                                                                                                                                  <  2:597841/221‑1 (MQ=255)
                                                       tatgGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGTTC‑TTT‑AACGGGTTGTAACCCGGCAATTGGGC‑GGGCTGGCGGTAA‑CTTTAATTATTGCAACAG‑t                                                                                                                                                                            >  1:200515/1‑239 (MQ=255)
                                                           ggCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTT‑AACGGGTGGGAACCCGGCAATGGGGC‑GGGGCGACGCTAA‑CTTTTATTATTTCAttcggtggg                                                                                                                                                                         >  1:654356/1‑231 (MQ=255)
                                                                      cacGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGCTGTATGGCGGGTTTTTTTTTTT‑TTG‑GGCGGG‑TGAAAACCCGCACTGGGGG‑GGGCTGCGCCTAACCTTTTTTTTTTGCA‑CAC‑TTGGGAGTGggggggg                                                                                                                                                             >  2:73413/1‑238 (MQ=255)
                                                                       acGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTTCTC‑TTT‑AGCGGGTTTGAACCCCGCAAAGGGGC‑GGGGTGACGCTAA‑CTTTTATTATTGCATCAG‑TTGGGATTTTgggggg                                                                                                                                                             >  1:872886/1‑237 (MQ=255)
                                                                       acGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTTCTC‑TTG‑AGCGGGTTGAAACCCGGCCATGTGGC‑GGGCTGACGCCAA‑CTTTAATTATTGCTTCCG‑TTTGGGTTGAgggggg                                                                                                                                                             >  1:675617/1‑237 (MQ=255)
                                                                       aagACCAACGTTGGAGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCGCCAATGTTCTGGAAAATTCTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTGCCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGgcg                                                                                                                                                            <  1:614786/237‑1 (MQ=255)
                                                                          accaccGGTTGCGATTTAACCTCTCATTACTGATTGGGCTTGGTGCGATCGCACCAATGTTGTGGAAAATCCTCCCGCATCGCGGTCGGTATATTGCGGCCTGGGCGGTGATTGACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTa                                                                                                                                                         <  2:394793/239‑1 (MQ=255)
                                                                            caccGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTAc                                                                                                                                                        <  1:62739/235‑1 (MQ=255)
                                                                               cgcgGGCGAATTAACCTCGCATTACTGTTTGGGCTTGGGTCGTGCGCCGCAATTTTCTGGAAAATACTCCCGCATCGCTGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCGCTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTg                                                                                                                                                    <  2:831346/239‑1 (MQ=255)
                                                                               cGCTGGAGAATTAACCTTGCATTACTGAATGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCTTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTg                                                                                                                                                    <  1:648790/239‑1 (MQ=255)
                                                                                  tGGCGAATTATCCTTGCATTACTGATTGGGTTTGGGTCCATCGCACCAATGTTCTGGAAAATTCTCCCGCAGCGCGTTCGCTTTATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTgggg                                                                                                                                                 <  2:204044/239‑1 (MQ=255)
                                                                                    gccaaTTACCCTCGCCTTACCGATTGGGCTTGGGGCGGTCGCACCAATGTTCTGGAAAATTCTCCCGCAGTGCGGTCGGTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGa                                                                                                                                               <  2:672324/236‑1 (MQ=255)
                                                                                         ttAACCCTGCATTTCCTAGTGGGCTTGTGTCGATCGCCCCAATGTTCTGGAAAATTCTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTa                                                                                                                                          <  1:773668/239‑1 (MQ=255)
                                                                                                  cATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTT‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGAAGCTAA‑CTTTAATTTTTGCATCAG‑TTTGGGTTGCGGGGGCGCTACCGTGGGGGGTCTTACTGGCGTTa                                                                                                                                 >  1:740351/1‑239 (MQ=255)
                                                                                                           ttGGGCTTGTGTCCATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTcgccg                                                                                                                         <  1:35881/238‑1 (MQ=255)
                                                                                                            tGGGTTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGGGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCt                                                                                                                       <  1:893135/239‑1 (MQ=255)
                                                                                                                 ttGTTTCGATCGTCCCAATGTTCTGGACAATCCTCCCCCATTGCGGTCGGTTTCTTGCGGCCTTGGCGTTGATTTCCCCTCTGCTTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGTTGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCa                                                                                                                   <  2:294499/238‑1 (MQ=255)
                                                                                                                 ttGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTTCTC‑TTG‑AGCGGGGTGAAACCCGGCAATGGGGC‑GGGGGGACGCGAA‑CTTTTTTTCTTGCATCCA‑GTGGGGTTTCGGGGGGGCTACCCTGTGGGGACTTACTGGCTCTAgggggcgcgccccac                                                                                                                  >  2:124415/1‑226 (MQ=255)
                                                                                                                  tgtgTCGCTCGGCCCAATGTTCGGGAAAAGTGCCCCGCCTCGTGGTCGCTATATGGCGGCCGGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTAAGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCata                                                                                                                 <  2:654356/239‑1 (MQ=255)
                                                                                                                    tgtCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGCGTTTTTTCTC‑CTT‑ACGCGGTTGAAACCCGGCAATGGGGCGGGGGGGACCCT‑A‑CTTTCATTATTTCGTCCGATGGGGCTTGGGGGGGCGC‑CTCGTGGGGGGTATTTCTGGCGTTTGGtgcccacgccccaatg                                                                                                               >  2:97354/1‑224 (MQ=255)
                                                                                                                       cgctCGCCCCATTGTTTGGGAAAATTCTCCCGCTTCGTGGTCGGTGTATTGCGGCCTGGGCGG‑GATTTTCCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTT‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGTCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCTCTCCCATA‑TGCCGa                                                                                                          <  2:557520/236‑1 (MQ=255)
                                                                                                                        gATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGa                                                                                                          <  1:591167/239‑1 (MQ=255)
                                                                                                                         aTCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGGTTTTTTGCTC‑TTT‑AGCGGGTTGAAACCCGGCAAATGGGC‑GGGGTGTGGCTCA‑CTTTTATTTTTTTATCCG‑TTGGGGTTTGGGGGGGGCGTCCCTGGGGGGTGTTACTTGCGTTTGGGCGGCGCTcccccaacgcggg                                                                                                          >  2:415144/1‑229 (MQ=255)
                                                                                                                          tCGCCCCCATGTTCGGGAAAATTCCCCCGCATCGGGGTGTCTATATGGGGGCCTGGGGGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TGGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGAtt                                                                                                        <  2:544201/239‑1 (MQ=255)
                                                                                                                                                   cccGCCTCGCGGTCTCTATTTTGGGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑GGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTa                                                                               <  2:740351/239‑1 (MQ=255)
                                                                                                                                                                     tgttGCGGCTTGGGCGGTGTTTTACCCAGTGATTGTTTGGTGGCTGATGTATGGC‑GTGTTTTTTGCGC‑TGG‑AGCGGGTTGAAACCCGGGAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCgcggcg                                                             <  2:675617/237‑1 (MQ=255)
                                                                                                                                                                       ttGCGGCCTGGGCGGTGATTTACCCACGGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTg                                                           <  1:415144/239‑1 (MQ=255)
                                                                                                                                                                       ttGCGGCCTGGGCGGTAATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGGTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTg                                                           <  1:853571/239‑1 (MQ=255)
                                                                                                                                                                                 gggggttttttACCCTTGGATTGTTTGGGGAGTGCTGTTTGGG‑GTTTTTTTTGCTC‑GTG‑GGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGTTGGTAA‑CTTTATTTATTGCATCAG‑GTGGGATTGCGGGGGCGTTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTa                                                 <  2:746695/231‑1 (MQ=255)
                                                                                                                                                                                 ggCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTTCTC‑TTT‑AGCGGGGTGAAACCCGGCAATGGGGC‑GGGGTGTCGCTAA‑CTTTTATTTTTGCGTCAG‑GTGGGGTTTCGGGGGCGCTACCCGGGGGGGGCTTTATGGGGGTAGGGGGCGCGTCCCCTATTGCCGTTTG‑GTGTATGTTTTCGGTAATTTTTTTTGAATTTTGGCGCGGGGGGCGGCGGAtt                                                  >  2:451267/1‑238 (MQ=255)
                                                                                                                                                                                    agttttttttaagagactggTCTGGTTGCTGGTGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGGAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCg                                              <  2:258634/220‑1 (MQ=255)
                                                                                                                                                                                       tagttACCCCGTGGTTGTCTGGTGGGTGATGGATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGtt                                            <  1:73413/235‑1 (MQ=255)
                                                                                                                                                                                       gttttACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTGTGCTC‑TTG‑AGCGGGTTGAAGCCCGGAAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGtt                                            <  1:784249/236‑1 (MQ=255)
                                                                                                                                                                                            aCCCACTGATTGTGTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CCTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGttt                                      <  1:679181/239‑1 (MQ=255)
                                                                                                                                                                                                  tGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCCG‑TTGGGATTGCGGGGGCGCTACCGGGGGGGATCTTACTGGCGTTAGGGCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGGGCCGCCGGTTACCGGTCTGTTTTTTTCt                                >  1:714287/1‑239 (MQ=255)
                                                                                                                                                                                                                  tgttggtGGG‑GGTTTTTTTTCTT‑GTT‑GGGGGGTTGAAACCGGGTCATGGGGG‑GGGGTGACGGTAA‑CTTTAAGTATTGCATGAT‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCg                 <  2:221606/232‑1 (MQ=255)
                                                                                                                                                                                                                  gATGTATGGC‑GTTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCg                 <  1:689867/238‑1 (MQ=255)
                                                                                                                                                                                                                  gATGTATGGC‑GGTTTTTTTGCTC‑TTT‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCCTA‑TGCCGATTGTGCGTATCTTTTCGGTCATTTTTATCGGATTCTGGCGGGGCGGGCCGCTGATTTCCGTTCTGTTTTTgcttgctcttataaacctt                 >  1:258634/1‑220 (MQ=255)
                                                                                                                                                                                                                  gATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAAGGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGACTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCg                 <  1:124415/238‑1 (MQ=255)
                                                                                                                                                                                                                   atttATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCg                 <  1:458296/234‑1 (MQ=255)
                                                                                                                                                                                                                              gTTTTTTTGCTC‑TTT‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTTA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTTGGTTGCCGCTCCCCTT‑TGCCGATTTGGCGGTTTTTATTGGGCATTTTTTTTCAATTTTGGGGCCGCGTGCCCCCGGTTACCCTGCTTTTTTTTTCTTTGGTGCTGCTGCCggtgttgtttggc     >  2:458296/1‑232 (MQ=255)
                                                                                                                                                                                                                               tttttttGCTC‑TTG‑AGCGGGTTGAAACCCGGGAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCGTTGTTTATGGCa     <  1:479991/238‑1 (MQ=255)
                                                                                                                                                                                                                               ttttttgGCTC‑TGG‑AGCGGGTTGACACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTAGTGCATCAG‑TTGGGATTGCGGGGGCGCTCACGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCGTTGTTTATGGCa     <  1:370508/238‑1 (MQ=255)
                                                                                                                                                                                                                                 ttgttGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCGTTGTTTATGGCAGaa  <  1:344363/236‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        
CACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGC‑GGTTTTTTTGCTC‑TTG‑AGCGGGTTGAAACCCGGCAATGGGGC‑GGGCTGACGCTAA‑CTTTAATTATTGCATCAG‑TTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATA‑TGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCGTTGTTTATGGCAGAA  >  NC_000913/3421532‑3421994

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