Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,524,603 C→T A578V (GCG→GTG)  mrcA → peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,524,6030CT96.4% 172.6 / ‑2.1 56A578V (GCG→GTG) mrcApeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (25/29);  total (26/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 = 7.72e-01

GCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGCAGCTGGA  >  NC_000913/3524356‑3524834
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       
gCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGCACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACCTTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCCATGCAGGTGGCGCGCCGCTACGCGGt                                                                                                                                                                                                                                      >  2:394065/1‑251 (MQ=255)
                aCTGGCAGCCGAAGAACTCTCCACCGCAGTATGCTGTTCCAATTCGCTTACGTCAGGGGCCGGGTCAGTCGAAACACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACggcggc                                                                                                                                                                                                                      <  1:286728/251‑1 (MQ=255)
                   ggCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTTCGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGGCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTc                                                                                                                                                                                                                   <  1:394065/251‑1 (MQ=255)
                      aGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCtg                                                                                                                                                                                                                <  1:5934/251‑1 (MQ=255)
                         cGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCtggtg                                                                                                                                                                                                             >  1:67955/1‑251 (MQ=255)
                                    ccaccGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGtt                                                                                                                                                                                                  <  1:230681/251‑1 (MQ=255)
                                     caccGCAGTATGCTGGTCCAAGTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGTGCGTCGACTCCGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGttt                                                                                                                                                                                                 <  2:171999/251‑1 (MQ=255)
                                        agCAGTATGCTGGTCGAATTCGCTTCCGTCAGGGGCTGGGTCAATCGCACAACGTGGAAATGGTACGCGCAATGAGGCCGATGGGCGTCGACTACGCTGCAGAATATCTGCAAAGCTTCGGCATCCCGGCACACAAAAATGGCCAAAGCGAATCGCTGGAGCTGGGTACAGCGTCCTTAACCCCAAAGCAGGTGGCGCGAGGCTACGTGCACATGGCGAACGGCGGATTACTggaggaagcggggttaatc                                                                                                                                                                                              >  2:30098/2‑238 (MQ=255)
                                          cAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATc                                                                                                                                                                                              >  1:101044/1‑249 (MQ=255)
                                          cAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAg                                                                                                                                                                                            <  1:381324/251‑1 (MQ=255)
                                                gCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAAt                                                                                                                                                                                      >  2:135865/1‑251 (MQ=255)
                                                   ggTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGGGGGCGATGGGCGTCGACTACGCTTCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTAGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAAtt                                                                                                                                                                                     <  2:106668/249‑1 (MQ=255)
                                                   ggTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTg                                                                                                                                                                                    >  2:251485/1‑250 (MQ=255)
                                                      ccAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACGACGCTGCAGAATATCTGCAACGCGTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGaaaa                                                                                                                                                                                <  1:180842/251‑1 (MQ=255)
                                                       cAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACCTCTGCTTGCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAAc                                                                                                                                                                               <  1:330979/251‑1 (MQ=255)
                                                                tACGTCAGGGGCTGGGTAAGTCGACCAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGACTATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCCTGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCCAAATTGAACACGATCAGggc                                                                                                                                                                      >  2:389242/1‑251 (MQ=255)
                                                                     cAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGgcggcg                                                                                                                                                                  <  2:113886/250‑1 (MQ=255)
                                                                         ggcggggTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGAtt                                                                                                                                                             <  1:37899/247‑1 (MQ=255)
                                                                               gTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGaa                                                                                                                                                       <  1:106673/251‑1 (MQ=255)
                                                                               gTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGa                                                                                                                                                        >  2:295200/1‑250 (MQ=255)
                                                                                 cAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGc                                                                                                                                                     >  2:16619/1‑251 (MQ=255)
                                                                                  attCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCg                                                                                                                                                    <  1:374352/249‑1 (MQ=255)
                                                                                     cgacaaaCGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGaa                                                                                                                                                  >  2:349324/5‑250 (MQ=255)
                                                                                       aacaaCGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCCTGGTTTCTCAGCAAAATTGAAAACGATCAGGACGGCGTGATTTTctcatcgaaccc                                                                                                                                               >  2:328622/4‑240 (MQ=255)
                                                                                        aacaCGTGGGGATGGTACGCTCAATGCGGACGATGGGCGTCGACTACACTGCAGAATATCTGCACCGCTTCGGCTTCCCTGCTCAAAACATTGTCCACACCGAATCGCTGGCGCTGGTTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGTCCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCg                                                                                                                                              <  2:182262/248‑1 (MQ=255)
                                                                                           aCGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGaa                                                                                                                                            >  1:189407/1‑250 (MQ=255)
                                                                                             gtggtgATGGTACGCGCAATGCGGGGGATGTGTGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTTGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGGGTGAGTTTCGAAGCGAAACCGAAAGt                                                                                                                                         <  2:380830/251‑1 (MQ=255)
                                                                                               ggtgATGGTACGCGCAATGCGGGCGATTGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTAATTTTCGAATCGAAACCGAAAGTa                                                                                                                                        <  2:76905/250‑1 (MQ=255)
                                                                                                     ggTACGCTCAATGCGGGCGATGGTCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACTCCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGcc                                                                                                                                 <  2:334932/251‑1 (MQ=255)
                                                                                                             cAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGAttt                                                                                                                                                            >  1:247658/1‑216 (MQ=255)
                                                                                                             cAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGAttt                                                                                                                                                            <  2:247658/216‑1 (MQ=255)
                                                                                                              aaTGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCg                                                                                                                        >  2:105838/1‑251 (MQ=255)
                                                                                                                    ggggATGGGCGTCGACTTCTCTGCAGAATATCTGCAACGCTTCGGCTGCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCTCGGCTACGTGGTCATGGCGAACGGCGGCTTACTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAGCCGAAAGTATCCTGCCCGGAATGCTATATTc                                                                                                                  <  2:168379/248‑1 (MQ=255)
                                                                                                                          gggCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCGACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTg                                                                                                             <  1:185080/250‑1 (MQ=255)
                                                                                                                          gggCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCCAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTg                                                                                                             >  2:288609/1‑250 (MQ=255)
                                                                                                                              gTCGACGACGCAGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCTCCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTa                                                                                                        <  2:382440/251‑1 (MQ=255)
                                                                                                                                   cTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGATCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTg                                                                                                   >  1:21617/1‑251 (MQ=255)
                                                                                                                                                      cAACGCGTCGGCTTCCCGGCACACAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGTGGCTACGTGGTCATGGCGAACGGCGGCGTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCTCGGAATGCGATATTCCGGTGGTTTACGGTGATACGCAGAAATCGAACGt                                                                                <  2:396604/251‑1 (MQ=255)
                                                                                                                                                           cTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTgg                                                                           >  2:177415/1‑251 (MQ=255)
                                                                                                                                                                  ttCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATa                                                                     >  2:226132/1‑250 (MQ=255)
                                                                                                                                                                    cccGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCATCCGTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAATCGATCAGGGCGGCGTGATTTTCGAAGCGAGACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGAAAAATAAc                                                                   <  2:189915/250‑1 (MQ=255)
                                                                                                                                                                       ggCACAAACCATTGTCCCCACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTAACCCAAATCCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGAACCGTGGTTTATCAGCAAAATTGAAAACGCTCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGAGATATTCCGGTGATTTACGGTGCTACGCAGAAATCGAACGAGCTGGAAAATAACGATg                                                               >  2:35766/1‑251 (MQ=255)
                                                                                                                                                                          acaAAACATTGTCCACACCGATTCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTCCGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGTAAAATTGAAAATGATCCCGGCGGCGTGATTTTCGAAGCTAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGGGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAATTAACGATGTTg                                                            >  2:41232/1‑251 (MQ=255)
                                                                                                                                                                                            cGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTccc                                          >  1:356408/1‑251 (MQ=255)
                                                                                                                                                                                             gAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGAGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGCTTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGTAAATAACGATGTTGCAGATGTCGCTATCTCCcg                                         >  1:4944/1‑251 (MQ=255)
                                                                                                                                                                                              aaTCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCcgc                                        <  1:16619/251‑1 (MQ=255)
                                                                                                                                                                                              aaTCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCcgc                                        >  1:234346/1‑251 (MQ=255)
                                                                                                                                                                                                               tCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGttt                        <  2:101044/250‑1 (MQ=255)
                                                                                                                                                                                                                cAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTc                       <  1:105838/250‑1 (MQ=255)
                                                                                                                                                                                                                  gCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGt                    >  2:326526/1‑251 (MQ=255)
                                                                                                                                                                                                                  gCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGt                    <  2:67955/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccTTCACCCCAATGCAGGGGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGATTCCTGGTGGACCCGTGGTTTAGCAGCAAAATTGAAAACGATCAGGGAGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTTCTGGATAATACCGAAGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCa                <  2:332096/251‑1 (MQ=255)
                                                                                                                                                                                                                          cACCCCAATGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGACTGTTTCTGTACCAATGc            >  2:401227/1‑251 (MQ=255)
                                                                                                                                                                                                                                aaTGCAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGAGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGCAGc      >  1:287031/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCTAGCAGCAGAATGTTTCTGTACCAATGCCGCAGCTGGa  <  1:18421/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cAGGTGGCGCGCGGCTACGTGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGCAGCTGGa  <  1:303676/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       
GCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGGCTGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCGGAATGCGATATTCCGGTGATTTACGGTGATACGCAGAAATCGAACGTGCTGGAAAATAACGATGTTGAAGATGTCGCTATCTCCCGCGAGCAGCAGAATGTTTCTGTACCAATGCCGCAGCTGGA  >  NC_000913/3524356‑3524834

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