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.8% 270.4 / NA 85A438A (GCC→GCTmrcApeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/51);  minor base G (1/0);  total (34/51)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGG  >  NC_000913/3523939‑3524430
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
gcCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTAt                                                                                                                                                                                                                                                   >  1:333757/1‑251 (MQ=255)
        ccAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTTGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAAGCCGCAAAACGGTGCTGTTATGGCGCTgg                                                                                                                                                                                                                                           <  1:454104/250‑1 (MQ=255)
         cAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGt                                                                                                                                                                                                                                          <  1:94318/251‑1 (MQ=255)
            cAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACGGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCgg                                                                                                                                                                                                                                       <  1:514236/251‑1 (MQ=255)
            cAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCgg                                                                                                                                                                                                                                       <  1:128108/251‑1 (MQ=255)
               gCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTgg                                                                                                                                                                                                                                    <  1:558575/251‑1 (MQ=255)
                    ggCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTg                                                                                                                                                                                                                               >  1:327081/1‑251 (MQ=255)
                        aTGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGAttt                                                                                                                                                                                                                           >  1:340533/1‑251 (MQ=255)
                        aTGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGAttt                                                                                                                                                                                                                           <  1:538915/251‑1 (MQ=255)
                          gCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTc                                                                                                                                                                                                                          >  1:132930/1‑250 (MQ=255)
                          gCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTc                                                                                                                                                                                                                          <  1:413370/250‑1 (MQ=255)
                               cGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCag                                                                                                                                                                                                                    <  1:261378/251‑1 (MQ=255)
                                ggACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATc                                                                                                                                                                                                                      >  1:250775/1‑248 (MQ=255)
                                   cggttCGACCGTCGCTGTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTGCTGCAAACGGGTCAGCAAATCGGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCGGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCa                                                                                                                                                                                                                <  2:503974/247‑1 (MQ=255)
                                    gggTCGACCGACGAATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCACGGACCGACGCCGCGTAATGTGACCGATGTTCTCCCAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTAGGAGCTGGGGTCGATCACTCCGAAAAACGGTGCGGTTATGGCGCTGGTAGGGGGGTTTGAtttaactcaaagca                                                                                                                                                                                                                >  2:261104/1‑239 (MQ=255)
                                      gTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGt                                                                                                                                                                                                             <  1:351632/251‑1 (MQ=255)
                                        cGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGttt                                                                                                                                                                                                           >  1:14059/1‑251 (MQ=255)
                                            cGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAACCGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCCCAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGGGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCCATCAGAGCAAGTTTAAcc                                                                                                                                                                                                       >  2:560039/1‑251 (MQ=255)
                                             gTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACcg                                                                                                                                                                                                      <  1:441219/251‑1 (MQ=255)
                                               cgcatagaatctggagagCGTGCGCTTGGCGCGTCCTGACCGTTCGGATACTCAGCAAGGACCGAAGCCGCGTACAGGGACCGATGTTCTGCAAACGGGTCAGCAAGTCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCGGGCACAAGTGCCGGAATTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGTTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACcgc                                                                                                                                                                                                     <  2:102196/233‑1 (MQ=255)
                                               cGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACcgcg                                                                                                                                                                                                    >  1:126560/1‑251 (MQ=255)
                                                 cATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGcc                                                                                                                                                                                                  >  1:446718/1‑251 (MQ=255)
                                                  aTTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGcc                                                                                                                                                                                                  >  1:532329/1‑250 (MQ=255)
                                                     gAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATATTCAGCAAGGACCGACGTCGCGTAAAGTGACCGATGGTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCATAGCCAGTTTATCCGCGCCAccc                                                                                                                                                                                              <  2:68622/251‑1 (MQ=255)
                                                             aaGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGc                                                                                                                                                                                                                                   <  2:418905/206‑1 (MQ=255)
                                                             aaGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGc                                                                                                                                                                                                                                   >  1:418905/1‑206 (MQ=255)
                                                                 cggtCGCTGGGCGCGTCCTTACCGTTCGGATACCCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAATCGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTc                                                                                                                                                                                  <  2:354614/248‑1 (MQ=255)
                                                                  tatCGCTGGGCGAGTCCTTACAGTTCGGATACTCAGCAAGGACCGACGACGCGTAAAGTGACCGCTGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGGGATGCATGGGGGCTGGCGCAAGTGCCCGAAGTGAACGCGGCGCTGGTGTCGATCAATCCGCAAAAAAGGGCTGTTATGGCGCTGGTAGGTGGCTTTGATTTCAAACCGCTCAAGTTTAACCACGCCACCCAGGCatttcgtc                                                                                                                                                                                  >  2:189582/3‑243 (MQ=255)
                                                                        tGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCAGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTggg                                                                                                                                                                           >  1:298183/1‑251 (MQ=255)
                                                                         gggCGCGGCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTggg                                                                                                                                                                           <  1:8657/250‑1 (MQ=255)
                                                                                 cattACCGTTCGGATACTCAGCAAGGCCCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGGAGTGAACTCGGCGCTGGTGTCGATCAATCCGCCAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCCCCACCCAGGCACTTAGTAAGGTGGAGTCCAACAt                                                                                                                                                                  >  2:529249/3‑251 (MQ=255)
                                                                                    tACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCa                                                                                                                                                                <  1:104062/250‑1 (MQ=255)
                                                                                     aCCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCaa                                                                                                                                                               >  1:411156/1‑250 (MQ=255)
                                                                                       cGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCACAcc                                                                                                                                                            >  1:247110/1‑251 (MQ=255)
                                                                                              aTACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACACGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGGGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCctc                                                                                                                                                     >  2:459344/1‑251 (MQ=255)
                                                                                                     gCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACcg                                                                                                                                              >  1:237329/1‑251 (MQ=255)
                                                                                                       aaGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACcgcg                                                                                                                                            >  1:409926/1‑251 (MQ=255)
                                                                                                       aaGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACcgcg                                                                                                                                            <  1:116381/251‑1 (MQ=255)
                                                                                                           aacGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAAcc                                                                                                                                                                                                       >  2:447467/3‑188 (MQ=255)
                                                                                                           aCCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAAcc                                                                                                                                                                                                       <  1:447467/188‑1 (MQ=255)
                                                                                                             cgacgaccagTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGCGTTCGTCCGGTTGGCGATGCATGGTGGCTGGCACAAGAGCCGGAAGTGAACTCGGAGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCTCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAACCCGTTCCCCTACACCGCGGCGATg                                                                                                                                      >  2:288200/10‑251 (MQ=255)
                                                                                                              gACGCCGCGTAAAGGGACCGATGTTCGGCAAACGGGTCAGCAAATCTGGTGTCGGCAGGTTGGGGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATgg                                                                                                                                     <  2:446718/251‑1 (MQ=255)
                                                                                                                cgccgcGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGAAGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGAt                                                                                                                                   <  1:476153/251‑1 (MQ=255)
                                                                                                                 tccgcGTAAAGTGACCGATGTTCTGCCAACGGGGCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTCGCACACCTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGAGCTGTTATGGCGCTGGTCGGCGGCTTTGATTTCAATCAGAGAAAGTTTCACCGCGCCACCCAGGCCCTGCGTCCGGTGGGTTCCACCATCCACCCGTTCCTCAACACagcgacgcaggat                                                                                                                                   >  2:228755/2‑241 (MQ=255)
                                                                                                                 gccgcGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATa                                                                                                                                  <  1:106075/251‑1 (MQ=255)
                                                                                                                   cgcgTAAAGTGACAGATGTTCTGCAAACGGGTCAGCAAATATGGGTTCGTCAGGTTGGCGATGCATGGCGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGGCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGTTGGCTTTGATTTAACTCGGAGCAGGTTTAACCACGCCACCCAGGACCTGCGTACGGTGGGTTCCAACATCAAACAGTTCCTCTACACCGCGGCGATggcaca                                                                                                                                 >  2:535485/1‑246 (MQ=255)
                                                                                                                    gcgTAAAGGGACAGATGTTCGGCAAACGGGTCAGCAAAGCTGGGTTCGTCAGGTGGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGACGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAg                                                                                                                               <  1:510576/251‑1 (MQ=255)
                                                                                                                      gTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGATTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGCACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGCTTTCAATCAGAGCAAGTTTAACCGCGCCCCCCAGGCACTGAGTCAGGTGGGTTCAAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGt                                                                                                                             >  2:341249/1‑251 (MQ=255)
                                                                                                                      gTAAAGTGACCGACGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCAAAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCCCCCAGGCACTGCGTACGGTGGGTTCCAACATCAACCAGTTCCTCTCCACCGCGGCGATGAATAAAGGt                                                                                                                             >  2:272189/1‑251 (MQ=255)
                                                                                                                                      tCTGCAAACGGGTAAGCGAATTCGGGTTCGTCAGTTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATTGCGCGGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCTTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGt                                                                                                              <  1:41268/250‑1 (MQ=255)
                                                                                                                                              cGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGCGGCTGGCACAAGGGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAAAGCAAGTTTAACCGCGCCACCCCGGCACTGCGTCAGGTGGTTACCAACATCAACCCGTTCCTCTACACCGCGGCGATGGGTAAAGGCCTGACGCTGGCACGTATGTTgcac                                                                                                     >  2:442785/1‑248 (MQ=255)
                                                                                                                                               gggTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGGTAAAGGCCTGCCCCTGGCAAGTATgctgca                                                                                                      >  2:176362/1‑244 (MQ=255)
                                                                                                                                                       aaaTCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCa                                                                                            <  1:403470/251‑1 (MQ=255)
                                                                                                                                                             gggTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTtc                                                                                       <  1:331035/250‑1 (MQ=255)
                                                                                                                                                                  cGTCGGGTTGGCTATGCATGGTGGGTGGCAAAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGGGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACGGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCt                                                                                  <  2:250775/250‑1 (MQ=255)
                                                                                                                                                                  cGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCt                                                                                  <  1:459344/250‑1 (MQ=255)
                                                                                                                                                                  cGTCAGGGTGGCGTTGCAGGGTGGCTGGCACAAGTGCCGGTAGTGAACTCGGCGCGGGTGTCGATCAATCCGCAAAACGGGTCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTg                                                                                 <  2:350123/251‑1 (MQ=255)
                                                                                                                                                                   gTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTgg                                                                                >  1:132118/1‑251 (MQ=255)
                                                                                                                                                                     cAGGTTGGCGATGCATGGTGTCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGGGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGAGTTCAATCAGAGCAAGTTTAACCGCGCCACACAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTggg                                                                               <  2:444793/250‑1 (MQ=255)
                                                                                                                                                                     cAGGTTGGCGATGCATGGTGGCTGTCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGa                                                                              <  1:160461/251‑1 (MQ=255)
                                                                                                                                                                      gggTTGGCGATGCATGCGTGCTGGCACAAGTGCCGGAAGTGAACGCGTCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATGTCAATCAGAGCAAGTTTAACCGCGCCTCCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAATGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATGTCTCGCTGGGAt                                                                             <  2:52513/250‑1 (MQ=255)
                                                                                                                                                                       gtttGGCGATGCATGGGGGCTGGCACAAGGGCCGGAAGTGAACTCGTCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCATAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAATCGTTCCTCTACACCGCGGCGATGGATAAAGGGCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATg                                                                            <  2:334347/249‑1 (MQ=255)
                                                                                                                                                                         gtggtgatggttggtggCTGGCACAAGGGCCGGAAGTGAACTCGGCGCTTGTTTCGATCAATCCGCAAAACGGTTCTGTTATGGCGCGGGTCGGTGGCTTTGATTTCAAGCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTAAGGTGGTTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTTAACGATGTGCCAATTTATCTCTGGGATGCa                                                                          <  2:316936/240‑1 (MQ=255)
                                                                                                                                                                          tGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCa                                                                          >  1:431216/1‑250 (MQ=255)
                                                                                                                                                                            gtgATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACGCGGCGCTGGGGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATGTCAATCAGAGCAAGCTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGt                                                                       <  2:128736/249‑1 (MQ=255)
                                                                                                                                                                            gagATGCGTGGCGGCTGGCACTAGTGCTGGTAGTGACCTCGGCGCTGGTGGTGAGTAATCCGCAAAGCGTGACTGTTATGGTGCTGGTTGGTTGCTTTGATTCTACACAGAGCAAGTTTAACCGCGCCACCCATGCTCTGCGTCAGGTGGGTTCCAACACCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTTTCGCTGGGATGCAAGt                                                                       <  2:414551/249‑1 (MQ=255)
                                                                                                                                                                                     ggtgggtgGCACCAGTGCCGACAGTGAACTCGGCGCTGGTGTCGTTCAAGCAGCAAAACGGTGCTGTTAGGGCGCTGGTCGGTAGCTGTGATTTCAATCAGAGTAAGTTTAACCGCGCCGCCCAGGCACTGCGTCAGGTGGGTTCCAACATCAATCCGTTCCGCTACACCGCGGCGATGGATAAAGGTCTTACTCGGGCAAGTAGGTTGAACGGTGGGCCAATTTCTCGCTGGGATGCAAGTGCCGGGTCt                                                              <  2:181005/251‑1 (MQ=255)
                                                                                                                                                                                       gtggCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGGTTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTg                                                             <  2:242585/251‑1 (MQ=255)
                                                                                                                                                                                             ggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTgg                                                        <  1:305924/250‑1 (MQ=255)
                                                                                                                                                                                             ggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGc                                                       <  1:212145/251‑1 (MQ=255)
                                                                                                                                                                                                     tGCCGGCAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGATTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCga                                                 >  2:513332/1‑249 (MQ=255)
                                                                                                                                                                                                       ccGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCgaagaa                                             <  1:197490/251‑1 (MQ=255)
                                                                                                                                                                                                        cGGAAGTGAATTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTTATTTCAATCAGAGCAAGTTTAACCGCGCCCCCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGAGGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAAc                                            <  2:434825/251‑1 (MQ=255)
                                                                                                                                                                                                               gAACTCGTCTCTGGGGTCTATCAAGCCGCAAAACGGTGCTGTTATGGCGCTGGGCGGTGGCGTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTcaccac                                     <  2:473396/251‑1 (MQ=255)
                                                                                                                                                                                                                        gcTGGTGTCGATCACTCCGCATAACGGGGCTGTGATGGCGCTGGTCGGGGGCTTTGATTTCAACCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTAGGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATg                            <  2:444652/251‑1 (MQ=255)
                                                                                                                                                                                                                            gtgtCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTgg                        <  1:238798/251‑1 (MQ=255)
                                                                                                                                                                                                                             tgtCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCg                                                                                    >  2:467240/1‑190 (MQ=255)
                                                                                                                                                                                                                             tgtCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTggg                       >  1:393283/1‑250 (MQ=255)
                                                                                                                                                                                                                             tgtCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTAGGTTGAACGATGTGCCAATTTCTCg                                                                                    <  1:467240/190‑1 (MQ=255)
                                                                                                                                                                                                                              gccggtCAATCCGCAAAACGGGGCTGTTATGGCGCTGGTCGGTGGCTTTGATTCCAAACCGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTc                      <  2:126560/246‑1 (MQ=255)
                                                                                                                                                                                                                                  aTCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAAt                  >  1:173809/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aTCAATCCGCAAAACGGTGCTGTTATGGCGCTGGGCGGTGGCGTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAAt                  <  1:173578/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aaattcGCAAAACTGTGCTGTTATGGCGATGGTAGGGGGCGTTGATTTCAATCTGATCAAGGATAACCGCGCCACCCAGGCACTGCGGCAGGTGGGTTCAAACATCCAACCGTTCCTATACACCGCGGCGAGGGATAAAGGTCTGACGCTGGCAATTATGTTGATCGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCCGCCGAAGAACTCACCACCGCAGTATGCTGGTCAAATTc                <  2:532329/246‑1 (MQ=255)
                                                                                                                                                                                                                                          gCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATGTCAATCATAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCGGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTAc          <  2:173809/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAgggg  <  1:498212/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 < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: