Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,787,449 G→T A204A (GCG→GCT)  envC → murein hydrolase activator EnvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,787,4490GT98.1% 169.1 / NA 53A204A (GCG→GCT) envCmurein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (23/29);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCGAGCATACCGGTATTCAGCTGATTCTCAGCGGTGAAGAAAGCCAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAGAGCGAGCAACAAACGCTGTTATATGAGCAGCGCGCCCAACAGGCGAAACTGACTCAGGCGCTGAACGAGCGTAAAAAGACGCTGGCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAG‑TTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGT  >  NC_000913/3787220‑3787693
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
gCGAACATACCTGTATTCAGCTGATTCTCAGCGGTGAAGAAAGCCAGCGTGGACAGCGTTTACAGGCTTATTGCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGTCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACg                                                                                                                                                                                                                                     <  1:185295/251‑1 (MQ=255)
  gaacATACCGGTATTCAACTGATTCTCAGCGGTGAAGAAAGCCAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCt                                                                                                                                                                                                                                   >  1:134608/4‑251 (MQ=255)
   aacATACAGGTATTCAGCTGATTCTCAGCGGTGAAGAAAGCCAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAATCGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGTCTCTGAATGAGCGTAAAAAGACGCTg                                                                                                                                                                                                                                  <  2:333125/249‑1 (MQ=255)
      aTACCGGTATTCAGCTGATTCTCAGCGGTGAAGAAAGCCAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGcg                                                                                                                                                                                                                               <  1:354771/251‑2 (MQ=255)
                            cAGCGGTGAAGAAAGCCAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTg                                                                                                                                                                                                                                                       >  1:96816/1‑205 (MQ=255)
                            cAGCGGTGAAGAAAGCCAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTg                                                                                                                                                                                                                                                       <  2:96816/205‑1 (MQ=255)
                            cAGCGGTGAAGAAAGCCAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTtccatcc                                                                                                                                                                                                            <  1:94086/248‑1 (MQ=255)
                                            cAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTcagcagcagct                                                                                                                                                                                         >  1:241088/1‑249 (MQ=255)
                                                  ggACAGCGTTTACAGGCTTCTTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGACCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGCATCTTCCATCCAACAAGGAcagcagcagccaagcg                                                                                                                                                                                    >  2:353806/1‑243 (MQ=255)
                                                      aGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTcagcagcagctaagtgaactg                                                                                                                                                                               <  1:92776/251‑13 (MQ=255)
                                                         gTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTcagcagcagctaagtgaactg                                                                                                                                                                               <  1:155836/248‑13 (MQ=255)
                                                             aCAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTcagcagcagctaagtgaactgcgtgcc                                                                                                                                                                         <  1:103042/250‑19 (MQ=255)
                                                                       ttCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAACGAATccc                                                                                                                                                              <  2:408784/250‑1 (MQ=255)
                                                                          ggCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAG‑CTAAGTGAACTGCGTGCCAAc                                                                                                                                                                     <  2:18480/240‑1 (MQ=255)
                                                                          ggCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCGGGGATCTTCCATCCAGCAAGGTcagcagcagccaaatgaactgcgtgcccac                                                                                                                                                                      >  1:18480/1‑219 (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)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
GCGAGCATACCGGTATTCAGCTGATTCTCAGCGGTGAAGAAAGCCAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAGAGCGAGCAACAAACGCTGTTATATGAGCAGCGCGCCCAACAGGCGAAACTGACTCAGGCGCTGAACGAGCGTAAAAAGACGCTGGCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAG‑TTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGG‑AAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TTCGCGACCGCCAGAAAGAAGCGACGCGC‑AAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGT  >  NC_000913/3787220‑3787693

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