Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,862,674 C→A 16.2% pseudogene (930/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,6740CA16.2% 79.6 / 2.2 73pseudogene (930/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 C (31/28);  new base A (2/10);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 2.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01

TGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTA  >  NC_000913/3862438‑3862912
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
tGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAACAACGGTGAAATGAACAGGAAAGTAAATTCCAGCGGTTCGGTAATTCCCCCCAGGATGGCGGGTACGGTGGAGGGAGTGAGCaa                                                                                                                                                                                                                                                      >  1:238890/1‑239 (MQ=255)
           catATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGGCCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGGCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGc                                                                                                                                                                                                                                                         <  1:280803/226‑1 (MQ=255)
             tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATACGACAAGCCCCCACtttt                                                                                                                                                                                                                                       >  1:1223493/1‑239 (MQ=255)
             tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGAATCAGCAAGCCCCCCCCtttt                                                                                                                                                                                                                                       >  1:441208/1‑239 (MQ=255)
                  gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGAAATCGCAAACACCCCCCTTTTAccc                                                                                                                                                                                                                                   >  1:953200/1‑238 (MQ=255)
                   cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGAACTAACCAACCCCCCCCCTTTACCCgg                                                                                                                                                                                                                                 >  1:1288243/1‑238 (MQ=255)
                   cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGGAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGCATTCCGAAGCCCCCCCCTTTTTCCCg                                                                                                                                                                                                                                  >  2:1200231/1‑238 (MQ=255)
                     tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTGGCGGGGAACCGGAAACCCGCCACCTTTTccccctt                                                                                                                                                                                                                               >  1:747996/1‑235 (MQ=255)
                       gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATACGACAGCCCGCCACCTTTATCCCAgttt                                                                                                                                                                                                                            >  2:1286358/1‑239 (MQ=255)
                         tGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGGTGGGGGAAGAAGAAACCCCCCCCCTTTTCCCCCAtttc                                                                                                                                                                                                                            >  2:1138748/1‑237 (MQ=255)
                           cTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTTCGGGGAACAGCCAACCCCCCACCTTTTCCCCGTTTTTcgg                                                                                                                                                                                                                         >  1:908361/1‑237 (MQ=255)
                                caccGGGATCAGGTTTTGCGATAAACCCGAGTCAACCAGCCCCCCCCCATTGTCGCCCCCAACCCAAAAGAGAAACATTACGGCCGACATGGAGGCCGCCAGCACGGCGTGTACCGCAAACAGAAACGGGAAAATGAACAGAAAGATAATTCCCGGCGGTTCGTTATTCCCCCCCTGCATGGGGTTTAAGGTTGGGGAAACAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCgg                                                                                                                                                                                                                   <  2:506774/236‑1 (MQ=255)
                                cATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGATGCGGGAATTAACCAACCCCCCCCTTTTACCCCACTTTTCGGTGCg                                                                                                                                                                                                                    >  1:368161/1‑238 (MQ=255)
                                  tCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGATTGCGGGATCAACAAACCCCCCCCCTTTACCCCATTTTTCCGTGCGGCa                                                                                                                                                                                                                 >  2:688805/1‑239 (MQ=255)
                                     ggATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGGTGGGGGAGCACGCAAACCCCCCCTTTTTACCCATTTTTCCgggcgggaggg                                                                                                                                                                                                               >  2:1118299/1‑230 (MQ=255)
                                      gcaccAGTTTGGCGGAAAACCCAATCCATACAGCCCGCCGCCCCTCTCCCCCCCCACACCAAAAAGATCCATTCCGGGAGACTTTGGGGCCGCAAGCCCGGGGGGTCCCGCAAAAAGCAGGGGTAAAATGAACAGGAGGTAAACTCCGCGCGGTCGGTAATTCCCCCCCAGCCGGGCGGTTAGGGTGGGGGCAAAAAGAAAGTCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTg                                                                                                                                                                                                              <  2:892977/235‑1 (MQ=255)
                                          gaGTTGTGCGATAAACCCTACTCAACCCGCCCCCCCCCCATTTCCCCCACCACAAAAAAGAGATACATAACGGTCGACATGGAGGCCGCCAGCAGGGCGTGTACCGCAAACAGCAACGGTGAAAGGAAAAGGAAGGAAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGGGGTTAAGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                                                                                                                                         <  2:908361/238‑1 (MQ=255)
                                          cAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTTCGGGGATACGCAAACCCGCCACTTTTACCCGATTTTTCGGTgggggcgcggaaga                                                                                                                                                                                                         >  1:1141593/1‑229 (MQ=255)
                                           aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGATCAGCAAGCCCGCCACTTTTTCCCCGTTTTTCCGTGCGGCAGGGAAGGAc                                                                                                                                                                                                        >  1:990994/1‑239 (MQ=255)
                                                               tCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATAAGCAAGCCCGCCACTTTTACCCGATTTTTCGGGGCGGCAGTGAAGTACATCGCTAAAGAAATACCCa                                                                                                                                                                                   >  1:353658/1‑238 (MQ=255)
                                                                     cAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGAACCGCAAACCCCCCACCTTTTCCCCATTTTTCCGTGCGGGAGGGAAGAACATCGCCAAAAAAATACCCACAGCa                                                                                                                                                                              >  2:270547/1‑238 (MQ=255)
                                                                                  aTGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAACCCGCCACTTTTACCCCCTTTTTCGGTGCGGCCGCGAAGTACATCGCTCAAAAAAATCCCCCGGCACCCAAGAATTTTGa                                                                                                                                                               >  1:312381/1‑239 (MQ=255)
                                                                                     ttGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGTTTCGGGGAAGAACAAACCCCCCACTTTTACCCGATTTTTCCGTCCGGCCGTGAAGAACATCGCTAAAGAAAATCCCCACGCACCCACACAttttttgtt                                                                                                                                                             >  2:505223/1‑235 (MQ=255)
                                                                                      tccccACCACACCAAAGAGAACAATGAGGGTGGAAATGGAGGCCGCCAGCACGGGGTGTACCGAAAACAGCAACGGTAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTa                                                                                                                                                           <  1:688805/237‑1 (MQ=255)
                                                                                       gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAATGGTGCGGGGATCAGCCAGCCCGCCACTTTTTCCCCATTTTTCCGGGCGGCCGGGGAGTACATCGCTAAAAAAAATCCCACGGCCCCAAAAAATTTTTAGTTa                                                                                                                                                           >  1:274649/1‑238 (MQ=255)
                                                                                       gCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAACCCGCCCCTTTTACCCGCTTTTTCCGGGCGGCAGTGGAGGTCAACCCCAAAAAAATGCCCACCGCCCCAAAGATTTTTTAGGTa                                                                                                                                                           >  1:320884/1‑238 (MQ=255)
                                                                                           accacaccaAAGAGAGACAGCAGGGTCGAAAGTGAGGGCGCGACAGGGGTGATGACCGCAAAAAGAAACGGTGAAGTAACGGGAAAGGAAAATGCGAGGGGTGCGATATTCCCCACCAGGAGGGTGGTTAGGGGTGCGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg                                                                                                                                                      <  2:368161/239‑1 (MQ=255)
                                                                                           accacaacaAAGAGGAACATTACGAGCGAAATGGGGGCCGCCAGCCCGGCGTGTACCGAAAACACAAACGGTGAAATGAAAAGGAAGGAAAATCCCAGCGGTTCGGTAACCCCAACCAGGAGGGCGGTTAGGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg                                                                                                                                                      <  2:1314339/239‑1 (MQ=255)
                                                                                               caccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAACCCCCCACCTTTTACCGATTTTTCCGTGCGGCCGTGGAGTACATCCCCAAAGAAAATCCCCCCGCACCAAAGAATTTTTTGTTAACGTGCAgg                                                                                                                                                 >  1:1200607/1‑239 (MQ=255)
                                                                                                          tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATTAGCCAGGCCGCCCCCTTTACCCCATTTTCCCGTGCGGCCGTGGAGGACAACGCTAAAGAAATGCCCCCGGGACCCAAACTTTTTGAGTTACCGCGCAAGGCGAAAAccc                                                                                                                                       >  2:1165510/1‑237 (MQ=255)
                                                                                                            cATTGCGGGCGAAGTGGAGGCCGCCAGCACGGCGTGTACCGAAAACAGCAACGGTGAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGCGGTAAAGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCtt                                                                                                                                   <  1:1118299/239‑1 (MQ=255)
                                                                                                                 cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGGGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCggg                                                                                                                              >  1:432679/1‑239 (MQ=255)
                                                                                                                   gTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGAACAGCAAGCCCCCCCCCTTTTACCCCTTTTTCGGGGCGGCAGTGGAGAACAACCCCAAAAAAAATCCCCCGGCCCCCAAGAATTTTTAGTTTCCCTGGCGGGCAAAAACCCCCTCCCGGGa                                                                                                                             >  2:1343247/1‑239 (MQ=255)
                                                                                                                    tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCGCCACTTTTACCCGCTTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGACCCCGGCCCCCAAGAATTTTGTGTTACCGTGCAGGGCAAAAACGCCCTCCGGGaac                                                                                                                           >  2:914463/1‑239 (MQ=255)
                                                                                                                    tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGATGGGGGGGCCAGCAAACCCCCCCCCTTTTCCCGAATTTCCCGGGGCGGCGGGGAAAACACCGCACAAACAAATCGCACCGTACAAAAAATTTTTGAGGTACACTGGATGGGCAAAACGGCTTCCCGGGaa                                                                                                                            >  2:49826/1‑239 (MQ=255)
                                                                                                                     cGACATGGAGGCCGCCACAGCGGCGTGCACCGAAAACAGCAACGTTGAAATGAACAGGAAGTAAAATTCCAGCGGTGCGGTAATCCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaac                                                                                                                           <  2:274649/238‑1 (MQ=255)
                                                                                                                     cGACAGTGAGGCCGCCAGCAGGGCGTGTACCACAAACAACAACGTAAAAATAAACGGGAAGGTAAGTCCCGGCGTTTCGTTAATCCCCCCCAGGATGGCGGTAAGGGTGGCGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaac                                                                                                                           <  2:284471/238‑1 (MQ=255)
                                                                                                                                  gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGAACAACAAAACCCCCCCTTTTTTCCCCTTTTTCCGGGGGGGCGTGGAGAGCAACGCTAAAAAAAAGCCCCCCGCCCCAAAGACTTTTTTGTTTACCGGCCGGGCAAAAACCCCCTTCCGGGACAAAGGAtttagagg                                                                                                              >  2:834755/1‑234 (MQ=255)
                                                                                                                                        agggggtgtCCCGCAAACAGAAACGGTAAAATGAACAGGAAGGTAAATCCCGGGGGTTCGGTAATTCCCCCCAGCAGGGCGGTTAAGGTTGCGGAAAACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTc                                                                                                          <  1:338102/231‑1 (MQ=255)
                                                                                                                                        acgtcgttcaCCACAAACAGCATCTGTAAAATGAACAGGAAGATAAATCCCAGCGGTCCGGTAATCCCCCCCAGCATGGCGTTTGAGGTTGGTGAAACCAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTc                                                                                                          <  2:159476/227‑1 (MQ=255)
                                                                                                                                                  ccGCAAACAACAGGGGAAAAAGGAACGGAAGGGAAAATTCCAGCGGTGCGGAATTCCCCACCAGCGGGGGGTTTAAGTTGGGGGGAACAAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTg                                                                                             <  1:444499/239‑1 (MQ=255)
                                                                                                                                                    caaaacgccAACGGTGAAAGGAACAGGAAGAGAAACTCAAGCGTTGGGCTAATTCCCACCAGCAGGGCGGTTAGGGTGGGGGAAAACAGCAATCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTCCCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                                           <  2:925732/231‑1 (MQ=255)
                                                                                                                                                    acAAACAGCAAGGGAAAAAAAAGAAGGAGGGTAAATCCCGGCGGTGAGGTAATCCCCCCGAGCAGGGTGTTTAAGGTGGCGGGAAACAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                                           <  2:234251/238‑1 (MQ=255)
                                                                                                                                                        aCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCCCCTTTACCCCATTTTCCCGTGGGGCCGGGAAGTACATCCCTAAAGAAAATTCCACCGCCCCCAAAATTTTTTTGTTTACCTGCAGGGCAAAAACGCCCTCCCGGGACAACGATTTCAACGGCCCCGCGGCCCAACTGAActct                                                                                       >  2:396004/1‑239 (MQ=255)
                                                                                                                                                             aaCGGAGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCAGGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAg                                                                                  <  1:1265340/239‑1 (MQ=255)
                                                                                                                                                                 gagaaaaaaaCGGGAAGGAAAACCCCGGCGTTTCGTAATTCCCCACCAGCAGGGGGGTTAGGGTGGGGGGAACCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                               <  2:441208/230‑1 (MQ=255)
                                                                                                                                                                     aaggacagtgAAGTTAATTCCCGGCGGTCCGCAAACCCCCACCGCCGTGGGGTTAAAGGTTGGGGAAAGCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                          <  2:1141593/230‑1 (MQ=255)
                                                                                                                                                                     aaggAACAGGAAGGAAAATCCCAGCGGTTCGTTAATTCTCACCAGCAGGGCGTTTAGGGTGGCGGGAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                          <  1:1006323/236‑1 (MQ=255)
                                                                                                                                                                     aagaaaCAGGAAGGAAAATTCCAGCGGTTCGTTAATCCCCCCAAGCAGGGCGGTTAAGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                          <  1:914463/235‑1 (MQ=255)
                                                                                                                                                                         aaCAGGAAGGTAAATTCCAGCGGTGCGGTACTCCCCACCAGCAGGGCGGTTAAGGTTGCGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                        <  1:1085209/237‑1 (MQ=255)
                                                                                                                                                                         aaCAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCATGGGGGTTAAGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGccc                                                                       <  1:165303/238‑1 (MQ=255)
                                                                                                                                                                           cAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGGGCGGCAGTGAAGTACATCGCCAAAGAAAAGCCCACGGCACCAAAGATTTTTTAGTTACCGTGCCGGGCAAAAACGCCCTCCCGGGACAACGATTTCAACGGCTCGGGGCCCAAAACTAACTCTTTTCGATGCTCGGccccgc                                                                    >  2:682382/1‑236 (MQ=255)
                                                                                                                                                                                                               ccccggCATGGGGGTTAGGGTTGGGGAAATCAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                                   <  1:505223/234‑1 (MQ=255)
                                                                                                                                                                                                               ccccggCATGGGGGTTAGGGTGGGGGGATACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTTCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                                   <  2:1219011/234‑1 (MQ=255)
                                                                                                                                                                                                               caccaGCATGGCGGTTAGGGTGGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                                   <  1:623841/239‑1 (MQ=255)
                                                                                                                                                                                                               cacaaGCGTGGGGGTGAAGGTTGCGAGAAACAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                                   <  2:1130817/239‑1 (MQ=255)
                                                                                                                                                                                                               cacaaGCATGGCGGTTAAGGTTGCGGGAAACAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                                   <  1:855493/239‑1 (MQ=255)
                                                                                                                                                                                                               cacaaGCAGGGCGATTAGAGTCGCGGGAATCAGCAAGGCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                                   <  1:1308047/239‑1 (MQ=255)
                                                                                                                                                                                                                 ccaGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg                                                                                                                                                      >  1:884646/1‑124 (MQ=255)
                                                                                                                                                                                                                      ggtgcggggtAAGGTGGGGGGAAACAGAAAGACCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAg                              <  2:747996/229‑1 (MQ=255)
                                                                                                                                                                                                                          cGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                         <  1:270547/238‑1 (MQ=255)
                                                                                                                                                                                                                          cGGTTAAGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                         <  1:857521/238‑1 (MQ=255)
                                                                                                                                                                                                                            gtgtaagggTTGGGAGAATAAGAACGCCCGCCACTTTTACCCGATTTGCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                         <  2:953200/231‑1 (MQ=255)
                                                                                                                                                                                                                             ggtAAGGGTGCGGGAAACCGCCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAAc                        <  2:320884/236‑1 (MQ=255)
                                                                                                                                                                                                                                 aaGGTGGCGGCAAACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTcc                   <  2:1288243/239‑1 (MQ=255)
                                                                                                                                                                                                                                   ggTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCACCAGCTGCCGGACCAAAGATAAACTGTCCGt                 >  1:564674/1‑239 (MQ=255)
                                                                                                                                                                                                                                    gTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt                 <  1:1343247/238‑1 (MQ=255)
                                                                                                                                                                                                                                      ttGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt                 <  1:886557/237‑1 (MQ=255)
                                                                                                                                                                                                                                         cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAACTGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACGACCTTCCGGGAACAACGATTTCACCGGCTCGGGGCTCAAACTGACCGCTTGTCGCTGCTGCGCCCAGGCCTTCTGAACTCCGCATTCAACACCTGCCGGAGCAACGAAAAACTCTGCGTAGATg            >  2:721404/1‑239 (MQ=255)
                                                                                                                                                                                                                                         cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCACAAGCTGCCGGACCCAAGATAAACTGTCCGTAGATg            >  2:125505/1‑239 (MQ=255)
                                                                                                                                                                                                                                            aaaTTAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa          <  1:1002779/237‑1 (MQ=255)
                                                                                                                                                                                                                                              aagcAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg         <  1:671317/235‑1 (MQ=255)
                                                                                                                                                                                                                                               tCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                                   >  1:675449/1‑178 (MQ=255)
                                                                                                                                                                                                                                               accAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAgtg       <  1:1217169/237‑1 (MQ=255)
                                                                                                                                                                                                                                                     cAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCAAGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTa  >  2:1251530/1‑239 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
TGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTA  >  NC_000913/3862438‑3862912

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