Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,056,114 A→C 22.8% G75G (GGT→GGG)  glnL ← sensory histidine kinase in two‑component regulatory system with GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,056,1140AC22.8% 145.1 / 8.6 69G75G (GGT→GGG) glnLsensory histidine kinase in two‑component regulatory system with GlnG
Reads supporting (aligned to +/- strand):  ref base A (48/3);  new base C (2/16);  total (50/19)
Fisher's exact test for biased strand distribution p-value = 6.94e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGT  >  NC_000913/4055880‑4056346
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
cTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATCCGAGCGCCCGTCGAT‑GACCAGCGTCCCTTCGTTATCGGTAAAACCt                                                                                                                                                                                                                                         >  2:1017083/1‑238 (MQ=255)
     ccgcgccACGTAAACGGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTACCGGATCCTG‑ACTTAAGCGGCGCTGGGTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCAGCAATACGCTGGGCCGTCCCAGAAAGGTTATACAGGCGCCCGCGGCA‑GACCCGGGTCCCTTCGATACCGGAAAAACCtacccgc                                                                                                                                                                                                                                   >  1:903190/1‑233 (MQ=255)
                                tttttctctcCTGGGCCGCGCCGCGCACTTAACCGCGGGGAGCCACCTGTTGGGGGGGGGGTTGGTGTTCCCGAACTGAGGCGGCGTTGTTTTCCCACCGGAGCC‑TCCCCAGAAGGGTCCTTGCCGCCCGGATTCCGCTGGGCCGTCACAGAAAGGATATGCGCGCCCCCGTCGAG‑GACCAGCGTCACTTCGTTGTCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTc                                                                                                                                                                                                        <  1:1000804/233‑1 (MQ=255)
                                ttttAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGa                                                                                                                                                                                                                         >  2:918680/1‑221 (MQ=255)
                                          aTGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGAGATTTAAt                                                                                                                                                                                              >  2:1000804/1‑238 (MQ=255)
                                                    ccGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGAAGCt                                                                                                                                                                                   >  2:199528/1‑239 (MQ=255)
                                                    ccGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCCCGTATTTTAATGGGAAGTAGcc                                                                                                                                                                                   >  2:534203/1‑238 (MQ=255)
                                                       cgcACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTaaa                                                                                                                                                                                >  2:979160/1‑239 (MQ=255)
                                                       cgcACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGAAGCTCac                                                                                                                                                                                >  2:868143/1‑238 (MQ=255)
                                                              aaTCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATTCGAGCGCCAGTCGAT‑GACAAGCGTCACTTCGTTCTCGGTAAACCCTTG‑CACCGCCTCCAGGCTTTCCTGCTTCAGCTCGTATTTTAATGGAAAGTAGccaaccaggtac                                                                                                                                                                         >  1:934926/1‑228 (MQ=255)
                                                                    gggCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCCCAAAAGTTCCGGTAAc                                                                                                                                                                   >  2:928321/1‑239 (MQ=255)
                                                                      gcagcaACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATCCCGTCCGGCATACGCTGGGCCGGCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCACCGTCACTTCGTTCTCGGTCAAACCTGG‑CCCCGCCTCAAGCCTTTCTTGCTTCAGCTCGATATTTAATCAGAAGGAGCCCAACAGtgcagggaaagg                                                                                                                                                                 >  1:751876/1‑228 (MQ=255)
                                                                        agcaACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAACGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTCCTTGCCTCACCTCCATCTTTAATTAACAGTATCCAAACCTTCCCGGTAActgg                                                                                                                                                                >  1:316124/1‑235 (MQ=255)
                                                                         gcaACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATAGTCGAGAGCCCGTAGAT‑GACCATCGTAACTACGTTCTCGGGAAACCCTTG‑CCCCACCTCCCGACTTTCTTGCATCTCCACCATAATTAATAAAAAGTGACCCAACAGTTCCGGAAACAgtgt                                                                                                                                                              >  1:762661/1‑239 (MQ=255)
                                                                               tgtgggcgcgtttGTATCTGTCCTTG‑ATTTACGCGGCGCTTGTTACCCAGCGGGCCCCTCTCCAGGAGTAC‑TCGGCCGCCCGCCATACGTCGGCCCTGCACAGAAAGGTTTGGCGCGCGCCCGCGGAG‑AACCAGCGTCCTTTCTTTATGGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCaaa                                                                                                                                                        <  1:1011286/227‑1 (MQ=255)
                                                                                        tGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTTCg                                                                                                                                               >  2:256158/1‑239 (MQ=255)
                                                                                                   ttCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAACACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCACCAGTTCCGGTAACGGTGTACCCAACCATTTTGGGGAGCTTTtgg                                                                                                                                    >  1:931295/1‑239 (MQ=255)
                                                                                                       tG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGCT‑GACCAGCGTCACTTCGTAATCGGTAAAACCTCG‑CCCCGCCTCCAGACTTTCTTGCATAAGCTCGATATTTAATGAGAAGTAGCGCAACAGTTCCGGTAAAGGTGAACCAAACAATTTACGGAGGCTTATGGCGa                                                                                                                                 >  1:211235/1‑238 (MQ=255)
                                                                                                              aaGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCCCAACAGTTCCGGCAACGGGTTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttg                                                                                                                          >  1:54711/1‑239 (MQ=255)
                                                                                                              aaGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGAAATTTATTGGGAAGTAGCCCAACATTCCCGGTAACGGGTAACCAACCAATTTGCGGGGGCTTTGGGCGGCCAgttt                                                                                                                          >  2:611810/1‑238 (MQ=255)
                                                                                                               aGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACGGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGGGGGAGCTTTGGGCGGGCGgttgt                                                                                                                         >  1:407989/1‑239 (MQ=255)
                                                                                                               aGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCCCTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCCGCCTTTCTTGACTCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACCGTGTACCAAACAATTTGCGGGAGCTTTGGGGcagcaggttg                                                                                                                         >  1:199076/1‑235 (MQ=255)
                                                                                                                  ggggttGGTTATCCAGCGGGCCCATCCCCAGCAGGAT‑CATGCCGCCCGGCATACGCTGGGCCGTCAAAGAAAGGATTGGCGAGCCCCCGTCGAT‑GCCCAGCGTCACTTCGTTACCGGAAAACCCTGG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttgtt                                                                                                                        <  2:962346/232‑1 (MQ=255)
                                                                                                                    cgcTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCCCTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGGGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAACCAATTTCCGGAAGCTTTGGGCGCGCAgttgttcgca                                                                                                                    >  2:707124/1‑235 (MQ=255)
                                                                                                                    cgcTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGCAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGATTTGGGGGAGCAGttgttgggg                                                                                                                    >  1:730879/1‑236 (MQ=255)
                                                                                                                      cTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGcgcg                                                                                                                   >  2:729690/1‑238 (MQ=255)
                                                                                                                            tccctgcGGGGCCATCTCCAGTCGGAC‑CCTGCCGTCCGGAATTAGTGGGGCCGTCCAAGAAAGGATACGCGACCGCCCGTCGTC‑GCCAAGCGTCACTTCGTTATCGGTAACCCCTTG‑CCCCGCCTCCAGACTGTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGt                                                                                                            <  1:912812/233‑1 (MQ=255)
                                                                                                                            acccagggggcccATCCCCAGAGGGAA‑CAGCCCGCCGGGCACGCGGTGGCCCGTCACAAAAGGGATAGGCGAGCGCCCGTCGAT‑GCCCAGCGTCACTTCGTTATCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGt                                                                                                            <  2:926469/228‑1 (MQ=255)
                                                                                                                                  cGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGCAACGGTGTACCAAACAATTTGCGGGGACTTTTGGCGGGAAGGTGTTGGGCGGCAgggggggggt                                                                                                      >  1:156882/1‑232 (MQ=255)
                                                                                                                                   ggAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCAAGTAGTTGTTgccaggaccggctggcgtc                                                                                                     >  2:218305/1‑222 (MQ=255)
                                                                                                                                       ccATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GCCCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCGGTTGTTGCGCGGCAGGGTTGGCGTATTgg                                                                                                 >  2:72864/1‑239 (MQ=255)
                                                                                                                                           ctcCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAACCAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGAGCGGCAGGGTTGGCGTAATGGATc                                                                                              >  2:360301/1‑238 (MQ=255)
                                                                                                                                                         aTGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTACCGGTGTACCAAACAATTTCCGGGAGCTTTGGGCGAGCCGTTTTTTCCCGCCCGGGTTGGCGTAATGGATCGCCA‑GGTTGtatc                                                                                >  1:374283/1‑235 (MQ=255)
                                                                                                                                                          tGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTGGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGGGAAGTAGCTCAACAGTTCCGGTAACGGTGTCCCAAACAATTTGCGGGAGCTTTGGGCGCGCAGTTGTTGCGCGGCGGGGTTGGGGAAATggtacgtccgggtgtgatgg                                                                               >  1:641903/1‑221 (MQ=255)
                                                                                                                                                            ccGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCCGGGTTGGCGTAATGGATCCCCA‑GGTTGTCATCgct                                                                            >  2:342809/1‑237 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGa                                                                                                                                                                                                                         >  2:898189/1‑96 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTTGGCGAGAAGTTGTGGCGCGGAGGGGTTGGGGTAATGGATCGCCG‑GGTTGTCATCGATTaa                                                                         >  1:873602/1‑239 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGTAGAGCAGTTGTTGCGCGGGAGGGTTGGCGTAATGGaagggcggggtggcatgatataa                                                                          >  1:878997/1‑217 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGAACGCCA‑GGGTGGCATGGaaaaa                                                                         >  1:404168/1‑235 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGCAACGGGGTACCAAACAATTTGCGGGAGCCTTGGGCGCGCAGTTGTTGAGCGGCCAGGTTGGCGTAATGGTTAGCCA‑GGGTCTAATCGATTaa                                                                         >  1:191457/1‑239 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCACGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTGGTTACACGGCGGGGGTGGCGTAAAGGAACGCCAGGGTTGT‑ACCGATTaa                                                                         >  1:920219/1‑239 (MQ=255)
                                                                                                                                                                     ctatcccgggcccgcccaGAAAAGATAAGCGGGCGCCCGGCGATAGCCCTGGGT‑TCCTCCTTATTGGTAAAACCCTGCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGTGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACaaaa                                                                    <  1:868143/223‑1 (MQ=255)
                                                                                                                                                                     attcgcgGGGCCTGCACAAAAATGATATGCGGGCGGCCGGCGAT‑GCCCAGCGTTACTTCGTTACCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACaaaa                                                                    <  2:211235/235‑1 (MQ=255)
                                                                                                                                                                     atccggtGGCCCGACACAAAAGGCATATGCGCGCGCCTGTCGAT‑GACCAGCGTCACTGCGTTACCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACaaaa                                                                    <  2:954066/232‑1 (MQ=255)
                                                                                                                                                                              gCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAAAACGGTTaa                                                           >  2:805975/1‑238 (MQ=255)
                                                                                                                                                                                 ttCACAAAAAGGATATGCGGGCCCCCGTCGAT‑GACCAGCGTCACTTCTTTACCGGAAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAg                                                       <  2:404168/238‑1 (MQ=255)
                                                                                                                                                                                 tccaaaaaaaGGATTTGCGTGCCTCCGTCGTT‑CACCCGCCCTCCGTTGTTCTCGATAAACCCTGC‑CCCCTCCCCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAGCAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAg                                                       <  1:12528/232‑1 (MQ=255)
                                                                                                                                                                                  caacaGAAGGGATTTGCGAGCGCCCGTCAGT‑CCCCAGCGTCACTTCGTTTCGGGTAAACCTTTC‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGc                                                      <  2:903190/237‑1 (MQ=255)
                                                                                                                                                                                   acaaaaaatgttctGCGCGCCCCCGTCGAT‑GCCCACCTCCATTTTGTTCTCGGTAAACCCTTC‑CCCCGCCTCCAGACTTTCTTTCATCAGCTCGATTTTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCg                                                     <  1:997891/226‑1 (MQ=255)
                                                                                                                                                                                               atGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATAAGCGAGTTGAGGATc                                          >  2:823411/1‑238 (MQ=255)
                                                                                                                                                                                                   gAGCGCCCGTGGAG‑GCCCGGCTTCATTTCGTTAGCGGAAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGcc                                      <  2:199076/238‑1 (MQ=255)
                                                                                                                                                                                                   gAGCGCCCGGCGGA‑GCCCAGCTTCACTTGGTTACCGGCAAACCTTTG‑CGCCGCCTCCAGGCTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGcc                                      <  2:751876/238‑1 (MQ=255)
                                                                                                                                                                                                      cgcCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTACTGGATAGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATAAGCGAGTTGAGGATATTCCCCGc                                  >  2:810765/1‑239 (MQ=255)
                                                                                                                                                                                                      cgcCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCggc                                                                                           >  2:764337/1‑181 (MQ=255)
                                                                                                                                                                                                      cgcCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGc                                  >  2:695440/1‑239 (MQ=255)
                                                                                                                                                                                                                    ccAGCGTCACTTCTTTATCGGTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGc                                                      <  2:731773/206‑1 (MQ=255)
                                                                                                                                                                                                                    ccAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGc                                                      >  1:731773/1‑206 (MQ=255)
                                                                                                                                                                                                                    ccAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATAAGCGAGTTGAGGATTTGCCCAGCATCGGGCTGCGt                      >  2:944907/1‑238 (MQ=255)
                                                                                                                                                                                                                       gCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGGTTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCACCATCGGGCTGCGTg                     >  1:205061/1‑236 (MQ=255)
                                                                                                                                                                                                                        cGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAATTACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGGGccgg                 >  2:856397/1‑237 (MQ=255)
                                                                                                                                                                                                                        cGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCGGCATCGGGCTggtgg                     >  2:910416/1‑231 (MQ=255)
                                                                                                                                                                                                                        cGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTg                     >  2:888148/1‑235 (MQ=255)
                                                                                                                                                                                                                         tcctcTTCGTTTCCGATAAACCCTGG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGt                <  2:513088/235‑1 (MQ=255)
                                                                                                                                                                                                                         gtcattctgttacgggTAAACCCTTA‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTg                 <  2:934926/224‑1 (MQ=255)
                                                                                                                                                                                                                          tttcttttttttgcggtaaaCCTTGCCCCCCCCTCCAGCCTTTCTTGCATCAGCTCGATATTTAATGAGAAGGAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGt                <  1:695440/222‑1 (MQ=255)
                                                                                                                                                                                                                             cTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑gg                                                                                       >  2:979483/1‑164 (MQ=255)
                                                                                                                                                                                                                                cGTTAGCGGTAACACCGTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATa         <  2:878997/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tgttaaaacCCTTC‑CCCCGCCTCCAATCTTCTTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGGTGCGGTAACGGTGTACCAAACTATTTGCGGGATCTTTGGGCGAGCAGGTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAATCa    <  1:72864/234‑1 (MQ=255)
                                                                                                                                                                                                                                     tCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGGTCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCa    >  2:605851/1‑239 (MQ=255)
                                                                                                                                                                                                                                      cGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCTTAAAGCAg   >  2:129075/1‑239 (MQ=255)
                                                                                                                                                                                                                                       ggTAAACCCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGt  <  2:407989/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
CTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTG‑ACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGAT‑CATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAT‑GACCAGCGTCACTTCGTTATCGGTAAAACCTTG‑CCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCA‑GGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGT  >  NC_000913/4055880‑4056346

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: ‑