Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,725,652 C→G 7.2% P438P (CCG→CCC)  opgH ← Glucans biosynthesis glucosyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154515,725,6520CG 7.2% 242.1 / 8.9 69P438P (CCG→CCC) opgHGlucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base C (23/41);  new base G (2/3);  total (25/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

ACCGCACGGTGAACCGGGTGCATGCCCTTGACCAGGAACAGGCGGAAGTTCATCAGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTCGGCCCGGTCTGGATGATGCCGGCGTCCGGGTTGGCCTCCATCAGGC  >  AE015451/5725422‑5725878
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
acCGCACGGTGAACCGGGTGCATGCCCTTGACCAGGAACAGGCGGAAGTTCATCAGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCg                                                                                                                                                                                                                                   >  2:747779/1‑232 (MQ=255)
   gCACGGTGAACCGGGTGCATGCCCTTGACCAGGAACAGGCGGAAGTTCATCAGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCg                                                                                                                                                                                                                                   <  2:669777/229‑1 (MQ=255)
         tGAACCGGGTGCATGCCCTTGACCAGGAACAGGCGGAAGTTCATCAGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAg                                                                                                                                                                                                                      >  1:346068/1‑236 (MQ=255)
                    cATGCCCTTGACCAGGAACAGGCGGAAGTTCATCAGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCt                                                                                                                                                                                                           >  1:531276/1‑236 (MQ=255)
                    cATGCCCTTGACCAGGAACAGGCGGAAGTTCATCAGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGc                                                                                                                                                                                                                                 >  2:376771/1‑214 (MQ=255)
                               ccAGGAACAGGCGGAAGTTCATCAGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGa                                                                                                                                                                                                     <  2:323803/231‑1 (MQ=255)
                                        ggCGGAAGTTCATCAGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGc                                                                                                                                                                                         <  2:392089/234‑1 (MQ=255)
                                         gCGGAAGTTCATCAGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGAt                                                                                                                                                                                                       <  1:60984/219‑1 (MQ=255)
                                                  catcaGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGt                                                                                                                                                                             <  1:36728/236‑1 (MQ=255)
                                                          ggCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATcgc                                                                                                                                                                               <  1:681179/225‑1 (MQ=255)
                                                                                gcgcTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCtt                                                                                                                                                                                              >  2:309534/1‑189 (MQ=255)
                                                                                   cTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCa                                                                                                                                                                                                                                >  1:232300/1‑152 (MQ=255)
                                                                                      gAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAgg                                                                                                                                                                                                                     <  1:88397/160‑1 (MQ=255)
                                                                                      gAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTtgt                                                                                                                                                                          <  2:336016/203‑1 (MQ=255)
                                                                                            gTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCa                                                                                                                                   <  1:747780/236‑1 (MQ=255)
                                                                                               aaGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGAt                                                                                                                                                                                                       <  2:538242/165‑1 (MQ=255)
                                                                                               aaGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGaa                                                                                                                                        <  1:577683/228‑1 (MQ=255)
                                                                                                     gTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACg                                                                                                                                                                                                  >  1:114786/1‑164 (MQ=255)
                                                                                                     gTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACg                                                                                                                                                                                                  <  2:114786/164‑1 (MQ=255)
                                                                                                        tGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGACGTAGGCAATCCACTCGCCCCAGCCGGCACGGCACATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCAccc                                                                                                                                                   <  1:206374/208‑1 (MQ=255)
                                                                                                                gTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAgg                                                                                                                                                                                                                     <  1:433491/134‑1 (MQ=255)
                                                                                                                   cttcGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGcc                                                                                                             <  1:336325/235‑1 (MQ=255)
                                                                                                                     tcGTGGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAg                                                                                                                                                                                                               <  2:19860/135‑1 (MQ=255)
                                                                                                                      cGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCg                                                                                                                              <  2:643062/215‑1 (MQ=255)
                                                                                                                            tGCCCGGCAAGTCGTAGGCAAGCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTaca                                                                                                        <  1:193610/231‑1 (MQ=255)
                                                                                                                                cGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTgg                                                                                                                                                                        <  2:635866/163‑1 (MQ=255)
                                                                                                                                 ggCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCAccc                                                                                                                                                   >  1:736790/1‑183 (MQ=255)
                                                                                                                                     aGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGaa                                                                                           <  2:586167/235‑1 (MQ=255)
                                                                                                                                      gTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGcc                                                                                                             <  1:376772/216‑1 (MQ=255)
                                                                                                                                          tAGGCAATCCACACACCCCGGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAActgct                                                                                      <  1:340291/235‑1 (MQ=255)
                                                                                                                                             gCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCAt                                                                                  <  1:202920/236‑1 (MQ=255)
                                                                                                                                              cAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCAt                                                                                  <  1:313967/235‑1 (MQ=255)
                                                                                                                                                    acacacCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAAGGCGCCTTTGCCCGGCAACGGCGCCAgg                                                                                                                                                                                                                     >  2:601906/1‑98 (MQ=255)
                                                                                                                                                    acacacCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCa                                                                                                                                   <  2:359155/180‑1 (MQ=255)
                                                                                                                                                         ccccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATAc                                                                        <  1:281397/234‑1 (MQ=255)
                                                                                                                                                          cccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAgg                                                                     >  1:118381/1‑236 (MQ=255)
                                                                                                                                                             aGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGCTGAACGGCTTCATGCGGATGATcgcg                                                                                                                                                                              >  2:190452/1‑128 (MQ=255)
                                                                                                                                                               ccGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGgtg                                                                  >  1:430751/1‑234 (MQ=255)
                                                                                                                                                                cGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGt                                                                                                                                                                             >  2:227894/1‑126 (MQ=255)
                                                                                                                                                                cGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTc                                                                >  1:677030/1‑235 (MQ=255)
                                                                                                                                                                    aCGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCg                                                                                                                              <  2:116825/169‑1 (MQ=255)
                                                                                                                                                                       gcgcATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCa                                                                                                                                   <  1:365922/161‑1 (MQ=255)
                                                                                                                                                                       gcgcATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCa                                                                                                                                   >  2:365921/1‑161 (MQ=255)
                                                                                                                                                                       gcgcATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACa                                                                                                                    >  1:111091/1‑176 (MQ=255)
                                                                                                                                                                       gcgcATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGccc                                                         >  1:581012/1‑235 (MQ=255)
                                                                                                                                                                                cgcAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTCg                                               <  1:692656/236‑1 (MQ=255)
                                                                                                                                                                                cgcAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTCg                                               <  1:176371/236‑1 (MQ=255)
                                                                                                                                                                                   aGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGGGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACtgctgc                                                                                    <  2:724249/196‑1 (MQ=255)
                                                                                                                                                                                    gCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGc                                                                                                                            <  1:37293/155‑1 (MQ=255)
                                                                                                                                                                                     cTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTggg                                                                                                               >  2:188478/1‑167 (MQ=255)
                                                                                                                                                                                     cTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAActgct                                                                                      >  1:101341/1‑192 (MQ=255)
                                                                                                                                                                                       tCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAActgct                                                                                      <  2:560727/190‑1 (MQ=255)
                                                                                                                                                                                       tCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTCg                                               <  2:677029/229‑1 (MQ=255)
                                                                                                                                                                                         gacggAGGCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTg                                                                                                                                              <  2:54378/132‑1 (MQ=255)
                                                                                                                                                                                             aaGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGc                                                                            <  1:649727/194‑1 (MQ=255)
                                                                                                                                                                                                  gTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGACGAACGGCTTCATGCGGATGATCGCGTGGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAActgct                                                                                      <  2:663637/179‑1 (MQ=255)
                                                                                                                                                                                                           gATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTCGGCCCGGTCTGGATGATGCCgg                          <  1:337719/230‑1 (MQ=255)
                                                                                                                                                                                                           gATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTCGGCCCGGTCTGGATGATGCCGGCGTCCg                    >  1:351458/1‑236 (MQ=255)
                                                                                                                                                                                                             tGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCagt                                                                                                                                                                 <  1:580264/93‑1 (MQ=255)
                                                                                                                                                                                                             tGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCagt                                                                                                                                                                 >  2:580263/1‑93 (MQ=255)
                                                                                                                                                                                                               gCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGccc                                                         >  1:368889/1‑195 (MQ=255)
                                                                                                                                                                                                                 aCCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACa                                                                       >  1:382047/1‑179 (MQ=255)
                                                                                                                                                                                                                  ccGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTCGGCCCGGTCTGGATGATGCCGGCGTCCgg                   <  2:531275/230‑1 (MQ=255)
                                                                                                                                                                                                                       cctggcgcCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGa                                                                                                                       <  1:258279/121‑1 (MQ=255)
                                                                                                                                                                                                                       cctggcgcCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGa                                                                                                                       >  2:258278/5‑125 (MQ=255)
                                                                                                                                                                                                                        aaCGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGc                                                                            <  1:12260/167‑1 (MQ=255)
                                                                                                                                                                                                                        aaCGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGc                                                                            >  2:12260/1‑167 (MQ=255)
                                                                                                                                                                                                                            cgcCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTCGGCCCGGTCTGGATGATGCCGGCGTCCGGGTTGGCCTCCATCAg    <  1:336517/235‑1 (MQ=255)
                                                                                                                                                                                                                                   gCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTc                                                <  1:422393/184‑1 (MQ=255)
                                                                                                                                                                                                                                    cccGGCAACGGCGCCAGGGCGCAGGGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTCGGCCCGGTCTGGATGATGCCGGCGTCCGGGTTGGCCTCCATCAGgc  <  2:515842/229‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
ACCGCACGGTGAACCGGGTGCATGCCCTTGACCAGGAACAGGCGGAAGTTCATCAGGTTGCCGTGGCACCAGCGGCGGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCATACAGGGTGTCCATGCCCGAGGCTTTCGGCCCGGTCTGGATGATGCCGGCGTCCGGGTTGGCCTCCATCAGGC  >  AE015451/5725422‑5725878

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑