Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,862,661 T→G 20.7% pseudogene (942/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,6610TG20.7% 93.5 / 3.6 67pseudogene (942/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 (34/17);  new base G (4/12);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 4.39e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAA  >  NC_000913/3862424‑3862888
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTg                                                                                                                                                                                                                                                  <  1:1213083/239‑1 (MQ=255)
     tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCgggg                                                                                                                                                                                                                                            >  1:207405/1‑238 (MQ=255)
     tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGa                                                                                                                                                                                                                                            <  2:684848/239‑1 (MQ=255)
     tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGCTTAGAGTGCgggg                                                                                                                                                                                                                                            >  1:986196/1‑238 (MQ=255)
          aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATTAAGGTTGGGGAATCac                                                                                                                                                                                                                                       >  2:248536/1‑238 (MQ=255)
          aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCa                                                                                                                                                                                                                                        >  1:185378/1‑238 (MQ=255)
                  cAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGAATCacaaaccc                                                                                                                                                                                                                                 >  1:502906/1‑231 (MQ=255)
                         atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGTTGCGGGAATCAACAAGccccccccctt                                                                                                                                                                                                                    >  2:1081795/1‑234 (MQ=255)
                                cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGCGGGGATCAGCAAACCCCCCACTTTTTCCCgg                                                                                                                                                                                                          >  1:1097435/1‑238 (MQ=255)
                                cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGAAAATTCCAGCGGTTCGGTAACTCCCACCAGCATGGCGGGTAAGCGTTGGGGGaacaacaacccacccacattcacac                                                                                                                                                                                                                    >  2:244337/1‑214 (MQ=255)
                                    gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGTGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGTTGCGGGAATCAGCAAACCCGCCACtttttacccttttt                                                                                                                                                                                                       >  1:109894/1‑229 (MQ=255)
                                    gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGTGGCGGGGATCAGCAACACCCCCCCTTTTTCCCCAttt                                                                                                                                                                                                       >  1:773953/1‑238 (MQ=255)
                                    gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCACATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAACAACGGAGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGAAAACTTGCAGGGACAAACAAGCTCGCCCCATTTAACCCAtttt                                                                                                                                                                                                      >  2:1040278/1‑239 (MQ=255)
                                          gaccttgggaacccAGTTGGGCGGAAACACGTAGTAAACCAGCCTCCCGCCCATGTGGCCCACCACCCAAAAGAGAAAAATCACGGTCGAAAGTGAGGCCACGAGCCGGGCGTGTCCCCCAAACAGCAACGTGAAAAAGAAGAGGAAGGAAAATCCCGGCGGTTGGGTAATCCCCCCCAGCAGGGCGGTTAAGGTGGGGGAAATAAGGAACCCCGCCACTTTTACCCGATTTTCCGGTg                                                                                                                                                                                                <  2:366749/227‑1 (MQ=255)
                                           aaCGTCGGGATCCGTTTTTGCAGAAAACCCTAGCCAATCAGCCCCCCCCCCATGTCCCCCACCACACCAAAGAAATACATGCCGGTCGACATGGAGGCCGCCAGCACGGGGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGTTAAATCCCGGCGGTTGGGTAATCCCCCCAAGCTGGGGGGTTAAGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTg                                                                                                                                                                                                <  2:835676/238‑1 (MQ=255)
                                            cgggcGGGATCCAGTTTCGCGAAAAACCCTACCAAACCAGCCCCCCCCCCATTGTCCCCCCCACACAAAAGAGATACATTACGGTCGACATGGAGGCCGCCGGCACGGCGTGACCCAAAAACAGCAACGGTGAAATGAACGGGAAGGAAAACCCCGGCGGTTCGGTAATCCCCACCAGAAGGGGGGTTAGGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTg                                                                                                                                                                                                <  2:639637/233‑1 (MQ=255)
                                                  ggATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGGGGGGAATAACGAAAGCGCCCCCTTTTACCAATTTTCTGGTgggggcgag                                                                                                                                                                                        >  2:1249681/1‑234 (MQ=255)
                                                       cattttgGCGGAAAAACATCGTAAACCACACCCCCCCCAATGTCCCCCACCACACAAAAGAAATACATCGGGGTCGACATTGAGCCCGCCAGCAGGGCGTGTCCCGCAAACAGCAGCGGTAAAAGGAAAGGGAAGGTAAATCCCAGCGGTTCGGTAACCCCCCCCAGCAGGGCGGTTAAGGTTGCGGAACTAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                                                                                                                  <  2:207405/236‑1 (MQ=255)
                                                       cagtttggtgaaaaaaCCTAACCAACCAGCCCCCCGCCCATGTCCCCCCCCACCCAAAATAGAAACATTCGGGTGTAAGTGGCGGCCACCAGCAGGGCGTGTGCCGAAAACAGAAAGGGTAAAAGGAACGGGATGATAAGTCCCAGCGGTGGGGAAATCCCCCCCAGTTTGCGGTTTAAGTTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                                                                                                                  <  2:524086/227‑1 (MQ=255)
                                                        aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGATCAGCAAGCCCCCCACTTTTACCCCATTTTTCTGTGCGGCAGTGAAGTAc                                                                                                                                                                                 >  1:762962/1‑239 (MQ=255)
                                                          ttttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGGGGGGAATAGCAAACCCGCCCCCTTTTCCCGATTTTCCGGTgggggcggggaggcg                                                                                                                                                                                  >  1:639637/1‑226 (MQ=255)
                                                          tttgggcaaaaaaCCTATACACTCAGCCCCCCGCCCAGGTCCCCCCCCACCCAAAAGAGAAACATTACGCGCGAAATGGAGGCCGCCAGCAGGGGGTGTACCGAAAACAGAAACGGGAAAAGGAACAGGAGGTTAAATCCCGGCGGTGCGGTAATCCCCACCAGCATGGGGGTTAGGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt                                                                                                                                                                               <  2:1041887/230‑1 (MQ=255)
                                                           tttGGGATAAACCCGGGTCAACCCGACCCCCCCCCTTTTGCCCCCCCCCACAAAGGACATACTTTAGGGTGGACATGGAGGCCGCCAGCCCGGCGGGACCCCAAAACAGACACGGAGAAATGATAAGGAGGGAAAGTCCAAGCGGAGCGTTACTTCCCCCCGGCATGGGGTTTAGGGTGGGGGGAAGAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAc                                                                                                                                                                                 <  2:874743/236‑1 (MQ=255)
                                                            tgcgggaaaaaCCTAGTAAACCAGCCCTCCGCCCATGTCGCCCCCCACACAAAAAAGATACATTACGGTGGAAATTGAGGCCGCCGGCCCGGCGTGAACCGAAACCACAAAGGGAGAAAGAAACAGGAAGTAAAATCCCAGCGGTTCGGTAATCCCCACCAGCATGGCGGTTAAGGTTGGGGAAAACAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCg                                                                                                                                                                             <  2:1097435/232‑1 (MQ=255)
                                                                    aaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGTGGGGATCAGCAAAACCGCCCCTTTTTTCCGATTTTTCCGTGCGGCGGTGGAGGACATCGCTaaaaaa                                                                                                                                                                     >  1:231935/1‑239 (MQ=255)
                                                                            tCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTTGGGGGAACAACAAAGCCCCCCCTTTTATCCCGTTTTTCCGGTGGGGCGGGAAAGAACTCGCCAAAGAAAAGcccc                                                                                                                                                              >  2:316051/1‑237 (MQ=255)
                                                                                tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGGGGGGATCAGCAAGCCCCCCCCCTTTTCCCCATTTTCCGGGGCGGGAGGGGAGTACAACCCCTAAAAAAAGCCCCCGGc                                                                                                                                                          >  1:136869/1‑238 (MQ=255)
                                                                                tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGCTAAGGATTCGGGGATCCGCACACCCCCCCCCTTTACCCCAATTTTCCGTGCGGCGGTGAATGACATCGCCAACAAAATGctcccgac                                                                                                                                                          >  2:517063/1‑231 (MQ=255)
                                                                                 cAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCCCCCCTTTTAACCGCTTTTCCGGGGGGGGAGTGAAGGACATCGCTAAAAAAATGCCCACGGCa                                                                                                                                                         >  1:968735/1‑238 (MQ=255)
                                                                                    aCCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTTCGGGGAACAGCAAACCCCCCCCCTTTTACCCATTTTTCGGTGCGGCAGTGGAGTACATCGcctgtctcttttacccatcctcgg                                                                                                                                                     >  1:366214/1‑216 (MQ=255)
                                                                                         cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGAACAGCAAACCCGCCCCTTTTTCCCCGTTTTTCGGTGCGGGCGGGGAGTACACCGCCAAAAAAAAGCCCCCCGGCCCCAAGAtt                                                                                                                                                >  1:535343/1‑239 (MQ=255)
                                                                                              gaaaaagaacaagccaaaaaGAGAAAAATGACGGTCGAAATGGAGGCCGCCAGCAGGCGGTGCCCCAAAAACAGAAACGGTAAAAAAAACAGAAGGAAAAATCCCGGCGGTTGGGTAATTCCCACCAGCATGGCGTTTAAGGTTGGGGAAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTg                                                                                                                                            <  2:890572/221‑1 (MQ=255)
                                                                                               tGTTGCCCACCACACCAAAGAAATACATTACGGACGACATTGAGGCCGCCAGCACGGCGTGTCCCGCAAACAGCAACGGTGAAATAAACAGGAAGGAAAATTCCAGCGGTTCGGTAATCCCCACCAGCATGGCGTTTAAGGTTGCGGGAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAg                                                                                                                                         <  1:1037355/239‑1 (MQ=255)
                                                                                               ggtcccccACCACCCAAAAGAGATACATTACGGTCGACATGGAGGCCGCCAGCACGGGGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGATAAATTCCAGCGGTTCGGTAACTCCCACCAGCAGGGCGGTAAGGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTg                                                                                                                                            <  1:1427858/232‑1 (MQ=255)
                                                                                                       accacaccaAGGAGCAACATTGCGGTAGGAATGGGGGCCGCCGGCACGGCGGGTCCCGCAAACAGCAACGGTAAAATGAAAAGGAAGGTAAATCCAAGGGGTTCGGTAATCCCCACCAGCAGGGCGTTTAGGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg                                                                                                                                <  2:1098257/239‑1 (MQ=255)
                                                                                                         cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTGCGGGGAACAGCAAGCCCCCCACTTTTACCCGATTTTTCGGGGGGGCGGGGAAGGAAAACGCTAAAAAAATGCCCACCGCACCACAGAATTTTGAGGTTACCGTga                                                                                                                               >  1:1151775/1‑238 (MQ=255)
                                                                                                            accaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCGCGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCCCCAGCAGGGCGGTTAAGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCAc                                                                                                                                                             <  1:44161/207‑1 (MQ=255)
                                                                                                              caAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGATCCGCAAACCCGCCCCTTTTTCCCGCTTTTTCCGGGCGGCAGGGAGATACATCGCCAAAGAAACGGCCCCCGCACCCAAACTATTTTAGGTACCCTGCAgggga                                                                                                                        >  2:1030318/1‑237 (MQ=255)
                                                                                                                     aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGGAGTACATCGCTAAAGAAATGCCCACGGCACCCAAGATTTTTGAGTTACCGTGCACGGGAAAACCg                                                                                                                >  2:770245/1‑238 (MQ=255)
                                                                                                                      tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAACCCGCCCCTTTTTCCCCGTTTTTCCGGGCGGGAGGGGAGGAAACCGCCCAAAAAAAGCCCCCCGGCCCCAAGATTTTTTTGTTTCCCTGGCAGGGCAAAAccc                                                                                                                >  1:696034/1‑237 (MQ=255)
                                                                                                                      tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGCGGGGATCAACAAACCCCCCACCTTTTCCCCCTTTTTCCGGGGGGCAGTGGAGAACAACGCTAAAGAAAAACCCACCACACCAAAAAATTTTTGGTTAACGGGcccgtttcttttaca                                                                                                                 >  1:890572/1‑225 (MQ=255)
                                                                                                                        cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGAACAGCAAACCCCCCCCTTTTACCCCATTTTTCGGGGCGGGAGTGAAAGACAACGCCAAAAAAATGACCACGGCACCAAAAATTTTTTAGTTACCCTGCACGGGCAAAAcccct                                                                                                              >  2:537825/1‑235 (MQ=255)
                                                                                                                             cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGCCGGGAACAACCAAGCCCCCCCCTTTATCCCATTTTTCGGGGCGGCGCGGAAAGAACAACCTCAAAAAAATCCCACCGGCCCCAAAAATTTTTTAGTTTCCGGGAGGGGACAAACCCCCTTcccg                                                                                                        >  2:491274/1‑237 (MQ=255)
                                                                                                                             cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATGAAGGTTGGGGAATTAGGAAACCCCCCACCTTTTCCCGATTTTTCCGGGGGGGAGTGGAGTAAAACGCTAAAAAAAAACCCACAGCACCAAAGAATTTTGAGATAACGTGGCAGGCAAAAACGCCTTCCggg                                                                                                       >  2:1382028/1‑239 (MQ=255)
                                                                                                                               gTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGATCAGCAAAGCCCCCCCTTTTTACCCATTTTTCCGTGCGGCCGTGGAGGAAAACGCCAAAAAAAAGCCCCCCGCCCCCAAGAATTTTTAGGTACCGGGCAGGGCCAAAACCCCCTCCCGGGa                                                                                                      >  1:1048976/1‑239 (MQ=255)
                                                                                                                                  gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGAATCAACAAAACCGCCCCCTTTTCCCCCATTTTCGGTGCGGGCGTGAAGGACAACCCTAAAAAAAAACCCACCGCACCCCAAATTTTTTAAGTAACGTTGAGGGCAAAACCCCCTCCCGGGaaaca                                                                                                  >  2:1424731/1‑236 (MQ=255)
                                                                                                                                  gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGTTGCGGGGAACAGCAAGACCCCCACCTTTTCCCGATTTTTCGGTGCGGCAGGGAAGTAACACGCTAAAGAAAAGCCCCCCGCACCAAAGATTTTTTAGGTAACGGGCCAGGGAAAACCCCCCTCCCGGaacaa                                                                                                  >  2:234796/1‑239 (MQ=255)
                                                                                                                                              gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGAACAGCAAGCCCCCCACATTTAACCGATTTTTCGGTGGGGCAGGGAAAGAAAATGCCAAAAAAAAGCCCCAGGCCCCAAAAATTTTTTGGTTCCCGGCCGGGGCAAAACCCCCTTCCGGGAAAAACAAttctaccg                                                                                       >  2:1072488/1‑232 (MQ=255)
                                                                                                                                              gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGTGGGAATCAGAAAAGCCCCCCCCTTTACCCCCATTTTCCGGGCGGGGGGGAAGGACAACGCCAAAAAAAAGCCCCCCGCACCCAAAATTTTTTGGTTATCCTGTATGGCACAAAACCCCTCCCGGGagaaaaatttataccg                                                                                       >  2:402350/1‑224 (MQ=255)
                                                                                                                                                                caaaaCAGCAAGGGTGAAATGAACAGGAAGGAAAATCCCGGCGGTTCGGAAATTCCCACCAGCATGGCGGTTAGGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                    <  1:994202/237‑1 (MQ=255)
                                                                                                                                                                aaaaaCAGCAGCGGTGAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCACCAGCATGGGGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                    <  1:1438721/237‑1 (MQ=255)
                                                                                                                                                                 cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTTCGGGGATTAACAAACCCCCCCCCTTTTCCCGGTTTTTCGGTGCGGCAGGGAAGTAAATCGCCACAAAAAAACCCCCAGCACCAAACACTTTTTAGTTACCCCTCAGAGCAAAACCGCCTGCCCGGGACAACCAATTTAGCGGAGCGGCGCtacacactca                                                                    >  2:1386523/1‑229 (MQ=255)
                                                                                                                                                                    aaaGCAACGGTAAAAGGAACATGAAGGTAATTCCCAGCGGTCCGGAATCCCCCCCCAGCGTGGGGGTTAAGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAActct                                                                <  2:136869/237‑1 (MQ=255)
                                                                                                                                                                         aaCGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTTCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCCCCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCCCCAACTGAACTCTTGccg                                                           >  1:309975/1‑237 (MQ=255)
                                                                                                                                                                           cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGACCCCCCCCTTTACCCGATTTTCCGGTGCGGCCGGGAAGTACATCGCTCAAGAAATTCCCACGGCACCCAAGATTTTTTAGGTAACCCGCAGGGGAAAAACGCCTTCCGGGAAAAACGATTTCCACGGCCCGGAGCCCAAACAGaaacacttcagat                                                         >  2:446498/1‑228 (MQ=255)
                                                                                                                                                                               gaaaagacagggagtgaaatctCCAGCGGTAGGGTAATCCCCCCCGCCGGGGGGTTAAGGTGTGCGGAAACAAAAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                      <  2:968735/217‑1 (MQ=255)
                                                                                                                                                                                     aaaagaaaGGAAAATCCCGGCGGTCGGGAAATTCCCCCCAGCATGGCGTTAAGGGGTGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                 <  1:530942/231‑1 (MQ=255)
                                                                                                                                                                                     aaaagaaaGGAAAATACCAGCGTTTGGTAAATCCCCACCAGCAGGGGGGTTAGGGTGGGGGAAAACAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                 <  2:1419186/231‑1 (MQ=255)
                                                                                                                                                                                     aaaaGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCACCAGCATGGCGATTAGGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                 <  1:770245/234‑1 (MQ=255)
                                                                                                                                                                                      agaggatgtaaaaTCCCAGCGGTACGGTACTCCCCACCGGCATGGCGGTTAAGGTTGCGGAAAGCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                              <  1:780352/229‑1 (MQ=255)
                                                                                                                                                                                      aCAGGAAGTTAAAGTCCAGCGGTTCGTTAATCCCCCCCAGCATGGCGGTTAAGGGTGCGGGAATAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                              <  1:650117/239‑1 (MQ=255)
                                                                                                                                                                                          gAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAAGTTGCGGGAATCAGCAAACCCCCCCCTTTTTCCCCATTTTCCGGTGCGGCCGTGAAGTACATCGCTAAAAAAATTCCCACGGCACCAAAAATTTTTTAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACCACGATTTCAGCGGCTCCGCGCTCCAACTTAAATCCTGCAGATGCTCTGccccctac                                           >  2:530942/1‑233 (MQ=255)
                                                                                                                                                                                                         ggtttggttAATCCCCCCCAGCAGGGCGGTTAGGGTTGCGGAAAACAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCt                            <  1:316051/230‑1 (MQ=255)
                                                                                                                                                                                                                   actccaccaagaggggCGGTAAGGGTGGCGGAAACAAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                 <  2:696034/225‑1 (MQ=255)
                                                                                                                                                                                                                        caccaGAATGGCGGTAAAGGTGGGGGGAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc            <  1:9660/239‑1 (MQ=255)
                                                                                                                                                                                                                        cacaaGGATGGGGGTTACGCGTGTGGGAATGCGAAAGCCCGCCACTTTTACCCGATTTGCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATTCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACATCTGCCGGAc            <  2:107189/239‑1 (MQ=255)
                                                                                                                                                                                                                                    gggTTAGGGTTGGGAGAAGCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa  <  1:84386/237‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAA  >  NC_000913/3862424‑3862888

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