Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,862,956 G→A 19.3% pseudogene (648/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,9560GA19.3% 92.3 / 2.2 53pseudogene (648/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base G (24/16);  new base A (8/2);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGAT  >  NC_000913/3862732‑3863183
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
aTGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGCAGTTGTGTAAACCGGTTGGGATCAGAATACGTTCTAGTACGGGTTAAACCCAAAc                                                                                                                                                                                                                                      >  1:304321/1‑227 (MQ=255)
aTGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCCGCGCCCAGGACATCTGAATTCCGCCTTCCACAGCTTCCGGACCAACGATAAACTTTCCGTAGCTGAGTTGGGGTAAACCGGTTGGGACAGAAAAACGCTCGCGGAAAGGGaaaaccaaaacc                                                                                                                                                                                                                                     >  2:19586/1‑228 (MQ=255)
       cGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGGGTAACCCGGTTGGGATCAAAATACGATCGCGGAAGGAGGAAACCCAAACCCCAAgcg                                                                                                                                                                                                                             >  1:258530/1‑226 (MQ=255)
        ggCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACCGAACTCCTGCAGATACTGCGCCCAGTACACCTCAATGCCGCCTTCAACAGCTGCCGGCACAAAGAAAAACATTCCGAAGACGAAGGGGGTTAAACCGGGGGGGACCCAGAAAAGTTCGAGGAACAGGGAAACCAAAAACCCAAGAGc                                                                                                                                                                                                                            >  2:250894/1‑228 (MQ=255)
                    ttttGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGACCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAAAACGTCACAGGAAAGTGCAAACCACAACCCCAAGGCGACCCGGCGaaac                                                                                                                                                                                                                >  1:443572/1‑226 (MQ=255)
                        gAGTTACCGGGCAGGGCAAACCCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCACAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAgg                                                                                                                                                                                                            <  1:41318/228‑1 (MQ=255)
                             aCCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTCGATGAAGTGGCGTAAACCGGTTGGGATCAGAAAACGTTCGGGGAAGGGGCAAACCCAAACCCCAAGGGCACACCCCGAACCCCGGAACGc                                                                                                                                                                                                       >  2:200850/1‑228 (MQ=255)
                             aCCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGCTGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGGGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACga                                                                                                                                                                                                       >  1:78242/1‑227 (MQ=255)
                             aCCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGGGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACg                                                                                                                                                                                                        >  1:650436/1‑227 (MQ=255)
                                      ggCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTAAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTTGGATCAGCATACGTTCGAGGAAGGTGTAAACCCAGACCCACAAGGCACCCGCCGAACTAAGGAAAGCTTCCagctg                                                                                                                                                                                              >  1:691349/1‑225 (MQ=255)
                                       gCAAAACCGCCTACCGGAAAAAACGATTCCAGCGGCTCGGCGCTCAAACGGACCCCTTGCAAATGCTGCGCCCAGTACATCTGAATGCCGCCTGCAACGGCTGCCGGACAAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATCCGTTCGAGGAAGGTTTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCa                                                                                                                                                                                                  <  1:660790/223‑1 (MQ=255)
                                           aaCCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGGGTAAACCCAGACCCCAAGGGCACCCGCCGAACCCCGGAAGGCTTTCaagaattca                                                                                                                                                                                          >  1:539806/1‑219 (MQ=255)
                                                ccTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCCAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAATCCACACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATcccc                                                                                                                                                                                    >  1:593727/1‑228 (MQ=255)
                                                ccTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATcccc                                                                                                                                                                                    >  1:484681/1‑228 (MQ=255)
                                                     cGGGAACAACGATTTCAGCGGCTCGGCGGTCAAACTGAACTCTTGCAGATGCGGCGCCCAGGACATCGGAATGCCGCCTTCAACCGCTGCCGGACCAAAGAAAAACTGTCCGTAGATGAAGTCGTGTAAACCGGTTGGGATCAGAAACAGGTCGAGGAAGGGGGAAACCAAAACCCCAAATGCACCCCCCAAATCCAGGAAGCCGTGCAAAGAGTCAACCCCCAATTg                                                                                                                                                                               >  2:666179/1‑228 (MQ=255)
                                                                                                tGCAGATGTCGCGCCCAGTAAATCGGAAGGCCGCCTTCAACAGCCGCCGGCCCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATTCGGTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGatc                                                                                                                                   <  1:83245/228‑1 (MQ=255)
                                                                                                       gctgcGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAAACCCCAAGTGCACCCGCCGAACTCAGGAAGGCTTGCAGGACTCCAATCCCCATTTGTACTTTTGGCCGGCCACGCACCGGCACCCAGGCCCACGGGCTCACCCCCa                                                                                                                            >  2:240749/1‑228 (MQ=255)
                                                                                                         tgcGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGGGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAAATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCACCCAGGCCCAGGGGATAAACCCCAga                                                                                                                          >  2:140695/1‑227 (MQ=255)
                                                                                                          gcgcCATGTACTTCGGAATGCCACCTGCAACAGGTGCGGGACCAAAAGAAAATGGTCCGTGAATGCATGGGTTTAACCCGGTTGGGCTCAGAATCCGTCCGAGGAAGGTGTAAACCCAGACCCCAAGTGCGCCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTCGGCCAGCCGAGCAACGTCACCCAGGCACAGTGGATCATCACCAGGa                                                                                                                         <  2:492142/228‑1 (MQ=255)
                                                                                                              ccAGTACATCGGAATGCCGCCTGCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTGGGGATCAGAATACGTTCGAGGAAGGGGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAgg                                                                                                                      <  1:240749/227‑1 (MQ=255)
                                                                                                                  tACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATa                                                                                                                 >  1:226492/1‑228 (MQ=255)
                                                                                                                       tctaaggccgacTTCAACAGTGGCCGGACAAAAGAAAAACGGTCCGAGAAGGAAGTGGTGTAAACCGGTTGGGATCAAAATACGTTCGAGGAAGGTGAAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATcac                                                                                                            <  2:539806/222‑1 (MQ=255)
                                                                                                                       cTGAATGCCCCTTCAAACAGCTGCCGTACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATcac                                                                                                            <  1:416901/228‑1 (MQ=255)
                                                                                                                            tGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGCGATTCAATCCCCATTTGTACTTTTGGCCAGCCCCGCAACGTCAGCCAGGACACGGGGGACATCACCAGTAAGGCGGTaaactttc                                                                                                          >  2:7287/1‑219 (MQ=255)
                                                                                                                            tGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTAAGCCAGGCCCAGGCGATAATCCCCAGGAAGGCGAAAAACacaacaa                                                                                                       >  2:143660/1‑227 (MQ=255)
                                                                                                                            tGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGGGCACCCGCCGAACGAAGGAACGCTTGCACAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTAAGCCAGGCACAGGGGATCATCAACAGGAAGACGAAAATcacacaaa                                                                                                       >  2:304323/1‑223 (MQ=255)
                                                                                                                             gccgccTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAAt                                                                                                               <  1:143660/219‑1 (MQ=255)
                                                                                                                             gccgccTTAAAAAGCCGCCGGACCAAAGATAAACTGGCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAgg                                                                                                                          <  1:7287/208‑1 (MQ=255)
                                                                                                                               cgccTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAAACCCCAAGGGCACCCGCCGAACGCAGGAACGCTTTCAGACATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTAAGCCAGGCACAGGGGAACATCACCAGGAAGGCGATAAACCCCACAAAAg                                                                                                    >  2:417431/1‑228 (MQ=255)
                                                                                                                                 ccTTAAAAAGGTGCCGGACCAAAGAGAAAGTGCCCGTAGATGAAGTGGTTTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGAc                                                                                                  <  1:140695/228‑1 (MQ=255)
                                                                                                                                    tCAACAGCTGCCGGACCAAATATAACCTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCCAGTGCACCCGCCGAACGCAGGAACGCTTGCCGAGAGTCAATCCCCATTTGTACTTTTGGCCACCCTATCACCGTCAGCCAAGCCCACCGGATCATCACCAGGACGGCGCTAAACACCAACAAAGACGt                                                                                                >  2:28628/1‑227 (MQ=255)
                                                                                                                                                                  tCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGCAAATGCCGAGAAAAACCGGCAg                                                                  >  1:144937/1‑228 (MQ=255)
                                                                                                                                                                  tCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAg                                                                  >  1:118168/1‑228 (MQ=255)
                                                                                                                                                                  tCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAAAGTCAGCCAGGACAATGGGATCATCACCAGGACGGAGATAATAACCACCTAAACAGGCCCCTGGAAAATACCGGGAAAAACCGGCAg                                                                  >  2:469407/1‑228 (MQ=255)
                                                                                                                                                                    cggaGATGAAGGGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGCGAAAAACCGGCAGtt                                                                >  1:454760/4‑228 (MQ=255)
                                                                                                                                                                              gtgggttAAACCTGTTGGGACCAGAACTCGTTGGAGGAAGGTGTAAACCCAGCCCCCAAGTGCCCCCGCCGAACGCAGGAACGCTTGCAGAGATTCAACCCCCATTTGTACCTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGa                                                       <  2:93611/221‑1 (MQ=255)
                                                                                                                                                                               tgtggtAAACGGGTTGGGTCCAGAATGGGTTCGAGGAAGTGGTAAACCACTACCCAAGGTGCCCCCGCCGAACGGAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCACTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGaa                                                      <  1:48086/223‑1 (MQ=255)
                                                                                                                                                                                ggccaaacccGGTTGGGCTCAGAATACGTTCGAGGAAGGGAGAACCCCAAACCCTAAGGTGACCGGCCGAAGGCAAGAACGCTTGCCGAGATTCAATCCCCGTTTGTACTGTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCAGCCCTTAAGGCGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACa                                                    <  2:593727/220‑1 (MQ=255)
                                                                                                                                                                                            ttGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGGCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCCCCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGGTTTTTATCGAACAAAcaggtaaaca                                        >  2:4722/1‑219 (MQ=255)
                                                                                                                                                                                                 gctCAGAATATATTGGAGAAGTGTTTAAACCCAGCCCCCACGTGCCCCCGCCGAACGCAGGGACGCTTGCAGAGACTCAATCCCCCTGTGTACTTTAGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAgcgc                                   <  2:258530/226‑1 (MQ=255)
                                                                                                                                                                                                             gTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCGGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCACCTGGACAATGCCGAGAAAAACCGGCGGTTTTTTATGGAACAGACGGTTATTCAGCGCCGTCCCAAGGc                        >  1:685818/1‑228 (MQ=255)
                                                                                                                                                                                                             gTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAAtcc                        >  1:78922/1‑226 (MQ=255)
                                                                                                                                                                                                                 gAGGAAGGGGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGa                    <  1:19586/228‑1 (MQ=255)
                                                                                                                                                                                                                    gAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGa                    <  1:304323/225‑1 (MQ=255)
                                                                                                                                                                                                                       ggTGTAAACCCAGACCCCAAGTGCACCCGCTGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAACTGCCGAGAAAAACCGCCAGTTTTTTATCGACCAGACGGTTATGCAGCGCCTTCACACTGCCGGAAATgca              >  1:419000/1‑226 (MQ=255)
                                                                                                                                                                                                                        gtgtAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACAACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGACAAGACGGTTATGAAACGCCGTCCCAAGACCGGAAATgcga             >  1:664218/1‑225 (MQ=255)
                                                                                                                                                                                                                           tAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTACCGAACAAACGGTTATGCAGCGCCGTACCAAATCCGGAAATAAAAAtt          >  2:465070/1‑228 (MQ=255)
                                                                                                                                                                                                                            aaaCCCAGACCCCAAGTGCCTCCGCCGAACGCAGGAACGGTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGATATGCAGCGCCGTCACAATGCCGGAAATGATAATTg         <  1:469407/228‑1 (MQ=255)
                                                                                                                                                                                                                            aaaCCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCCTTTGTACTTTTGGCCAGCAGAGCAACGTCAGCCAGGCACAGGGGTTCATCACCAGGAAGGCGATAATCACCCGAACAGCCGTCCCCTGGAAAATGCCCAGAAAAAACGGCCGTTTTTTATAGAACGAACGGTCATTCCCCGCCGTCACAATGACGAAAATTAGAAttc         >  2:34214/1‑227 (MQ=255)
                                                                                                                                                                                                                             aaCCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCGGGAAAAAGCCGGAAAAAACCGGCCGTTTTTTTCCGAACAGCCGGGTAGCAACAGCCGAACCAACGCCGGAAATgtaaaatga        >  2:429124/1‑220 (MQ=255)
                                                                                                                                                                                                                              ccccAGCCCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTgc        <  1:200850/226‑1 (MQ=255)
                                                                                                                                                                                                                              aCCCACACCCCAAGCCCACCCGGCGAACGCAGGAACGCTTGCAGAATTTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTgc        <  1:305673/227‑1 (MQ=255)
                                                                                                                                                                                                                              aCCCAACCCCCAGGGGCACCCGCCGAACGCAGGAACGGTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGGCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGG‑‑AATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTgcgc      <  1:250894/227‑1 (MQ=255)
                                                                                                                                                                                                                                  aGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGACAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACAACATAAAACGTCCCCTGGTAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTgacacag   >  1:269540/1‑222 (MQ=255)
                                                                                                                                                                                                                                   gACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGGACCGACGGGTAAGCAGCGCCGTCAAAAAGCCGGAAATGATAATTGCGCCGaa  >  2:110218/1‑227 (MQ=255)
                                                                                                                                                                                                                                   gACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCACCTGGAAAATGCGGAGAAAAACCGGCAGTTTTTTATCGAAACACCGGGTAGGCAGCGCCGCAACAATGCCGGGAATGATAATTAGGCCGAt  >  2:456038/1‑228 (MQ=255)
                                                                                                                                                                                                                                   gACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCAAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAAGGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTAACAAGGCCGGAAATGATAAtttgaccgaa  >  1:656309/1‑220 (MQ=255)
                                                                                                                                                                                                                                   gACCCCAAGTGCACACGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAAGGAATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTACAAAAGCCGAGAATGCTAATTGCGCCGAt  >  1:192181/1‑228 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
ATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCGGTTGGGATCAGAATACGTTCGAGGAAGGTGTAAACCCAGACCCCAAGTGCACCCGCCGAACGCAGGAACGCTTGCAGAGATTCAATCCCCATTTGTACTTTTGGCCAGCCGAGCAACGTCAGCCAGGCACAGGGGATCATCACCAGGAAGGCGATAATCACCACATAAGACGTCCCCTGGAAAATGCCGAGAAAAACCGGCAGTTTTTTATCGAACAGACGGTTATGCAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGAT  >  NC_000913/3862732‑3863183

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