Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,862,670 T→A 28.7% pseudogene (934/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,6700TA28.7% 67.2 / 13.3 72pseudogene (934/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 T (24/20);  new base A (7/10);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

ATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAA  >  NC_000913/3862435‑3862904
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCTCACCAGCATGGCGGTTTAGGTTTGGGGAATCAg                                                                                                                                                                                                                                                   >  2:1006323/1‑239 (MQ=255)
aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAAtta                                                                                                                                                                                                                                                     >  1:117052/1‑236 (MQ=255)
   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)
                                                                                                                                                                                                                                                  tCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                           >  1:675449/1‑178 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
ATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAA  >  NC_000913/3862435‑3862904

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: