Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,787,513 T→C L226L (TTG→CTA)  envC → murein hydrolase activator EnvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,787,5130TC97.8% 134.3 / ‑6.3 45L226L (TTG→CTA) envCmurein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (20/24);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAGAGCGAGCAACAAACGCTGTTATATGAGCAGCGCGCCCAACAGGCGAAACTGACTCAGGCGCTGAACGAGCGTAAAAAGACGCTGGCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAG‑TTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGC  >  NC_000913/3787291‑3787746
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            
ttCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATccc                                                                                                                                                                                                                   <  2:408784/250‑1 (MQ=255)
   ggCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAAc                                                                                                                                                                                                                          <  2:18480/240‑1 (MQ=255)
    gCTATCTCAACCAGGCGAGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGCGTTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGGGATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCACCCTGCAAGGTCAGCAGCAG‑CCAAGTGAGCTGCGTGCCAACGAATCCCGTCTg                                                                                                                                                                                                              <  2:357750/251‑1 (MQ=255)
         ctcAACCAGGCGCAACAAGAACCCATTGCACAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAAACTTCCATCCCGCAAGGGCAGCAGCAG‑CTAAGTGAACTGCGTGACAACGAAACCCGTCTGCGTa                                                                                                                                                                                                          >  2:418400/1‑250 (MQ=255)
                       aCAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCTGCAGCAG‑CTAAGTGAACGGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTg                                                                                                                                                                                           <  2:193804/251‑1 (MQ=255)
                        aaaGAATCCATTGCCCAGGTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTGATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAATACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTTAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGc                                                                                                                                                                                          <  1:279612/250‑1 (MQ=255)
                         aaGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCg                                                                                                                                                                                         <  2:241088/251‑1 (MQ=255)
                                cATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAACATTGCCCGTGCGG‑AAACcg                                                                                                                                                                                 >  1:174083/1‑251 (MQ=255)
                                              aGCAAACGCGTGAAGAAGTCGCCATGCAGCGAGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTACATCCCGCAAGGTCAGCAGCAG‑CTAAGGGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGGAAGCGCGt                                                                                                                                                                   >  1:154021/1‑251 (MQ=255)
                                               gCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCATCCAACCAGCCACGGCAGCAGCAG‑CTACGTGAACTGCGGGGCAACGGATCCCGACCGCGTCACAGCCATGCCCGTTCGGGAAGC‑GCGGCGGAAGcgcggg                                                                                                                                                                  >  2:313958/1‑249 (MQ=255)
                                                cAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGGATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACCGCATTTCCCGCGCGG‑AAGCCGCGGCGAAAGcgcgggg                                                                                                                                                                 >  2:252267/1‑248 (MQ=255)
                                                cAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGGGCCCACGAATCCCGTCTGCGGAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGCAAGcgcgggc                                                                                                                                                                 >  2:89589/1‑248 (MQ=255)
                                                 aaaCGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCCTCCAGCAAGGGCAGCAGCAG‑CTAAGAGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCACGTGc                                                                                                                                                                 >  2:124466/1‑250 (MQ=255)
                                                      cgTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCCTCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACg                                                                                                                                                           >  2:169004/1‑251 (MQ=255)
                                                      cgTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑CAGCCGCGGCGAAAGCGCGTGCAGAACg                                                                                                                                                           >  1:45545/1‑251 (MQ=255)
                                                       gTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAACAGACGCTGGCGGGGCTGGTATCTTCCATTCAGCAAGGGCAGCACCAG‑CTAAATTAAATGCGTGCCAACGAATACCGTTTGCGTAACAGCATTGCCCGTGTGG‑AAGCCGCGCGGAAAGCGCGTgtggaacgg                                                                                                                                                          >  2:107404/1‑243 (MQ=255)
                                                         gaagaaGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGATGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATATTTCCTCCAGCAAGGGCAACCGCAGCTTAAG‑GAACTGCGGGCCAACGAATCCCGTCTGCGTAAAAGCATTGTCCCTGCGG‑AAACCGCGGCGAAAGCGCGGGCAGAACg                                                                                                                                                           >  2:393278/1‑248 (MQ=255)
                                                         gaagaaGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGAAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGGCAGCAGCAG‑CTAAGTGAACTGCGGGCCAACGAATCCCGTCTGCGTAACCGCATTGCCCGTGCGG‑AAACCGCGGCGGAAGCGCGTGCAAAACg                                                                                                                                                           >  1:409280/1‑248 (MQ=255)
                                                          aagaagTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGaa                                                                                                                                                       <  1:393278/251‑1 (MQ=255)
                                                            gaagTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGc                                                                                                                                                     <  1:332296/251‑1 (MQ=255)
                                                                 cGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTg                                                                                                                                                <  2:99533/251‑1 (MQ=255)
                                                                    cATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAgg                                                                                                                                             <  1:185236/251‑1 (MQ=255)
                                                                                 gAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCGCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑tt                                                                                                                                 <  2:232203/249‑1 (MQ=255)
                                                                                       gAAGAGAAACAAAGCGAGCAACAAACGCTGTTAGATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CGAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCg                                                                                                                         <  1:293532/251‑1 (MQ=255)
                                                                                               acaaagcgagcAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCcgc                                                                                                                                                                                <  1:245533/185‑1 (MQ=255)
                                                                                                  aagcgagcAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCCAGGTCAGCAGCAG‑CTAAACGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGGGCAAAACGAGAAGCACGTGAGGACCAGGCGG‑GTCGCGACCGcccggaaaa                                                                                                                >  1:82437/2‑242 (MQ=255)
                                                                                                    gcgagcAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCCGGAATCTTCCCTCCAGCAAGGTCCGCAGCAG‑CTAAGTGAAATGCGTGCCAACGAATCCCGTCTGCGTAACAGCCTTGCCCGTGCGG‑GAGCCGCGGCGAAAACGCGTGCAGAACGGGAAGCACGGGAGGCCCAGGCGGGTTCGGGAcccccgaaaagagcg                                                                                                             >  2:348655/1‑238 (MQ=255)
                                                                                                       agcAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTTCATCCACCAAAGTCAGCATCCA‑CCAAGTGAACTGCGGGCCCAAGAATCCCGTCCGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCTGCGAAAGCGCGCGCCGAACGGGCAGCACGCGAGGCCCAGGCTG‑TTtgcggccccccgacaaaagccacc                                                                                                         >  1:167138/1‑231 (MQ=255)
                                                                                                       agcAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGAcg                                                                                                         <  1:416484/251‑1 (MQ=255)
                                                                                                              aaCGCGGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAg                                                                                                 <  2:338698/251‑1 (MQ=255)
                                                                                                                 gCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGca                                                                                              <  1:348655/251‑1 (MQ=255)
                                                                                                                    gTTATAGGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATGGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGcacc                                                                                            <  2:192064/250‑1 (MQ=255)
                                                                                                                        aatGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCCGCCG‑CTCAGTGAACTGCGTGCCAACGAATCCCGGCTGCGGAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTTCCGAACGAGAAGAACGTGAGGCCCCGGAGG‑GTCGCGAACGCCCGAAAGAAGGGGAGCCC‑AAAGGcaccccccc                                                                                       >  2:10733/2‑246 (MQ=255)
                                                                                                                              cagcGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAacc                                                                                 <  1:344913/251‑1 (MQ=255)
                                                                                                                                  gcgcCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAaccgacc                                                                             <  1:13138/251‑1 (MQ=255)
                                                                                                                                          cAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAACCGACCGAAAGCGa                                                                     <  1:145293/251‑1 (MQ=255)
                                                                                                                                              cGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCTTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAACCGACCGAAAGCGaaaaa                                                                 <  1:407864/251‑1 (MQ=255)
                                                                                                                                                    tGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCCG‑CTAAGGGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCCGGCGG‑GTCGCGACCCCCCGAAAGAAACGACGCGC‑CAAGGCACCCCCTACAAACCCACCGAAAACGGAAAAACGCTg                                                           >  1:413677/1‑251 (MQ=255)
                                                                                                                                                      aCTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGGATCTTCAATCCAGCAAGGTCAGCCGCAG‑CTAAGTGAACTGCGTGGCAAAGAATCCCGTCTGCGAAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑GTCCGGAACGCCAGAAAAAAACGAGGCGCCAAAGGC‑CCACCTACAAACCGCACGAAAAGGAAAAATCGCTGAt                                                         >  1:138582/1‑251 (MQ=255)
                                                                                                                                                        ccAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGt                                                       <  2:155294/250‑1 (MQ=255)
                                                                                                                                                                   cgacgctctgccgatcGCTGGCGGGGCTGGAATCTTCCATCCATCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAAGCCCGGCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCTCTGAGGTCCCgaactggt                                            <  2:295339/236‑8 (MQ=255)
                                                                                                                                                                                 aCGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTTAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGAACTGGTGGCCTGGGTGcgcc                              <  1:141324/251‑1 (MQ=255)
                                                                                                                                                                                    cTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAAATGCGTGCCAACGAATCCCGTCTGCGTAACCGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCCCGGGAGGCCCCGGCGG‑TTCGCGACCGTCCGAAAGAAGAGGCGCGC‑AAAGGCACCACCTCCCAACCGACCGAAAGCGAAAAATGGCTGGTGTCCCgccctggt                                            >  1:295339/1‑227 (MQ=255)
                                                                                                                                                                                                     ttCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCACAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCCCCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGAACTGGTGGCCTGGGTGCGCCGCGTGGTCAGGCATTCTGGc          >  2:154654/1‑251 (MQ=255)
                                                                                                                                                                                                             aGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAACGCGCGTGCAGAACGAGAAGCACGTGAGACCCAGGCGG‑TTCGCGCCCGCCAGAAAGCAGCGACGCGC‑AAAGGCACCACCTACACACCGACCGAAAGCGAAAAATCGCTGATGTCCCGAACTGTTGGCCTGGGTGCGCCGCGTGGTCAGGCCTTCTGGCCGGTTcgc  >  1:438730/1‑251 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            
TTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAGAGCGAGCAACAAACGCTGTTATATGAGCAGCGCGCCCAACAGGCGAAACTGACTCAGGCGCTGAACGAGCGTAAAAAGACGCTGGCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAG‑TTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGC  >  NC_000913/3787291‑3787746

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