Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,524,184 C→T A438A (GCC→GCT 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,1840CT98.6% 241.5 / NA 73A438A (GCC→GCTmrcApeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (27/44);  minor base G (1/0);  total (28/45)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGG  >  NC_000913/3523939‑3524430
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
gcCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTAt                                                                                                                                                                                                                                                   <  1:102432/251‑1 (MQ=255)
 ccAATCCGCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATg                                                                                                                                                                                                                                                  <  2:378300/251‑1 (MQ=255)
  cAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCCTGGTGGCTGGCCCAAGTGCCGGAAGTGAACTCGGCGCTGGGGTCGATCAATCCGCAAAACGGTGCTGTTATgg                                                                                                                                                                                                                                                 >  2:329001/1‑251 (MQ=255)
         cAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGt                                                                                                                                                                                                                                          <  1:106766/251‑1 (MQ=255)
            cAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGGGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATGCTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAAGCCGCAAAACGGTGCTGTTATGGCGCTGGTCgt                                                                                                                                                                                                                                       <  2:218017/251‑2 (MQ=255)
             aaGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAACGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGACCTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGt                                                                                                                                                                                                                                      >  1:354418/1‑251 (MQ=255)
                        aTGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGAttt                                                                                                                                                                                                                           <  1:316894/251‑1 (MQ=255)
                                      gTCGACAGGCGCATTGGGTAAGAAAGGCGTTCTCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAGGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGTCGATGCATGGGGGCTGGAACAAGTGCCGAAAGAGCAATCGGCGCAGGTGTCGATCAATCCGCAAAAAGGTGCTGGTATTGCGCTGGTCGGGGGCGTTGAATTCCATCatggcatgt                                                                                                                                                                                                             >  2:136577/1‑243 (MQ=255)
                                                  aTTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTCCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGGCAGGTTGGCGATGCATGGTGGCTGTCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCGGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCa                                                                                                                                                                                                 <  1:36079/251‑1 (MQ=255)
                                                           ggAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCAc                                                                                                                                                                                        >  1:348206/1‑251 (MQ=255)
                                                            gAAGGCGTTCGCTGGGCGCGTCCTTACAGTTCGGATACTCAGCAAGGACCGACGCCGCGTACAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACt                                                                                                                                                                                       >  2:253761/1‑251 (MQ=255)
                                                             aaGGCGTTCGCGGGGCGCTGCCTTACCGTTCGGATACTCAGCAAGGCCCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGGCAGGTTGGAGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCGGCCCAGGCACTg                                                                                                                                                                                      <  2:112643/251‑1 (MQ=255)
                                                                          ggCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGtt                                                                                                                                                                         >  2:134594/1‑251 (MQ=255)
                                                                               gTCCTTACCGTTCGGATACACAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTGTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAAc                                                                                                                                                                    <  2:43730/251‑1 (MQ=255)
                                                                                    tACCGTTCTGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCACATCTGGGTTCGTCAGATTGGCGATGCATGGTGGCTGGCACAAGTGCCGGACGTGAACTCGGCGCTGGGGTCGATCAATCCGCAAAACGGTGCGGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCa                                                                                                                                                                >  2:176693/1‑250 (MQ=255)
                                                                                    tACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCACCATCa                                                                                                                                                                >  2:12027/1‑250 (MQ=255)
                                                                                    tACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATc                                                                                                                                                                 <  2:391383/249‑1 (MQ=255)
                                                                                       cGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCAACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAAcc                                                                                                                                                            >  2:240735/1‑251 (MQ=255)
                                                                                                                  ccgcGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCa                                                                                                                                                                                                 >  2:105429/1‑187 (MQ=255)
                                                                                                                  ccgcGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCa                                                                                                                                                                                                 <  1:105429/187‑1 (MQ=255)
                                                                                                                  ccgcGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATaa                                                                                                                                 <  1:283723/251‑1 (MQ=255)
                                                                                                                   cgcgTAAAGTGACCGCTGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATaa                                                                                                                                 >  2:335726/1‑250 (MQ=255)
                                                                                                                   cgcgTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATaa                                                                                                                                 >  2:209429/1‑250 (MQ=255)
                                                                                                                    gcgTAAAGTGACCTATTTTCTGCAACAGGGTCATCAAATCTGGGTTCGTCAGGTTGGCGATGCATTGTGGCTGGCACAAGTGCCGTAAGTGAACTCGGCGCTGGTGTCGATCAAGCCTCAAATCGGTGGTGTTATGCCGCTGGTCGGTGGCTTTGATTACAATCAGGGCAAGTTTAACCGCTCCCCCCAGGCTCTTCGTCAGGTGGGTTCCAACATCAAACCGTTCCTTTACACCCCGGCGATGGATaaa                                                                                                                                <  2:226968/250‑1 (MQ=255)
                                                                                                                     cgTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAgg                                                                                                                              <  1:12027/251‑1 (MQ=255)
                                                                                                                      gTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGCTGGa                                                                                                                                    >  2:83787/1‑244 (MQ=255)
                                                                                                                      gTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGa                                                                                                                                    <  1:83787/244‑1 (MQ=255)
                                                                                                                       tAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTc                                                                                                                            >  1:270306/1‑251 (MQ=255)
                                                                                                                             gACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGc                                                                                                                      >  1:239210/1‑251 (MQ=255)
                                                                                                                             gACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAAGCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACg                                                                                                                       <  1:329001/250‑1 (MQ=255)
                                                                                                                                cGATGCTCTGAAAACGGGTCAGGAAATCAGGGCTCGCAAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTAGGCGCTGGTGTCGAACAATCCGCAAAACGGTGATCTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGACACACAAGCAATGCGTCAGGTGGGTTCCAACATCAAACCGTATCTCTACCACGCGGCGATGGACACAGGTCTGAGGCTgg                                                                                                                   >  2:279848/1‑251 (MQ=255)
                                                                                                                                    gTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAg                                                                                                               <  2:313124/251‑1 (MQ=255)
                                                                                                                                      tcggCAAACGGGGCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGGGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGGGCTGTTATGGCGCGGGTCGGTGGCTGTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGGCAGGTGGGTGCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAATGTCTGACGCTGGCAAGt                                                                                                              <  2:201949/247‑1 (MQ=255)
                                                                                                                                       cTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTAt                                                                                                            >  2:31406/1‑251 (MQ=255)
                                                                                                                                          cAAACGGGTCAGCAAATCTGGTTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCTCTGCGTCAGGTGGGTTCCAACATCAAACCGTGCCTCTACACCGCGGCGATGGATAAAGGTCTGTCGCTGGCAAGTATGtt                                                                                                         <  2:403609/251‑1 (MQ=255)
                                                                                                                                           aaaCGGGTAGGCAAATCGGGGTTCTGCAGGTTGGCGATGCAGTGGGGCTGGCACGAGTGCGGAAAGTGAACTCGGCGCTGGTGTCGATCAAGCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTGTGAGTTCAATCAGAGCAAGTTTAACCGTGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAACCCGTGCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGtt                                                                                                         <  2:36165/250‑1 (MQ=255)
                                                                                                                                                 gTCAGCAAATCTGGGTACGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCAGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGATGGCAAGTATGTTGAACg                                                                                                    >  2:126080/1‑249 (MQ=255)
                                                                                                                                                    aGCAAATCTGGGTTCGTCAGGTTGGGGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGAtgtg                                                                                               <  2:79661/251‑1 (MQ=255)
                                                                                                                                                       aactCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAAGCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGTCAAGTATGTTGAACGATGTGCCa                                                                                            <  2:24653/248‑1 (MQ=255)
                                                                                                                                                         aTCTGGGTTAGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAAt                                                                                          >  1:325924/1‑251 (MQ=255)
                                                                                                                                                           cTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAAttt                                                                                        >  1:113609/1‑251 (MQ=255)
                                                                                                                                                             gggTTCGCCGGGGTGGCTATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTTGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGCACGATCTGCCAATTtct                                                                                      <  2:335287/251‑1 (MQ=255)
                                                                                                                                                               gTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCg                                                                                    <  1:80670/251‑1 (MQ=255)
                                                                                                                                                                   gTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTGCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTgg                                                                                <  2:224944/251‑1 (MQ=255)
                                                                                                                                                                      aGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGAt                                                                             <  2:348206/251‑1 (MQ=255)
                                                                                                                                                                             cGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTg                                                                      <  1:31406/251‑1 (MQ=255)
                                                                                                                                                                              gATGCATGGTGGCCGGCACAAGTGCCGGACGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGATGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGc                                                                     >  1:128350/1‑251 (MQ=255)
                                                                                                                                                                               aTGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGcc                                                                    <  1:3321/251‑1 (MQ=255)
                                                                                                                                                                                 gCATGTTGGCCGCCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGGTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGGTCCAACATCAAACCTTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCGGGCAAGTATGTTGAACGATGTGCCAATTTCGCTCGGGGATGCAAGTGCCg                                                                   <  2:354418/250‑1 (MQ=255)
                                                                                                                                                                                   aTGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGtt                                                                <  1:176693/251‑1 (MQ=255)
                                                                                                                                                                                       gtggCTGGCACAAGTGGCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTg                                                             <  2:77725/251‑1 (MQ=255)
                                                                                                                                                                                          gctggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACt                                                          <  1:63666/251‑1 (MQ=255)
                                                                                                                                                                                             ggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTg                                                         <  1:358458/249‑1 (MQ=255)
                                                                                                                                                                                                acaAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCGGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGc                                                    <  2:325924/251‑1 (MQ=255)
                                                                                                                                                                                                    gTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCgaa                                                <  2:321191/251‑1 (MQ=255)
                                                                                                                                                                                                     tGCCTGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTg                                                                                                                          <  1:29541/176‑1 (MQ=255)
                                                                                                                                                                                                     tGCCGGAAGTGAACTCGGCGCTGGTGTAGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTg                                                                                                                          >  2:29541/1‑176 (MQ=255)
                                                                                                                                                                                                      gcctgAAGTGAACTCGGCGCTGGTGTCGATCAAACCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGCTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATGAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTTCCAATTTCTCGCTGGGATGCAATTGCCGGTTCTGACTGTCAGCCgaaga                                              <  1:371279/247‑1 (MQ=255)
                                                                                                                                                                                                       ccGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCACTTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCgaaga                                              >  1:299663/1‑250 (MQ=255)
                                                                                                                                                                                                               gAACTCGGCTCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACAGCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTTACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTcacca                                      <  1:50500/250‑1 (MQ=255)
                                                                                                                                                                                                                              gtCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTc                      <  1:209429/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCCACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAAtt                 >  1:74271/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAAtt                 >  1:45226/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aaTCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCg               <  1:98488/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aaTCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGGGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCg               >  2:101581/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aaaCGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTc       <  1:319345/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCGTTGCTGTTATGGTGCTGGTCGGTTGCTTTGATTTCCAGCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTCCACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGGTGTGCCAATTTCGCGCTGGGATTCAAGTCCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGTAGTATGCTGGTCCAATTCGCTTACGTCa      <  2:326896/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTc       >  1:15153/1‑250 (MQ=255)
                                                                                                                                                                                                                                              aaCGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCa      <  1:134594/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aCGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTATGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTAACCACCGCAGTATGCTGGTCCAATTCGCTTAAGTc       >  1:398754/1‑249 (MQ=255)
                                                                                                                                                                                                                                                cGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAgg    <  1:216818/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAgggg  <  1:240735/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTGCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAgggg  <  1:395837/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGG  >  NC_000913/3523939‑3524430

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