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→GGGglnLsensory 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

CTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGT  >  NC_000913/4055880‑4056346
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
cTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATCCGAGCGCCCGTCGATGACCAGCGTCCCTTCGTTATCGGTAAAACCt                                                                                                                                                                                                                                         >  2:1017083/1‑238 (MQ=255)
     ccgcgccACGTAAACGGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTACCGGATCCTGACTTAAGCGGCGCTGGGTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCAGCAATACGCTGGGCCGTCCCAGAAAGGTTATACAGGCGCCCGCGGCAGACCCGGGTCCCTTCGATACCGGAAAAACCtacccgc                                                                                                                                                                                                                                   >  1:903190/1‑233 (MQ=255)
                                tttttctctcCTGGGCCGCGCCGCGCACTTAACCGCGGGGAGCCACCTGTTGGGGGGGGGGTTGGTGTTCCCGAACTGAGGCGGCGTTGTTTTCCCACCGGAGCCTCCCCAGAAGGGTCCTTGCCGCCCGGATTCCGCTGGGCCGTCACAGAAAGGATATGCGCGCCCCCGTCGAGGACCAGCGTCACTTCGTTGTCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTc                                                                                                                                                                                                        <  1:1000804/233‑1 (MQ=255)
                                ttttAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGa                                                                                                                                                                                                                         >  2:918680/1‑221 (MQ=255)
                                          aTGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGAGATTTAAt                                                                                                                                                                                              >  2:1000804/1‑238 (MQ=255)
                                                    ccGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGAAGCt                                                                                                                                                                                   >  2:199528/1‑239 (MQ=255)
                                                    ccGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCCCGTATTTTAATGGGAAGTAGcc                                                                                                                                                                                   >  2:534203/1‑238 (MQ=255)
                                                       cgcACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTaaa                                                                                                                                                                                >  2:979160/1‑239 (MQ=255)
                                                       cgcACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGAAGCTCac                                                                                                                                                                                >  2:868143/1‑238 (MQ=255)
                                                              aaTCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATTCGAGCGCCAGTCGATGACAAGCGTCACTTCGTTCTCGGTAAACCCTTGCACCGCCTCCAGGCTTTCCTGCTTCAGCTCGTATTTTAATGGAAAGTAGccaaccaggtac                                                                                                                                                                         >  1:934926/1‑228 (MQ=255)
                                                                    gggCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCCCAAAAGTTCCGGTAAc                                                                                                                                                                   >  2:928321/1‑239 (MQ=255)
                                                                      gcagcaACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATCCCGTCCGGCATACGCTGGGCCGGCACAGAAAGGATATGCGAGCGCCCGTCGATGACCACCGTCACTTCGTTCTCGGTCAAACCTGGCCCCGCCTCAAGCCTTTCTTGCTTCAGCTCGATATTTAATCAGAAGGAGCCCAACAGtgcagggaaagg                                                                                                                                                                 >  1:751876/1‑228 (MQ=255)
                                                                        agcaACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAACGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTCCTTGCCTCACCTCCATCTTTAATTAACAGTATCCAAACCTTCCCGGTAActgg                                                                                                                                                                >  1:316124/1‑235 (MQ=255)
                                                                         gcaACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATAGTCGAGAGCCCGTAGATGACCATCGTAACTACGTTCTCGGGAAACCCTTGCCCCACCTCCCGACTTTCTTGCATCTCCACCATAATTAATAAAAAGTGACCCAACAGTTCCGGAAACAgtgt                                                                                                                                                              >  1:762661/1‑239 (MQ=255)
                                                                               tgtgggcgcgtttGTATCTGTCCTTGATTTACGCGGCGCTTGTTACCCAGCGGGCCCCTCTCCAGGAGTACTCGGCCGCCCGCCATACGTCGGCCCTGCACAGAAAGGTTTGGCGCGCGCCCGCGGAGAACCAGCGTCCTTTCTTTATGGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCaaa                                                                                                                                                        <  1:1011286/227‑1 (MQ=255)
                                                                                        tGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTTCg                                                                                                                                               >  2:256158/1‑239 (MQ=255)
                                                                                                   ttCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAACACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCACCAGTTCCGGTAACGGTGTACCCAACCATTTTGGGGAGCTTTtgg                                                                                                                                    >  1:931295/1‑239 (MQ=255)
                                                                                                       tGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGCTGACCAGCGTCACTTCGTAATCGGTAAAACCTCGCCCCGCCTCCAGACTTTCTTGCATAAGCTCGATATTTAATGAGAAGTAGCGCAACAGTTCCGGTAAAGGTGAACCAAACAATTTACGGAGGCTTATGGCGa                                                                                                                                 >  1:211235/1‑238 (MQ=255)
                                                                                                              aaGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCCCAACAGTTCCGGCAACGGGTTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttg                                                                                                                          >  1:54711/1‑239 (MQ=255)
                                                                                                              aaGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGAAATTTATTGGGAAGTAGCCCAACATTCCCGGTAACGGGTAACCAACCAATTTGCGGGGGCTTTGGGCGGCCAgttt                                                                                                                          >  2:611810/1‑238 (MQ=255)
                                                                                                               aGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACGGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGGGGGAGCTTTGGGCGGGCGgttgt                                                                                                                         >  1:407989/1‑239 (MQ=255)
                                                                                                               aGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCCCTTCGTTATCGGTAAAACCTTGCCCCGCCTCCCGCCTTTCTTGACTCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACCGTGTACCAAACAATTTGCGGGAGCTTTGGGGcagcaggttg                                                                                                                         >  1:199076/1‑235 (MQ=255)
                                                                                                                  ggggttGGTTATCCAGCGGGCCCATCCCCAGCAGGATCATGCCGCCCGGCATACGCTGGGCCGTCAAAGAAAGGATTGGCGAGCCCCCGTCGATGCCCAGCGTCACTTCGTTACCGGAAAACCCTGGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttgtt                                                                                                                        <  2:962346/232‑1 (MQ=255)
                                                                                                                    cgcTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCCCTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGGGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAACCAATTTCCGGAAGCTTTGGGCGCGCAgttgttcgca                                                                                                                    >  2:707124/1‑235 (MQ=255)
                                                                                                                    cgcTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGCAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGATTTGGGGGAGCAGttgttgggg                                                                                                                    >  1:730879/1‑236 (MQ=255)
                                                                                                                      cTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGcgcg                                                                                                                   >  2:729690/1‑238 (MQ=255)
                                                                                                                            tccctgcGGGGCCATCTCCAGTCGGACCCTGCCGTCCGGAATTAGTGGGGCCGTCCAAGAAAGGATACGCGACCGCCCGTCGTCGCCAAGCGTCACTTCGTTATCGGTAACCCCTTGCCCCGCCTCCAGACTGTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGt                                                                                                            <  1:912812/233‑1 (MQ=255)
                                                                                                                            acccagggggcccATCCCCAGAGGGAACAGCCCGCCGGGCACGCGGTGGCCCGTCACAAAAGGGATAGGCGAGCGCCCGTCGATGCCCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGt                                                                                                            <  2:926469/228‑1 (MQ=255)
                                                                                                                                  cGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGCAACGGTGTACCAAACAATTTGCGGGGACTTTTGGCGGGAAGGTGTTGGGCGGCAgggggggggt                                                                                                      >  1:156882/1‑232 (MQ=255)
                                                                                                                                   ggAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCAAGTAGTTGTTgccaggaccggctggcgtc                                                                                                     >  2:218305/1‑222 (MQ=255)
                                                                                                                                       ccATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCGGTTGTTGCGCGGCAGGGTTGGCGTATTgg                                                                                                 >  2:72864/1‑239 (MQ=255)
                                                                                                                                           ctcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAACCAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGAGCGGCAGGGTTGGCGTAATGGATc                                                                                              >  2:360301/1‑238 (MQ=255)
                                                                                                                                                         aTGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTACCGGTGTACCAAACAATTTCCGGGAGCTTTGGGCGAGCCGTTTTTTCCCGCCCGGGTTGGCGTAATGGATCGCCAGGTTGtatc                                                                                >  1:374283/1‑235 (MQ=255)
                                                                                                                                                          tGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTGGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGGGAAGTAGCTCAACAGTTCCGGTAACGGTGTCCCAAACAATTTGCGGGAGCTTTGGGCGCGCAGTTGTTGCGCGGCGGGGTTGGGGAAATggtacgtccgggtgtgatgg                                                                               >  1:641903/1‑221 (MQ=255)
                                                                                                                                                            ccGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCCGGGTTGGCGTAATGGATCCCCAGGTTGTCATCgct                                                                            >  2:342809/1‑237 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGa                                                                                                                                                                                                                         >  2:898189/1‑96 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTTGGCGAGAAGTTGTGGCGCGGAGGGGTTGGGGTAATGGATCGCCGGGTTGTCATCGATTaa                                                                         >  1:873602/1‑239 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGTAGAGCAGTTGTTGCGCGGGAGGGTTGGCGTAATGGaagggcggggtggcatgatataa                                                                          >  1:878997/1‑217 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGAACGCCAGGGTGGCATGGaaaaa                                                                         >  1:404168/1‑235 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGCAACGGGGTACCAAACAATTTGCGGGAGCCTTGGGCGCGCAGTTGTTGAGCGGCCAGGTTGGCGTAATGGTTAGCCAGGGTCTAATCGATTaa                                                                         >  1:191457/1‑239 (MQ=255)
                                                                                                                                                               tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCACGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTGGTTACACGGCGGGGGTGGCGTAAAGGAACGCCAGGGTTGTACCGATTaa                                                                         >  1:920219/1‑239 (MQ=255)
                                                                                                                                                                     ctatcccgggcccgcccaGAAAAGATAAGCGGGCGCCCGGCGATAGCCCTGGGTTCCTCCTTATTGGTAAAACCCTGCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGTGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACaaaa                                                                    <  1:868143/223‑1 (MQ=255)
                                                                                                                                                                     attcgcgGGGCCTGCACAAAAATGATATGCGGGCGGCCGGCGATGCCCAGCGTTACTTCGTTACCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACaaaa                                                                    <  2:211235/235‑1 (MQ=255)
                                                                                                                                                                     atccggtGGCCCGACACAAAAGGCATATGCGCGCGCCTGTCGATGACCAGCGTCACTGCGTTACCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACaaaa                                                                    <  2:954066/232‑1 (MQ=255)
                                                                                                                                                                              gCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAAAACGGTTaa                                                           >  2:805975/1‑238 (MQ=255)
                                                                                                                                                                                 ttCACAAAAAGGATATGCGGGCCCCCGTCGATGACCAGCGTCACTTCTTTACCGGAAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAg                                                       <  2:404168/238‑1 (MQ=255)
                                                                                                                                                                                 tccaaaaaaaGGATTTGCGTGCCTCCGTCGTTCACCCGCCCTCCGTTGTTCTCGATAAACCCTGCCCCCTCCCCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAGCAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAg                                                       <  1:12528/232‑1 (MQ=255)
                                                                                                                                                                                  caacaGAAGGGATTTGCGAGCGCCCGTCAGTCCCCAGCGTCACTTCGTTTCGGGTAAACCTTTCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGc                                                      <  2:903190/237‑1 (MQ=255)
                                                                                                                                                                                   acaaaaaatgttctGCGCGCCCCCGTCGATGCCCACCTCCATTTTGTTCTCGGTAAACCCTTCCCCCGCCTCCAGACTTTCTTTCATCAGCTCGATTTTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCg                                                     <  1:997891/226‑1 (MQ=255)
                                                                                                                                                                                               atGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATAAGCGAGTTGAGGATc                                          >  2:823411/1‑238 (MQ=255)
                                                                                                                                                                                                   gAGCGCCCGTGGAGGCCCGGCTTCATTTCGTTAGCGGAAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGcc                                      <  2:199076/238‑1 (MQ=255)
                                                                                                                                                                                                   gAGCGCCCGGCGGAGCCCAGCTTCACTTGGTTACCGGCAAACCTTTGCGCCGCCTCCAGGCTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGcc                                      <  2:751876/238‑1 (MQ=255)
                                                                                                                                                                                                      cgcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTACTGGATAGCCAGGTTGTCATCGATTAACAAAATACTGTTAATAAGCGAGTTGAGGATATTCCCCGc                                  >  2:810765/1‑239 (MQ=255)
                                                                                                                                                                                                      cgcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCggc                                                                                           >  2:764337/1‑181 (MQ=255)
                                                                                                                                                                                                      cgcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGc                                  >  2:695440/1‑239 (MQ=255)
                                                                                                                                                                                                                    ccAGCGTCACTTCTTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGc                                                      <  2:731773/206‑1 (MQ=255)
                                                                                                                                                                                                                    ccAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGc                                                      >  1:731773/1‑206 (MQ=255)
                                                                                                                                                                                                                    ccAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATAAGCGAGTTGAGGATTTGCCCAGCATCGGGCTGCGt                      >  2:944907/1‑238 (MQ=255)
                                                                                                                                                                                                                       gCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGGTTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCACCATCGGGCTGCGTg                     >  1:205061/1‑236 (MQ=255)
                                                                                                                                                                                                                        cGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAATTACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGGGccgg                 >  2:856397/1‑237 (MQ=255)
                                                                                                                                                                                                                        cGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCGGCATCGGGCTggtgg                     >  2:910416/1‑231 (MQ=255)
                                                                                                                                                                                                                        cGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTg                     >  2:888148/1‑235 (MQ=255)
                                                                                                                                                                                                                         tcctcTTCGTTTCCGATAAACCCTGGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGt                <  2:513088/235‑1 (MQ=255)
                                                                                                                                                                                                                         gtcattctgttacgggTAAACCCTTACCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTg                 <  2:934926/224‑1 (MQ=255)
                                                                                                                                                                                                                          tttcttttttttgcggtaaaCCTTGCCCCCCCCTCCAGCCTTTCTTGCATCAGCTCGATATTTAATGAGAAGGAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGt                <  1:695440/222‑1 (MQ=255)
                                                                                                                                                                                                                             cTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAgg                                                                                       >  2:979483/1‑164 (MQ=255)
                                                                                                                                                                                                                                cGTTAGCGGTAACACCGTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATa         <  2:878997/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tgttaaaacCCTTCCCCCGCCTCCAATCTTCTTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGGTGCGGTAACGGTGTACCAAACTATTTGCGGGATCTTTGGGCGAGCAGGTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAATCa    <  1:72864/234‑1 (MQ=255)
                                                                                                                                                                                                                                     tCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGGTCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCa    >  2:605851/1‑239 (MQ=255)
                                                                                                                                                                                                                                      cGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCTTAAAGCAg   >  2:129075/1‑239 (MQ=255)
                                                                                                                                                                                                                                       ggTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGt  <  2:407989/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
CTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGT  >  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: