Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,831,764 C→G 7.8% A436A (GCC→GCG yicH → putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,831,7640CG 7.8% 201.7 / 8.4 66A436A (GCC→GCGyicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base G (3/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

CAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGG  >  NC_000913/3831464‑3832056
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                    
cAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccgcc                                                                                                                                                                                                                                                                                                      >  2:546414/1‑301 (MQ=255)
      ttttccgatcTGGCTTCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCa                                                                                                                                                                                                                                                                                                <  1:75537/292‑1 (MQ=255)
         aGGGAACTGTCTGCGTGACGGGAAAACGTTGATTCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATGTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTc                                                                                                                                                                                                                                                                                             <  2:539086/301‑1 (MQ=255)
            ggacggggtGCGTGACGGGAAAACGTTGATTCTGGAGGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAAt                                                                                                                                                                                                                                                                                          <  1:524715/293‑1 (MQ=255)
              aCTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAATGGCACGGCGAATCTGAATGCCGCCCCAGCGACCTTCaaccc                                                                                                                                                                                                                                                                                        >  1:566539/1‑298 (MQ=255)
                 ggCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCaagatcggaaga                                                                                                                                                                                                                                                                                     >  2:75627/1‑290 (MQ=255)
                     ccGTGATGGAAAAACTTTTTTCGTTTATGATGCGGAAATTGCCGGGCTTGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGAt                                                                                                                                                                                                                                                                                 <  1:89044/300‑1 (MQ=255)
                            gggAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtcgtc                                                                                                                                                                                                                                                                          <  1:545741/301‑1 (MQ=255)
                                     gtggtCCGTGATGGGGGGGGAATTGCCGGGTTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                                                                                                                                                                                                                                 <  1:527022/297‑1 (MQ=255)
                                         ttcgttgtgatgCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATGGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTGACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgc                                                                                                                                                                                                                                                             <  1:629842/294‑1 (MQ=255)
                                                     cGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGCTTTAGCGCCAGCCGCAATCTGCTCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTCCCGATCATAAATGGGGCGTCTGGAATGGCTCGGCGCATCTGAATGCCCCCGCCGCGACCTTCACACCACTTGATGTTCGTCCCCCACCCCTCGCCCAGACCGCCac                                                                                                                                                                                                                                                 >  1:397610/1‑300 (MQ=255)
                                                      ggCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCACCCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCCCGGCGAATCTGAATGCCGCCGCCGCGACATTCCAGCGTGGTGATGGTCGTCGCCCGTcccccgccgcgccggcacac                                                                                                                                                                                                                                                >  1:14572/1‑282 (MQ=255)
                                                                  gCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAACCTGACTGCCGCCGCCCCGACATTCAATCGTGTTGATGTTCGTCGCCCGTAGCTGGCGCTGACCGCCAACAACAGCACggag                                                                                                                                                                                                                                    >  1:223192/1‑299 (MQ=255)
                                                                        atatACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAataaa                                                                                                                                                                                                                              >  2:28396/1‑299 (MQ=255)
                                                                             ccTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGa                                                                                                                                                                                                                         >  2:475216/1‑301 (MQ=255)
                                                                                gcttttccgttctgCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGAGTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACt                                                                                                                                                                                                                      <  1:364396/288‑1 (MQ=255)
                                                                                  ccGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa                                                                                                                                                                                                                    >  2:312646/1‑301 (MQ=255)
                                                                                   cGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg                                                                                                                                                                                                                   <  2:505973/301‑1 (MQ=255)
                                                                                     aaaaaCTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg                                                                                                                                                                                                                 <  1:474615/301‑1 (MQ=255)
                                                                                          cTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAttt                                                                                                                                                                                                            <  2:41155/301‑1 (MQ=255)
                                                                                            ggCAACAGTTGTGGATGGAAACGACACCCTGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCTCCATCTGATCATTTATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc                                                                                                                                                                                                          >  2:16308/1‑301 (MQ=255)
                                                                                             gCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACtagataggaagaga                                                                                                                                                                                                         >  2:364880/1‑288 (MQ=255)
                                                                                                     tgtgGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGCGAATAACAGCGAACTGAGTGCATTTACTGaaaaaaa                                                                                                                                                                                                 >  2:48294/1‑299 (MQ=255)
                                                                                                       tggttggaAACGACTCCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCa                                                                                                                                                                                               <  2:484588/301‑1 (MQ=255)
                                                                                                                                   aCAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGATTcaca                                                                                                                                                                   >  2:47911/1‑301 (MQ=255)
                                                                                                                                        cTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGAGTTTCACAAAcaa                                                                                                                                                              >  2:481952/1‑299 (MQ=255)
                                                                                                                                         gtcAAGTGAGGTGATTTAGCGCCAGCCGCAGGTTGATCATTGATATCGACCGTGGCTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGcc                                                                                                                                                             <  1:390496/299‑1 (MQ=255)
                                                                                                                                                gAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCCGACCGCCACCAGCAGCACGGTGAATATCCGCGAACTGAGTGCATTTACTGAAAAAGGCCTTCTGGAAACCACCGCCAGTGTTTCCCAACCACCACAACgc                                                                                                                                                      >  1:492908/1‑300 (MQ=255)
                                                                                                                                                       tttAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACACCACAACGTCAGacac                                                                                                                                               >  1:144864/1‑301 (MQ=255)
                                                                                                                                                        ttAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCAAAACGTCAGAcaca                                                                                                                                              >  2:528506/1‑301 (MQ=255)
                                                                                                                                                            cgcCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATc                                                                                                                                          <  2:144990/301‑1 (MQ=255)
                                                                                                                                                                       ttCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATgg                                                                                                                               <  2:92723/300‑1 (MQ=255)
                                                                                                                                                                               tttGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGgtg                                                                                                                       <  2:567253/300‑1 (MQ=255)
                                                                                                                                                                                   ttatCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGcc                                                                                                                   <  2:223435/300‑1 (MQ=255)
                                                                                                                                                                                                    ccGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGAGCCGGTGAATATTTTGCaaa                                                                                                  >  2:410597/1‑300 (MQ=255)
                                                                                                                                                                                                      ttGGCAGCTCACCGCTCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAg                                                                                                <  2:472834/300‑1 (MQ=255)
                                                                                                                                                                                                       tGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGt                                                                                               <  2:493554/301‑1 (MQ=255)
                                                                                                                                                                                                          cAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTggg                                                                                            >  2:333020/1‑301 (MQ=255)
                                                                                                                                                                                                             cTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCCCAAACGCCACAACGTCAGACACATAACAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACCATAGGGAt                                                                                         >  1:223818/1‑301 (MQ=255)
                                                                                                                                                                                                             cTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGAt                                                                                         <  1:28345/301‑1 (MQ=255)
                                                                                                                                                                                                              ctaCCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATg                                                                                        <  1:312242/299‑1 (MQ=255)
                                                                                                                                                                                                                             gTTGCGGTGCCTACCTGATGTTGTTTACCGATCATAAATGGGGCGTCTGGAGGGGCTCGGCGAAGCGGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATAGCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGtt                                                                         <  1:16257/301‑1 (MQ=255)
                                                                                                                                                                                                                              ttACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGTACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTg                                                                        >  2:429354/1‑301 (MQ=255)
                                                                                                                                                                                                                                 tcttGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACt                                                                     <  2:170464/298‑1 (MQ=255)
                                                                                                                                                                                                                                  ggtgcagtggtgtcGCTTGTTACTGATCATAAGTGGGGCGTCTGGGGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCGGGCGCTGACCGCCAACAGCAGCACGGGGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTg                                                                    <  1:138339/288‑1 (MQ=255)
                                                                                                                                                                                                                                    tGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGGGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTAACtaaa                                                                  >  1:170267/1‑298 (MQ=255)
                                                                                                                                                                                                                                       cAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCACTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGAAGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTCGACTGGAGAc                                                               >  2:420930/1‑301 (MQ=255)
                                                                                                                                                                                                                                            ggcggCTGGTTACGGGTGATAAATGGGGGGTGTGGAGTGGCTCGGCGAAGCTGAAGTCCGCCGCCGCGACATTCAAGCGTGTTGAGGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCaa                                                          <  1:374665/297‑1 (MQ=255)
                                                                                                                                                                                                                                                 ctggtggccgggcATAAAGGGGGCGTCGGTGGGTGCTCGGCGAATCTGATTGCCGCCGCCGCGACATTCATTCGTGTTGATTGTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGATTGCATTTACTTAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTc                                                     <  1:47858/289‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCa                                                    <  2:14572/301‑1 (MQ=255)
                                                                                                                                                                                                                                                    gTTATCGGTGATAAATGGGGCGTCTGGAGTGGCTTGGTGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAACTACCGTTGACTGGCGACGGCAATATTCAGc                                                  <  1:152766/301‑1 (MQ=255)
                                                                                                                                                                                                                                                          gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCCACGTCAATATTCAGCTTcccc                                            >  1:381519/1‑300 (MQ=255)
                                                                                                                                                                                                                                                            tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATTGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGATTACAGca                                          >  2:115210/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                   gggggCGTTTGGAGTGGCTCGGCGAATCTGAATGCCGGCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCg                                   <  1:410034/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccgctcttccgttctaagttGATTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAgg                          <  1:420153/282‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cgtgtgctcttccgatcttGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGc                         <  2:160374/283‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCCGCTACCCGCCAGTGGCAATAATCAGGCCAAtt                    >  1:581467/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                    gttggggagtCTGAATGCCGCCGGTGCGGTATTCTATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGACAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTc                  <  1:527833/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGAAGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGACCAATGTcc                 >  2:371671/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          gattcttagtGCCGCCTCCGCTACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTTAGTGCATGTACTGAAGAAGGCATGCTGGAAGCCACCGCCAGTGTTTCACAAACGACACAACGTCAGACACAAGTCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTg            <  1:420367/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          gAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTg            <  2:398124/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAggagatcggaagagcaa           >  2:420716/1‑286 (MQ=255)
                                                                                                                                                                                                                                                                                            aTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGCTGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAAGCCAATGTCCCGTTGaa          >  1:471131/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               tGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGAGGCGAAATTCAGGcagatcggaagagcgacgg       >  1:160223/1‑283 (MQ=255)
                                                                                                                                                                                                                                                                                               tGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAAcc       >  2:528268/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               tGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTCAAAcc       >  1:247090/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                gAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCACCAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAAAGTCAGACCCATATCAGCCTGAATGGACGCTGTGTGCCGGTGAATATTTTGCAACAGTGGGGATCGCCATAATTCACTTTGACTGGCCACAGCCATATTCCGCTTACCACCAGAGACCATAGTCAGACCCATCTCCCGttacacacc      >  1:352357/1‑294 (MQ=255)
                                                                                                                                                                                                                                                                                                    gccgccgccGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGCCCAATGCCCCGTTCAAACCTACgc  >  1:598830/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                    
CAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTAAACAGCCTGCAACTGAAGAGATTTAGCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGG  >  NC_000913/3831464‑3832056

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