Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,017,393 C→T L21L (CTG→TTG)  rmuC → putative recombination limiting protein RmuC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,017,3930CT98.5% 214.3 / NA 67L21L (CTG→TTG) rmuCputative recombination limiting protein RmuC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (36/29);  minor base G (0/1);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGCTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCACGCGGATGGAAGCCGCACAGCAACATGCTGACGATAAAATTC  >  NC_000913/4017146‑4017636
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
gCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGt                                                                                                                                                                                                                                                  >  2:197851/1‑251 (MQ=255)
   aTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTg                                                                                                                                                                                                                                               <  1:150277/251‑1 (MQ=255)
      cGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGcc                                                                                                                                                                                                                                             >  2:160263/1‑250 (MQ=255)
      cGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGACTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGcc                                                                                                                                                                                                                                             >  1:485748/1‑250 (MQ=255)
         tGAAAATCGTGGTGGCAGCGCCGCGTCGTCTGCTGTAACTCTCTTCTCCTGTCTGCCGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGCCGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTACGTCAGGCGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGtt                                                                                                                                                                                                                                         >  2:64049/1‑251 (MQ=255)
         tGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGtt                                                                                                                                                                                                                                         >  1:576975/1‑251 (MQ=255)
                cGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTTCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAAc                                                                                                                                                                                                                                   <  1:357559/250‑1 (MQ=255)
                         aGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCgcagca                                                                                                                                                                                                                         <  1:118781/251‑1 (MQ=255)
                                cgtcgtCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGcc                                                                                                                                                                                                                  >  2:474214/1‑251 (MQ=255)
                                 gtcgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCg                                                                                                                                                                                                                 <  1:234717/251‑1 (MQ=255)
                                   cgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAg                                                                                                                                                                                                               <  2:324298/251‑1 (MQ=255)
                                       tgctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAAt                                                                                                                                                                                                           >  1:615852/1‑251 (MQ=255)
                                        gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGTCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAAtt                                                                                                                                                                                                          >  1:53469/1‑251 (MQ=255)
                                         cggTAATTCTCATCACCTGACTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGACTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAACACAAGTCTATTTTGTGCGGGTAAGTTGTGGCGTCAGGAGGCGTTGTGGAGTTCGCAATCATTGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAAATGGCGGGTTGTTTGCCAGTTATCTACAAGCGCGGAAAAGAGCCGAGCAATTa                                                                                                                                                                                                         <  2:318575/249‑1 (MQ=255)
                                                      ctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGa                                                                                                                                                                                                                >  2:478924/1‑231 (MQ=255)
                                                      ctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGa                                                                                                                                                                                                                <  1:478924/231‑1 (MQ=255)
                                                         cTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAag                                                                                                                                                                                            >  2:624067/1‑248 (MQ=255)
                                                               gAAGGCCGACGCGTTCGGCCTTTTGCATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGCGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATCGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCACTTAGCTGAACGTGAAGAGATGGTcg                                                                                                                                                                                   >  2:104723/1‑251 (MQ=255)
                                                                    ccGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAg                                                                                                                                                                              <  1:624067/251‑1 (MQ=255)
                                                                             tCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACTCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGGTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTATgcgc                                                                                                                                                                      <  2:100170/250‑1 (MQ=255)
                                                                              cGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGGGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGACGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCgg                                                                                                                                                                    >  2:285550/1‑251 (MQ=255)
                                                                              cGGCCTTTGGTATTTATGCGTAGCTCCGCGCAGGAAATGCATTTCCAACTGGGCGTCTGTACAGCACAATTCTATTTCGTGCGTGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTTAATCATGGTTTACGCAGTTAGTGCGTTGCTGTGTGTGGCAATTGGCTGGTGGTTTGCCAGTTAGCAAAATGCGTAGCAAAAAGCCTAGGAATTAGCTGAACGTGAAAAGTTGGTCGTGGAGTTAAgcgccg                                                                                                                                                                    <  2:401639/251‑3 (MQ=255)
                                                                                gCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCACCTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCAAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGc                                                                                                                                                                   >  1:263276/1‑250 (MQ=255)
                                                                                gCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAACGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGc                                                                                                                                                                   >  2:253814/1‑250 (MQ=255)
                                                                                    tttGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCaaaa                                                                                                                                                               <  1:637871/250‑1 (MQ=255)
                                                                                    tttGTATTTTTGCGAAGCGCATCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAaac                                                                                                                                                              >  1:198565/1‑251 (MQ=255)
                                                                                      tGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAAt                                                                                                                                                                                                           <  1:487295/204‑1 (MQ=255)
                                                                                      tGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGACGAGCAAt                                                                                                                                                                                                           >  2:487295/1‑204 (MQ=255)
                                                                                       gTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCCGAACGTGAAGAGAGGGTCGCGGAGTTAAGCGCGTCAAaacaa                                                                                                                                                            <  2:67994/250‑1 (MQ=255)
                                                                                       gTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGGGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTGATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGTAGTTAAGCGCGGCAAaacaac                                                                                                                                                           <  2:123433/251‑1 (MQ=255)
                                                                                               gCGTAGCCCCTCGCAGGAAATGCCTTTCCAACTGGACGGTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGTTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCATCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTAc                                                                                                                                                    <  2:530163/250‑1 (MQ=255)
                                                                                                cGTAGCGCCTCGCAGGACCAGCCTTTCCAACTGGACGTTTATACAGCACAATTCTATTTTGTGCGGGCACGTTGTCGCGTCAGGAGGCGTACTGGATTTATCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTTGCTGCATGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGTACGTGCAACCATGGTCTCCGAGTTAACAGAGGCAAAACTACAAATTaacc                                                                                                                                                  >  2:573218/1‑248 (MQ=255)
                                                                                                      gcCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCg                                                                                                                                            >  2:282977/1‑251 (MQ=255)
                                                                                                          cGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAg                                                                                                                                          <  2:660236/249‑1 (MQ=255)
                                                                                                                 aaTGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGAGCGGGCAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGATGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGTGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAACGAGCACTGGCGt                                                                                                                                 >  1:289296/1‑251 (MQ=255)
                                                                                                                  aTGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTg                                                                                                                                >  2:202340/1‑251 (MQ=255)
                                                                                                                  aTGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGAGGATAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTTTCAATCATGGTTTAAGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTGAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTg                                                                                                                                >  2:290497/1‑251 (MQ=255)
                                                                                                                         tCCAACTGGACGTTTGCACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCGGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAAATACCCAAAGCGAGCACTGGCGTGCAGAGTg                                                                                                                         >  2:71314/1‑251 (MQ=255)
                                                                                                                               tGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGtt                                                                                                                   >  2:270077/1‑251 (MQ=255)
                                                                                                                                 gACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACTCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTAc                                                                                                                 <  1:422329/251‑1 (MQ=255)
                                                                                                                                      ttGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTAc                                                                                                                 <  1:581856/246‑1 (MQ=255)
                                                                                                                                      ttGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTAc                                                                                                                 >  2:581856/1‑246 (MQ=255)
                                                                                                                                            aGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCaaa                                                                                                                                                                >  1:109590/1‑193 (MQ=255)
                                                                                                                                            aGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCaaa                                                                                                                                                                <  2:109590/193‑1 (MQ=255)
                                                                                                                                                 aaTTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATGAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTgcg                                                                                                 <  1:168143/251‑1 (MQ=255)
                                                                                                                                                  aTTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTgcgc                                                                                                >  1:152797/1‑251 (MQ=255)
                                                                                                                                                              gCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGt                                                                                    <  2:161377/251‑1 (MQ=255)
                                                                                                                                                                           tgCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACctctc                                                                       >  2:166134/1‑251 (MQ=255)
                                                                                                                                                                           tgCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACctctc                                                                       >  2:355060/1‑251 (MQ=255)
                                                                                                                                                                             cGTCAGGAGGCGGTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACTAAGTGCGCAGCCTGCAAAGTATTAACACCGCTCTg                                                                     <  2:29673/251‑1 (MQ=255)
                                                                                                                                                                                  ggaggCGTTGTGGATTTCTCAATCATGGATTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCACTTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATACCGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGc                                                                >  2:380892/1‑251 (MQ=255)
                                                                                                                                                                                   gaggCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGcc                                                               >  2:613024/1‑251 (MQ=255)
                                                                                                                                                                                     ggCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGa                                                             <  1:656159/251‑1 (MQ=255)
                                                                                                                                                                                      gCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGAt                                                            >  2:212622/1‑251 (MQ=255)
                                                                                                                                                                                          tgtgGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTg                                                         <  1:576810/250‑1 (MQ=255)
                                                                                                                                                                                                        tCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCCATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGCTAAGCGCGGACCAACAACAAATTACCCAAAGCGAGCACTGGCGCGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACATCTCTGGAGGCCGATTGGCGTGAAGTAACAAcg                                          >  2:421865/1‑251 (MQ=255)
                                                                                                                                                                                                         cATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCAcg                                          <  1:244619/250‑1 (MQ=255)
                                                                                                                                                                                                             gTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGCGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCACGCGGAt                                     >  1:401641/1‑251 (MQ=255)
                                                                                                                                                                                                                    cagatATTGAGTTGGTGGGAGTGGCAATTGGCTGGTTGTTTGGCAGTTATCAACATGCGCGGCAAAAAGCCGAGCAATTAGCTGAACGGGAAGAGATGGACGCGGAGTTAAGCGCGGCAAAACAAAAAATTACACAAAGCGAGAACTGGCGTGCAAAGTGCGAGTTACTTAATAACAAATGGAGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGGAACAAGGCGGATGGANgac                              >  2:600029/5‑249 (MQ=255)
                                                                                                                                                                                                                      gTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCACGCGGATGGAAGCCGc                            <  2:401641/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gtgGCAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCACGCGGATGGAAGCCGCGCAGCAACATGCTGACg           <  1:565628/250‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCAATTGGCTGGTTGTTTTCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAATAGATGGTCGCGGAGGTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCACGCGGATGGAAGCCGCACAGCAACATGCTGACGATa        <  2:193262/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cAATTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCACGCGGATGGAAGCCGCACAGCAACATGCTGACGATaaa      <  2:53469/251‑1 (MQ=255)
                                                                                                                                                                                                                                             aaTTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCACGCGGATGGAAGCCGCACAGCAACATGCTGACGATaaa      >  1:164427/1‑250 (MQ=255)
                                                                                                                                                                                                                                              aTTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCTCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGc                                                                                         <  2:497227/166‑1 (MQ=255)
                                                                                                                                                                                                                                              aTTGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGc                                                                                         >  1:497227/1‑166 (MQ=255)
                                                                                                                                                                                                                                               ttGGCTGGTTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCACGCGGATGGAAGCCGCACAGCAACATGCTGACGATAAAAtt   >  1:459633/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 ggctggttGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCACGCGGATGGAAGCCGCACAGCAACATGCTGACGATAAAATTc  <  1:282977/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
GCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTGGTGGGTGTGGCAATTGGCTGGCTGTTTGCCAGTTATCAACATGCGCAGCAAAAAGCCGAGCAATTAGCTGAACGTGAAGAGATGGTCGCGGAGTTAAGCGCGGCAAAACAACAAATTACCCAAAGCGAGCACTGGCGTGCAGAGTGCGAGTTACTCAATAACGAAGTGCGCAGCCTGCAAAGTATTAACACCTCTCTGGAGGCCGATCTGCGTGAAGTAACCACGCGGATGGAAGCCGCACAGCAACATGCTGACGATAAAATTC  >  NC_000913/4017146‑4017636

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