Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,862,672 G→A 25.9% pseudogene (931/1617 nt) glvC ← pseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6720GA25.9% 42.8 / 12.5 59pseudogene (931/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base G (19/22);  new base A (11/2);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 2.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-02

GGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAG  >  NC_000913/3862439‑3862904
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
ggTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGAAACTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTGCGGAATCAgacag                                                                                                                                                                                                                                                     >  1:764632/1‑234 (MQ=255)
      catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGGATCAACAAGCCCGcc                                                                                                                                                                                                                                             >  1:402749/1‑238 (MQ=255)
           tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTTCGGGAATCAGCCAACCCCCCCCcttt                                                                                                                                                                                                                                     >  1:209819/1‑239 (MQ=255)
                gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGAAACAACAAACCCCCCACTTTTTcccc                                                                                                                                                                                                                               >  1:260670/1‑238 (MQ=255)
                 cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCCAACAGCAACGGCGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGAAATTCCCACCAGCATGGCGGGTAAGGGTGCCGGGATCCACAAGCCCGCCACATTTAcccc                                                                                                                                                                                                                               >  2:1061460/1‑237 (MQ=255)
                 cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTGCGGGGAACAAGAAGACCGCCCCtttttaccca                                                                                                                                                                                                                               >  1:1064912/1‑233 (MQ=255)
                 cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTGAGAGTTCTGGAAACAGCaaaactcccccactcacccca                                                                                                                                                                                                                                  >  2:596831/1‑220 (MQ=255)
                     gatGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAACTTCCAGCGGTGCGGACATTCACACCAGCATGGCGGGAAAAAGTTGCGGGAAACAGCCAGCCCGCACCttttcaccgtgtt                                                                                                                                                                                                                            >  2:768914/3‑230 (MQ=255)
                        gCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGCTAAGGGTTGCGGGAATACCAAAACCCCCCCCTTTATCCGCTTTTCCg                                                                                                                                                                                                                       >  1:430932/1‑239 (MQ=255)
                          tGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGAACTGCAAAACCCCCCCTTTTTACCGGATTTCCCgt                                                                                                                                                                                                                      >  1:1206568/1‑238 (MQ=255)
                                tCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGGAAAGGTGCCGGGCAACACCAACCCCCCCCCTTTACCCCGATTTTCCgtgccgcc                                                                                                                                                                                                                >  2:442019/1‑232 (MQ=255)
                                 cGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTCGGGGATAAAACAAACCGCCCCCTTTTACCCATTTTTCTGGGGGGGcg                                                                                                                                                                                                               >  1:142706/1‑238 (MQ=255)
                                  gcaaagaaGTTTTGCGATAAAGCGTAGCGAATCACGCCCACGCCGGTGTTCCCCCCCACAAAGAAGAGAAAAAGCACGCGCGAGATTGAGGCCACGAGAAGGGCGTGTACCGAAAACGAAAGCGGTAAAATGAGGAGGAGGATAAATTCCATGGGTGCGGAATCCCCCACCAGCAGGGTGGTTAAGGTGGGGGAATTCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCgg                                                                                                                                                                                                                 <  2:222549/228‑1 (MQ=255)
                                      tCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGATCAACAAGACCGCCACTTTTTCCCCATTTTTCGGTGCggggaggaag                                                                                                                                                                                                         >  1:712708/1‑231 (MQ=255)
                                             ttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTTAGGGGATAAAGAAAACCGCCAACTTTTCCCCATTTTCCGGGGCGGGCGGGAAGTaaacc                                                                                                                                                                                                   >  1:595969/1‑234 (MQ=255)
                                              tgcgcaaaaaCCTAGTCAATCAACCCCCCCCCCATGTCCCCCACCACACCAAAGAGAAACATTACGGACGACAATGGGGCCGCCACCACGGCGTGCACCGAAAGCAGCAACGGTAAAATGAAGAGGAAGGTAATTCCCAGCGGTTGGGTACCCCCCACCAGCAGGGGTGTTGAGGTGGGGGAAATAACCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGc                                                                                                                                                                                                <  2:402749/233‑1 (MQ=255)
                                                cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGGGGGGATCAACAAAACCCCCACTTTTTACCGGGTTTCCCGTGCGGGAGGGAAGTACACCGccc                                                                                                                                                                                              >  1:926009/1‑237 (MQ=255)
                                                cGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGGAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATTAGCAAGCCCCCCACTTTTATCCCGTTTTCCGTTGCGGCatgggaagacaatccca                                                                                                                                                                                               >  1:347422/1‑223 (MQ=255)
                                                             tCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGGATCAGCAAGCCCCCCACCTTTTCCCCGTTTTTCCGTGCGGCCGTGAAGGAAAACGCTCAAGAAAAGCCCa                                                                                                                                                                               >  1:785806/1‑238 (MQ=255)
                                                                      ccgcggcgcgtGTTGCCCCCCACACAGAGGAGATAGATTGCGTTGCATGTTGGGGCCCCCATCGCGGGCTATGCCGCAAACAGAAGGGGGGAAGTGATCAGGAAGTGAAGTCCCAGCGGCGCGGTATGTCCCCCGAGCAGGGGTTTTAAGGTGGAGGAAATCAGCATCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCCCGGCGCCaa                                                                                                                                                                    <  2:112315/228‑1 (MQ=255)
                                                                          cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGAATCAACAAGCCCGCCACATTTTACCCGTTTTTCCGGGCGGCCGGGAAAGACATAGCTAAAAAAAAACCCACCGGCACAAAAGAt                                                                                                                                                               >  1:547921/1‑239 (MQ=255)
                                                                          cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGATGGGGGGATCAGCAAACCCCCCCCCTTTACCCCGTTTTTCGGTGCGGGAGTGGAGGAAAACGCGTAAGGAAATCCCACgggccccaagcat                                                                                                                                                                >  2:291853/1‑228 (MQ=255)
                                                                          cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGATGCGGGGATCAGCAAACCCCCCCCTTTTTCCCGATTTTTCCTGGGGGCAGGGAAGGACATCGCCTAAGAAAAACCCCCCGCCCCCAAGaat                                                                                                                                                               >  1:1054497/1‑237 (MQ=255)
                                                                                       caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTCGGGAATCAGCAAACCCACCACTTTTACCCGATTTTTCGGTTGGGCAGTGGAGTACAACACTAAAAAAAAGCCCACGGCTCCCAAGACTTTTTAGGTaaccg                                                                                                                                                 >  1:1003611/1‑235 (MQ=255)
                                                                                                  agagATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGATGCGGGGAACCGCAAACACCCCACCTTTTACCCACTTTTCGGTGCGGGACGGAAAGACAACCACAAAGAAAACCCCACAGCACCCAAGACTTTTTGGGTACCTATCCAGGGCaaa                                                                                                                                      >  2:1127342/1‑239 (MQ=255)
                                                                                                      aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTTCGGGAGACAGCAAACCCCCCCCCTTTACCCCCTTTTCCGCTGGGGCAGGGCAGTACAACCCCCAAGAACAACCCCCCGCCACCAAGACTTTTTTGGTACCATGCCGGGGCAAAAcc                                                                                                                                   >  2:178141/1‑238 (MQ=255)
                                                                                                       tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTTCGGGGAGCAGCAAGCCCGCCCCCTTTTACCCAGTTTTTGGTGCGGCAGGGAAAGAAAACGCCCAAACAATGCCCCCGGCCACAAAAATTTTTTAGTTAccagccaagagcaaaccgc                                                                                                                                  >  2:265192/1‑221 (MQ=255)
                                                                                                            tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGGATCAACAAACCCCCCCCTTTTACCCGATTTTTCCGTGCGGGAGGGAAGAACATCGCTAAAAAAATGCCCCCGGCACCAAAGATTTTTGTGTTACCGCGCAGGGCaaaaaacccctccc                                                                                                                           >  1:1025253/1‑229 (MQ=255)
                                                                                                            tACGGTCGACATTGAGGCCGCCAGCACGGCGTGGACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAAATCCCACCAGCATGGCGGATAAGGTTGGGGGGATGCGCCAGCCCGCCCACTTTAACCCATTTTACGGGGGGGGAGGGGAGGACACTACAACAAAAAAGCCCCACGCCACCAAGAGTTTTGAAAGACCGTGCGCGGAAAAAcacccctcac                                                                                                                           >  2:404026/1‑230 (MQ=255)
                                                                                                              cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAAGTCCAGCGGTTCGGTAATTCCCACCAGAATGGCGGATAAAGGTTCCGGAAATACCAACCGCCCCACCTTTACCCCGCTTTTCTGGTAGGCGGTTTAATAGAATCGCAAAAGAAATACACCGGCACCCACGGATTTTTGAGTTATCTTCCGGGCaaaaacacccctccgg                                                                                                                         >  2:503344/1‑227 (MQ=255)
                                                                                                                 tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCAACAGCATGGCGGGTAAAGGTTGGGGGAATAACAAAACCCCCCCTTTTACCCCATTTTTCGGGGCGGGGGGGGAGCAAATCGCTAAAAAAATATCCACCGCCCACCAGAGTTTTGTGGTACCGCGCAGGGCaaaaacacactcccgggac                                                                                                                      >  2:560639/1‑224 (MQ=255)
                                                                                                                    aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCCCCCCTTTTACCCGATTTTTCGGTGCGGCAGGGGAGTACATCGCTAAAGAAAAGCCCCCGGCACCCAAGATTTTTTGGTTAACGTGCAGGGGAAAACCGCCCTTCCGGaaaaa                                                                                                                   >  1:767567/1‑236 (MQ=255)
                                                                                                                             ccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTTGGGGAATCAGCAAACCCGCCCCCTTTAACCGGTTTTTCGGGGCGGCAGGGAAGTACATCGCTCAAAAAAAGCCCACGGTCCCCAAGGTTTTTTGGGTAACGGGCAAGGCAAAAACCCCTTCCCGGGaaaaagcttttcgc                                                                                                         >  2:1194356/1‑230 (MQ=255)
                                                                                                                                                 aaaaaCCGCAGCGGTAAAATGAAAGGAAAGTTAAATCCCAGCGGTTCGTTAATCCCCCCCAGCAGGGCGGTTAGGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                                    <  1:47914/237‑1 (MQ=255)
                                                                                                                                                  cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCACCTTTACACGCTTTTCAGGGGCGGCAGTGAAGTACATCGCTAAAAAAAAGCCCACGGCACCcaaaa                                                                                                                                                                 >  1:1041919/1‑163 (MQ=255)
                                                                                                                                                  cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGAATCAGCAAACCCCCCACCTTTAACCGATTTTTCGGTGCGGCAGTGAAGGACACCGCCAAAAAAAATCCCCCGCCCCCCAAGATTTTTTTGTTACCGTGCAGGGCAAAAACCCCCTCCGGGAACAAAGATTTTAGCGGCTCGGCGCGCAAAccgcac                                                                                   >  2:834236/1‑234 (MQ=255)
                                                                                                                                                   aaaCAGCAACGGTGAAGTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCACCGGCATGGCGGTTAAGGGTGCGGAAATCAGCAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAAct                                                                                  <  1:1127342/239‑1 (MQ=255)
                                                                                                                                                     aCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCGCCACTTTTACCCGGTTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCa                                                                                                                                                                     >  1:183287/1‑158 (MQ=255)
                                                                                                                                                     aCAGCAACGGTAAAAGAAAAAGGAAGGAAAATCCCAGGGGTCCGGTAATCCCCCCCAGCTGGGCGGTTAGGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAActct                                                                                <  1:291853/239‑1 (MQ=255)
                                                                                                                                                                gAAAGGAACGGGAAGAAAAATTCCAGGGGTCCGGTTATTCCCACCAGCATGGCGGTAAAGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                      <  1:1106069/238‑1 (MQ=255)
                                                                                                                                                                 aaaTAAACAGGAAGTAAAATTCCAGCGGTTCGGAAATCCCCACCAGCATGGCGGTAAGGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                      <  1:8900/237‑1 (MQ=255)
                                                                                                                                                                   aataaCAGGAAAGTAAATTCCCGGCGGTCGGGTAATCCCAACCAGCAGGGCGGTAAGGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                   <  1:660197/236‑1 (MQ=255)
                                                                                                                                                                    tGAACGGGAGGTAAAATCCCAGGGTGTGGGTAATTCCCCCCAGCTGGCGGTTTAGGGTTGCGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                 <  1:768931/239‑1 (MQ=255)
                                                                                                                                                                    tGAACAGGAATGAAAATCCCAGCGTTTCGATAATCCCCCCCAGCATGGCGTGTAAGGTTGGGGGAATCAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                 <  1:1194356/239‑1 (MQ=255)
                                                                                                                                                                      aaCAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGACCCCCACTTTTACCCGATTTTTCGGTGCGGCAGGGGAGTACATCGCCAAAAAAAATCCCCCCGCACCCAAAATTTTTGAGGTACCGTGCAGGGCAAAAAACCCTTCCGGGAACAACCATTTTAGCGGCCCGGGGCTCAAACGGAACTCtttccgcagcggggccc                                                                >  2:665205/1‑223 (MQ=255)
                                                                                                                                                                         aGGAAGTAAATTTCCAGCGTTGCGGTAATCCCCACCAGCTGGGCGTTTAAGGTTGTGGAACTCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                            <  1:106471/239‑1 (MQ=255)
                                                                                                                                                                                    atcccAGCGGTGCGTTACTCCCCCCCAGAAGGGCGGTTAGGGTTGCGGAAATCAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                                  <  1:740045/236‑1 (MQ=255)
                                                                                                                                                                                       ccagtggtttgtaacatCCCACCAGCAGGGCGGTCAGTGTGGGGGGAACAGGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccg                                               <  1:768914/223‑1 (MQ=255)
                                                                                                                                                                                                  ggacaccaccccacggcggggggTTAAGGGTGGGGGTAACCGGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACa                                      <  2:260670/218‑1 (MQ=255)
                                                                                                                                                                                                          cccccgCATGGCGGTTAAGGTGGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                            <  1:665205/234‑1 (MQ=255)
                                                                                                                                                                                                            cgagatgggCGGTTAAGGTTGCGGGAAGCAGCAGCCACGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                          <  1:657388/232‑1 (MQ=255)
                                                                                                                                                                                                            ccaGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                          <  1:721084/239‑1 (MQ=255)
                                                                                                                                                                                                              ctctggcagggggttgagtgtggCGGAAATCAACAACCCCGCCACTTCTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                            <  2:933716/217‑1 (MQ=255)
                                                                                                                                                                                                              ccaggttgctgTTTAGGCTGGGGGAAATCAGAACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAATTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                          <  1:1182449/229‑1 (MQ=255)
                                                                                                                                                                                                                   ggCGGTTAAGGTGGCGGGAATAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATa                    <  1:442019/238‑1 (MQ=255)
                                                                                                                                                                                                                     cGGTTAAGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                  <  1:842811/238‑1 (MQ=255)
                                                                                                                                                                                                                      gggTTAAGGTGGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                  <  2:1025253/237‑1 (MQ=255)
                                                                                                                                                                                                                       ggTTAAGGTTGCGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAAc                 <  1:834236/238‑1 (MQ=255)
                                                                                                                                                                                                                                 ggtGCGGGAACAACCAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt          <  2:209819/236‑1 (MQ=255)
                                                                                                                                                                                                                                          aaTCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg  <  1:86328/238‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
GGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAG  >  NC_000913/3862439‑3862904

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