Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,787,509 A→G Q224Q (CAA→CAG envC → murein hydrolase activator EnvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,787,5090AG97.7% 114.5 / ‑6.1 45Q224Q (CAA→CAGenvCmurein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (19/24);  minor base T (1/0);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

CAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAGAGCGAGCAACAAACGCTGTTATATGAGCAGCGCGCCCAACAGGCGAAACTGACTCAGGCGCTGAACGAGCGTAAAAAGACGCTGGCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGC  >  NC_000913/3787264‑3787746
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
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)
                           ttCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATccc                                                                                                                                                                                                                   <  2:408784/250‑1 (MQ=255)
                              ggCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAAc                                                                                                                                                                                                                          <  2:18480/240‑1 (MQ=255)
                              ggCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCGGGGATCTTCCATCCAGCAAGGTcagcagcagccaaatgaactgcgtgcccac                                                                                                                                                                                                                           >  1:18480/1‑219 (MQ=255)
                               gCTATCTCAACCAGGCGAGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGCGTTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGGGATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCACCCTGCAAGGTCAGCAGCAGCCAAGTGAGCTGCGTGCCAACGAATCCCGTCTg                                                                                                                                                                                                              <  2:357750/251‑1 (MQ=255)
                                    ctcAACCAGGCGCAACAAGAACCCATTGCACAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAAACTTCCATCCCGCAAGGGCAGCAGCAGCTAAGTGAACTGCGTGACAACGAAACCCGTCTGCGTa                                                                                                                                                                                                          >  2:418400/1‑250 (MQ=255)
                                                  aCAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCTGCAGCAGCTAAGTGAACGGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTg                                                                                                                                                                                           <  2:193804/251‑1 (MQ=255)
                                                   aaaGAATCCATTGCCCAGGTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTGATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAATACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTTAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGc                                                                                                                                                                                          <  1:279612/250‑1 (MQ=255)
                                                    aaGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCg                                                                                                                                                                                         <  2:241088/251‑1 (MQ=255)
                                                           cATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAACATTGCCCGTGCGGAAACcg                                                                                                                                                                                 >  1:174083/1‑251 (MQ=255)
                                                                         aGCAAACGCGTGAAGAAGTCGCCATGCAGCGAGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTACATCCCGCAAGGTCAGCAGCAGCTAAGGGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGGAAGCGCGt                                                                                                                                                                   >  1:154021/1‑251 (MQ=255)
                                                                          gCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCATCCAACCAGCCACGGCAGCAGCAGCTACGTGAACTGCGGGGCAACGGATCCCGACCGCGTCACAGCCATGCCCGTTCGGGAAGCGCGGCGGAAGcgcggg                                                                                                                                                                  >  2:313958/1‑249 (MQ=255)
                                                                           cAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGGATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACCGCATTTCCCGCGCGGAAGCCGCGGCGAAAGcgcgggg                                                                                                                                                                 >  2:252267/1‑248 (MQ=255)
                                                                           cAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGGGCCCACGAATCCCGTCTGCGGAACAGCATTGCCCGTGCGGAAGCCGCGGCGCAAGcgcgggc                                                                                                                                                                 >  2:89589/1‑248 (MQ=255)
                                                                            aaaCGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCCTCCAGCAAGGGCAGCAGCAGCTAAGAGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCACGTGc                                                                                                                                                                 >  2:124466/1‑250 (MQ=255)
                                                                                 cgTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCCTCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACg                                                                                                                                                           >  2:169004/1‑251 (MQ=255)
                                                                                 cgTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGCAGCCGCGGCGAAAGCGCGTGCAGAACg                                                                                                                                                           >  1:45545/1‑251 (MQ=255)
                                                                                  gTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAACAGACGCTGGCGGGGCTGGTATCTTCCATTCAGCAAGGGCAGCACCAGCTAAATTAAATGCGTGCCAACGAATACCGTTTGCGTAACAGCATTGCCCGTGTGGAAGCCGCGCGGAAAGCGCGTgtggaacgg                                                                                                                                                          >  2:107404/1‑243 (MQ=255)
                                                                                    gaagaaGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGATGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATATTTCCTCCAGCAAGGGCAACCGCAGCTTAAGGAACTGCGGGCCAACGAATCCCGTCTGCGTAAAAGCATTGTCCCTGCGGAAACCGCGGCGAAAGCGCGGGCAGAACg                                                                                                                                                           >  2:393278/1‑248 (MQ=255)
                                                                                    gaagaaGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGAAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGGCAGCAGCAGCTAAGTGAACTGCGGGCCAACGAATCCCGTCTGCGTAACCGCATTGCCCGTGCGGAAACCGCGGCGGAAGCGCGTGCAAAACg                                                                                                                                                           >  1:409280/1‑248 (MQ=255)
                                                                                     aagaagTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGaa                                                                                                                                                       <  1:393278/251‑1 (MQ=255)
                                                                                       gaagTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGc                                                                                                                                                     <  1:332296/251‑1 (MQ=255)
                                                                                            cGCCATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTg                                                                                                                                                <  2:99533/251‑1 (MQ=255)
                                                                                               cATGCAGCGTGCTGAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAgg                                                                                                                                             <  1:185236/251‑1 (MQ=255)
                                                                                                            gAACTGGAAGAGAAACAAAGCGAGCAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCGCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGtt                                                                                                                                 <  2:232203/249‑1 (MQ=255)
                                                                                                                  gAAGAGAAACAAAGCGAGCAACAAACGCTGTTAGATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCGAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCg                                                                                                                         <  1:293532/251‑1 (MQ=255)
                                                                                                                          acaaagcgagcAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCcgc                                                                                                                                                                                <  1:245533/185‑1 (MQ=255)
                                                                                                                             aagcgagcAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCCAGGTCAGCAGCAGCTAAACGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGGGCAAAACGAGAAGCACGTGAGGACCAGGCGGGTCGCGACCGcccggaaaa                                                                                                                >  1:82437/2‑242 (MQ=255)
                                                                                                                               gcgagcAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCCGGAATCTTCCCTCCAGCAAGGTCCGCAGCAGCTAAGTGAAATGCGTGCCAACGAATCCCGTCTGCGTAACAGCCTTGCCCGTGCGGGAGCCGCGGCGAAAACGCGTGCAGAACGGGAAGCACGGGAGGCCCAGGCGGGTTCGGGAcccccgaaaagagcg                                                                                                             >  2:348655/1‑238 (MQ=255)
                                                                                                                                  agcAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTTCATCCACCAAAGTCAGCATCCACCAAGTGAACTGCGGGCCCAAGAATCCCGTCCGCGTAACAGCATTGCCCGTGCGGAAGCCGCTGCGAAAGCGCGCGCCGAACGGGCAGCACGCGAGGCCCAGGCTGTTtgcggccccccgacaaaagccacc                                                                                                         >  1:167138/1‑231 (MQ=255)
                                                                                                                                  agcAACAAACGCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGAcg                                                                                                         <  1:416484/251‑1 (MQ=255)
                                                                                                                                         aaCGCGGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAg                                                                                                 <  2:338698/251‑1 (MQ=255)
                                                                                                                                            gCTGTTATATGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGca                                                                                              <  1:348655/251‑1 (MQ=255)
                                                                                                                                               gTTATAGGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATGGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGcacc                                                                                            <  2:192064/250‑1 (MQ=255)
                                                                                                                                                   aatGAACAGCGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCCGCCGCTCAGTGAACTGCGTGCCAACGAATCCCGGCTGCGGAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTTCCGAACGAGAAGAACGTGAGGCCCCGGAGGGTCGCGAACGCCCGAAAGAAGGGGAGCCCAAAGGcaccccccc                                                                                       >  2:10733/2‑246 (MQ=255)
                                                                                                                                                         cagcGCGCCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAacc                                                                                 <  1:344913/251‑1 (MQ=255)
                                                                                                                                                             gcgcCCAACAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAaccgacc                                                                             <  1:13138/251‑1 (MQ=255)
                                                                                                                                                                     cAGGCGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGa                                                                     <  1:145293/251‑1 (MQ=255)
                                                                                                                                                                         cGAAACTGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCTTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGaaaaa                                                                 <  1:407864/251‑1 (MQ=255)
                                                                                                                                                                               tGACTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCCGCTAAGGGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCCGGCGGGTCGCGACCCCCCGAAAGAAACGACGCGCCAAGGCACCCCCTACAAACCCACCGAAAACGGAAAAACGCTg                                                           >  1:413677/1‑251 (MQ=255)
                                                                                                                                                                                 aCTCAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGGATCTTCAATCCAGCAAGGTCAGCCGCAGCTAAGTGAACTGCGTGGCAAAGAATCCCGTCTGCGAAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGGTCCGGAACGCCAGAAAAAAACGAGGCGCCAAAGGCCCACCTACAAACCGCACGAAAAGGAAAAATCGCTGAt                                                         >  1:138582/1‑251 (MQ=255)
                                                                                                                                                                                   ccAGGCTCTGAATGAGCGTAAAAAGACGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGt                                                       <  2:155294/250‑1 (MQ=255)
                                                                                                                                                                                              cgacgctctgccgatcGCTGGCGGGGCTGGAATCTTCCATCCATCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAAGCCCGGCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCTCTGAGGTCCCgaactggt                                            <  2:295339/236‑8 (MQ=255)
                                                                                                                                                                                                            aCGCTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTTAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGAACTGGTGGCCTGGGTGcgcc                              <  1:141324/251‑1 (MQ=255)
                                                                                                                                                                                                               cTGGCGGGGCTGGAATCTTCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAAATGCGTGCCAACGAATCCCGTCTGCGTAACCGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCCCGGGAGGCCCCGGCGGTTCGCGACCGTCCGAAAGAAGAGGCGCGCAAAGGCACCACCTCCCAACCGACCGAAAGCGAAAAATGGCTGGTGTCCCgccctggt                                            >  1:295339/1‑227 (MQ=255)
                                                                                                                                                                                                                                ttCCATCCAGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCACAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCCCCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGAACTGGTGGCCTGGGTGCGCCGCGTGGTCAGGCATTCTGGc          >  2:154654/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aGCAAGGTCAGCAGCAGCTAAGTGAACTGCGTGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAACGCGCGTGCAGAACGAGAAGCACGTGAGACCCAGGCGGTTCGCGCCCGCCAGAAAGCAGCGACGCGCAAAGGCACCACCTACACACCGACCGAAAGCGAAAAATCGCTGATGTCCCGAACTGTTGGCCTGGGTGCGCCGCGTGGTCAGGCCTTCTGGCCGGTTcgc  >  1:438730/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
CAGCGTGGACAGCGTTTACAGGCTTATTTCGGCTATCTCAACCAGGCGCGACAAGAAACCATTGCCCAGTTGAAGCAAACGCGTGAAGAAGTCGCCATGCAGCGTGCTGAACTGGAAGAGAAACAGAGCGAGCAACAAACGCTGTTATATGAGCAGCGCGCCCAACAGGCGAAACTGACTCAGGCGCTGAACGAGCGTAAAAAGACGCTGGCAGGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCGACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGC  >  NC_000913/3787264‑3787746

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