Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,749 A→G A18A (GCA→GCG creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,636,7490AG98.3% 181.0 / NA 60A18A (GCA→GCGcreCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (23/36);  total (24/36)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCG  >  NC_000913/4636504‑4636995
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
aCTCAAGTCGCCAGGCCGCGTCTGGTCCCTCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTg                                                                                                                                                                                                                                                   <  1:273640/251‑1 (MQ=255)
   caaggcgccaggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGtt                                                                                                                                                                                                                                                <  1:143138/241‑1 (MQ=255)
    aaGTCGCCAGGCTGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCTGGTCACCTACGATCGCACCGTTGATTCCCACATTAAAACGTTGCGCGCCAAGCTGCGCGCCATCAACCCCGATCTTGCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTc                                                                                                                                                                                                                                               <  2:578310/251‑1 (MQ=255)
         gccaggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACt                                                                                                                                                                                                                                          <  1:24130/247‑1 (MQ=255)
          ccaggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTg                                                                                                                                                                                                                                         >  1:455781/4‑251 (MQ=255)
          ccaggCCGCGTCTGGTACCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTg                                                                                                                                                                                                                                         >  1:223290/4‑251 (MQ=255)
                   gtcgggTCCCGCCAGCAACTAATGGATAGCGGGTGGGAAGATGCGCAGGACACCGACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTAttttt                                                                                                                                                                                                                                <  2:92914/247‑1 (MQ=255)
                          cccGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATTCTTATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCaaa                                                                                                                                                                                                                          <  2:220967/250‑1 (MQ=255)
                                   aaCTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCTGCATGGGATATCGCCTGAGGGGCCTGTCATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGTCGGCCTGGTTCGTACCGGCTATTTTTGTCTAAGAAGTTAAAc                                                                                                                                                                                                                >  1:305374/1‑251 (MQ=255)
                                       aaTGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCggg                                                                                                                                                                                                            <  2:204090/250‑1 (MQ=255)
                                               gCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATCTAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCgc                                                                                                                                                                                                     >  2:429324/1‑249 (MQ=255)
                                                cgtgtGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAAACGGGCGTGCGa                                                                                                                                                                                                     >  1:581836/5‑249 (MQ=255)
                                                               gcAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAAcc                                                                                                                                                                                                               <  1:452607/224‑1 (MQ=255)
                                                               gcAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTGACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGAc                                                                                                                                                                                    <  1:536207/251‑1 (MQ=255)
                                                                       ccTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCg                                                                                                                                                                            >  1:588898/1‑251 (MQ=255)
                                                                             aTCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCGACGGAGGGGACGTTAATCGATACCg                                                                                                                                                                      >  1:283836/1‑251 (MQ=255)
                                                                             aTCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCg                                                                                                                                                                      <  1:478890/251‑1 (MQ=255)
                                                                                 cACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATGTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAAc                                                                                                                                                                  <  2:312264/251‑1 (MQ=255)
                                                                                  cccGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCGATTTTTGTCAAAGAAGTTAAACCGTGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACg                                                                                                                                                                 <  2:290044/250‑1 (MQ=255)
                                                                                     gttgATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGCTATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCTGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTg                                                                                                                                                              >  1:279495/4‑251 (MQ=255)
                                                                                        gATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCt                                                                                                                                                            <  1:350513/250‑1 (MQ=255)
                                                                                         aTACCCACCTTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTgg                                                                                                                                                          >  1:518781/1‑251 (MQ=255)
                                                                                         aTACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCTGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTgg                                                                                                                                                          <  2:455781/251‑1 (MQ=255)
                                                                                          tACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGc                                                                                                                                                         >  1:158250/1‑251 (MQ=255)
                                                                                             ccACACTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGGAGAGCAACGGGGGGGACGTTAATCGCTACCGCAACGTTGATGGCggg                                                                                                                                                      >  2:417533/1‑250 (MQ=255)
                                                                                               acaTTAAAACGTTGCGTGCCAGGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCAGGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTAGTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACTTTGCTGGCGGAGc                                                                                                                                                    <  2:223290/251‑1 (MQ=255)
                                                                                                caTTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCt                                                                                                                                                   <  1:166902/251‑1 (MQ=255)
                                                                                                    aaaaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTgg                                                                                                                                                 <  1:282133/249‑1 (MQ=255)
                                                                                                    aaaaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGgc                                                                                                                                                <  1:131973/250‑1 (MQ=255)
                                                                                                    aaaaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAAACGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATTGATACCGCAACGTTGCTGGCGGAGCTgg                                                                                                                                                 >  2:282133/1‑249 (MQ=255)
                                                                                                       acgttGCGTGCCAAGCTGCGCGCCATCAACCCCGCTCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGGGCAACGGGGGGGACGTTAATCGATCCCGCAACGTTGCTGGCGGAGCTGGcgcgg                                                                                                                                            >  2:573140/5‑250 (MQ=255)
                                                                                                        cgttGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTc                                                                                                                                           <  2:158250/248‑1 (MQ=255)
                                                                                                               ggCCAAGCTGCGCGCCATCAACCCCGAGCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCGGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGAttt                                                                                                                                    <  2:328963/250‑1 (MQ=255)
                                                                                                                gCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATAGATAACGCAACGTTGCTGGCGGAGCTGGCGCGTCCCCATTTg                                                                                                                                   >  1:445098/1‑251 (MQ=255)
                                                                                                                                tCAACCCCGCTCTTACACCGATTAATACTCCTCCCGGCGTGGGATATAGCCTGAGGGGCCTGTAATGCGAATCGGCCTGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGATTGTACTGGCTATTTTTGTCAAAGCAGTGCAAACGGGCGAGCGAAGAGAAACGGATGGGAAGTTACTCGATACCGCAACGTTTCTGGGGGAGCTGGCGCGGCCCGATTTGttctcgaggggcaca                                                                                                                    >  2:430539/1‑240 (MQ=255)
                                                                                                                                   aCCCCGCTCTTTCACCGATTAATACTCATCGCGGCCTGCGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCAATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGACGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGGGCTGGCGCGTCCCGATTTGCTCTCTGGGGGCCCAACGc                                                                                                                >  2:270469/1‑251 (MQ=255)
                                                                                                                                     cccGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCAt                                                                                                              <  1:149188/251‑1 (MQ=255)
                                                                                                                                          tCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGc                                                                                                          >  1:12961/1‑250 (MQ=255)
                                                                                                                                               cACCGATTAATCCTCATCGCTGCATGGTATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCATCTg                                                                                                     <  2:431702/250‑1 (MQ=255)
                                                                                                                                                    aTTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTGTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAg                                                                                               <  1:337207/251‑1 (MQ=255)
                                                                                                                                                                      aTGGGATATAGCCTGAGGGGACTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGACTGGTTCGAACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTTGGAAGAGCACCGGAGGGGGCGTTAAACGATACCGCAACGTTGCTGGCGGGGCAGGGGCGTCACGATTTGCTCTCTGGGGGCCCAAAGCAAGGGCAACGGGCGGCGGCGTTTAATCGGCTaca                                                                              >  2:461050/1‑250 (MQ=255)
                                                                                                                                                                          gATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATGTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATc                                                                         <  1:627592/251‑1 (MQ=255)
                                                                                                                                                                             ataGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGcc                                                                      <  2:130873/251‑1 (MQ=255)
                                                                                                                                                                                    gaagggCCTGTAATGCGTATCGGCAAGCGGTTGCTGCTGGGGTATTTTTTACTGGTGGCGTTGGCGGCCTGGTGCGTGCTGGCTAGTTTTGTCAAAGAGGGTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACATCGTCGCCCGTTTcg                                                               <  2:529579/248‑1 (MQ=255)
                                                                                                                                                                                    gAGGGGCCTGTAATGCTTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcg                                                               <  2:12238/251‑1 (MQ=255)
                                                                                                                                                                                    gAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGttt                                                                 <  1:346002/249‑1 (MQ=255)
                                                                                                                                                                                         gCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGcc                                                           >  1:43555/1‑250 (MQ=255)
                                                                                                                                                                                               aaTGCGTATCGGCATGCGTGTGCTGCTGGGCTATTTTTTACGGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCGCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCATCTACAACATCGCCATTTTCGCGCCAATATc                                                     <  2:445098/250‑1 (MQ=255)
                                                                                                                                                                                                          gCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAAc                                         >  1:9941/1‑251 (MQ=255)
                                                                                                                                                                                                           cATGCGGTTGTTGCTGGGCTATTTTTTTCTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATGTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATTCCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCGGGGGACCCAACGCTTGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCTATTCGCGCCAATATCTGTGTCATTAACa                                        <  2:329279/251‑1 (MQ=255)
                                                                                                                                                                                                                     ctgCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGGTCGTACTGGCTATTTTTGTCAAAGGAGTTAAACCGGGCGGGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGGGCGTCCCGATTTGCTCTCCGGGGACCCAACGCATGGGCAAATTGGGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTAGGGCCAATATCGGTGGCATTAACAAAGAGCGCa                               >  2:515361/2‑250 (MQ=255)
                                                                                                                                                                                                                          gggCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGATGAGCAACGGACGGGACCTTAATCGGTACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGAACCACCGCCTGGGCAACTGGGGCAGGCGTTTAATCAGCCACAACATCGCCCGTTTCGCCCCACAATCGGAGGCATTCACAAAGTGCGCAATggct                         >  1:291462/1‑248 (MQ=255)
                                                                                                                                                                                                                              gaTTTTTTACTGGTGGCGGTGGCGGCCGGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGGTTAATCAGCTACAACATCGCCCGGTTCGCGCCAATATCGGTGGCGTTAACAAAGTGCGCAATGAATATCa                     <  2:89430/250‑1 (MQ=255)
                                                                                                                                                                                                                               aTTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAATGCGCAATGAATATCAt                    >  1:113635/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ttACTGGGTGCGGTGGCGGCCTGGGTCGGACTGGCTATTTTTGGCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATGCCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCGCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCGACAACATCGCCCGTTTCTCGCCAATATCGGTGGCATTAACAAATTGCGCAATGAATATCATGTCt                <  2:560012/250‑1 (MQ=255)
                                                                                                                                                                                                                                     tACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCta               <  1:573140/250‑1 (MQ=255)
                                                                                                                                                                                                                                      aCTGGTGGAGGTGGCGGCCTGGGTCGTACTGGCTATTTGTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGTGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCtat              <  2:249341/250‑1 (MQ=255)
                                                                                                                                                                                                                                         ggtggcggtggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATg           <  1:283526/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ggcggtggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAcc        <  1:150659/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 tggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGAttcg  >  1:20482/1‑248 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
ACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCG  >  NC_000913/4636504‑4636995

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