Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,435,815 A→T 11.2% L50Q (CTG→CAG)  PP_1256 ← putative alpha‑ketoglutarate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,435,8150AT11.2% 257.3 / 20.3 80L50Q (CTG→CAG) PP_1256putative alpha‑ketoglutarate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base T (5/4);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

GATCCGCACGTCCAGCCATTCGCCGGCCCGCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGATTACGACCCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGA  >  AE015451/1435583‑1436046
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       
gATCCGCACGTCCAGCCATTCGCCGGCCCGCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgc                                                                                                                                                                                                                                       <  1:68568/235‑1 (MQ=255)
     gCACGTCCAGCCATTCGCCGGCCCGCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAg                                                                                                                                                                                                                                 <  1:1280460/236‑1 (MQ=255)
        cGTCCAGCCATTCGCCGGCCCGCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgcg                                                                                                                                                                                                                                      <  1:951878/228‑1 (MQ=255)
         gTCCAGCCATTCGCCGGCCCGCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTgg                                                                                                                                                                                                                                >  2:1365764/1‑233 (MQ=255)
             aGCCATTCGCCGGCCCGCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTgg                                                                                                                                                                                                                           <  2:520017/234‑1 (MQ=255)
              gCCATTCGCCGGCCCGCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGc                                                                                                                                                                                                                         <  1:988820/235‑1 (MQ=255)
                 aTTCGCCGGCCCGCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTgg                                                                                                                                                                                                                           <  1:78124/230‑1 (MQ=255)
                     gCCGGCCCGCACTGTACGGGCAAAGGGGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgc                                                                                                                                                                                                                                       <  2:400642/214‑1 (MQ=255)
                      ccGGCCCGCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTc                                                                                                                                                                                                                >  1:276879/1‑236 (MQ=255)
                             gCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTgg                                                                                                                                                                                                                           <  1:1341867/218‑1 (MQ=255)
                                cTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTgc                                                                                                                                                                                                        <  1:1365767/234‑1 (MQ=255)
                                   tACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCa                                                                                                                                                                                                                                  <  1:1366853/205‑1 (MQ=255)
                                      gggCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGa                                                                                                                                                                                                  >  2:1309868/1‑234 (MQ=255)
                                      gggCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAg                                                                                                                                                                                                 >  1:527317/1‑235 (MQ=255)
                                                gCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgc                                                                                                                                                                                                                                       <  2:51628/187‑1 (MQ=255)
                                                gCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCt                                                                                                                                                                                      >  1:1233099/1‑236 (MQ=255)
                                                 cGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTc                                                                                                                                                                                     <  1:1032800/236‑1 (MQ=255)
                                                     gCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCa                                                                                                                                                                                                                                         >  2:13810/1‑180 (MQ=255)
                                                          gTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTg                                                                                                                                                                            <  1:1072222/236‑1 (MQ=255)
                                                              aGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgc                                                                                                                                                                                                                                       <  2:227410/173‑1 (MQ=255)
                                                                    ggCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgc                                                                                                                                                                                                                                       <  2:385745/167‑1 (MQ=255)
                                                                        aCGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTggg                                                                                                                                                               >  1:605557/1‑235 (MQ=255)
                                                                        aCGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTggg                                                                                                                                                               >  1:1109241/1‑235 (MQ=255)
                                                                          gCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTc                                                                                                                                                             <  1:1186759/235‑1 (MQ=255)
                                                                              gCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCggtg                                                                                                                                                                     <  2:605555/223‑1 (MQ=255)
                                                                                       gcgcgcCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGc                                                                                                                                                                                       <  2:574418/196‑1 (MQ=255)
                                                                                       gcgcgcCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGc                                                                                                                                                                                       >  1:574420/1‑196 (MQ=255)
                                                                                           gcCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGCTg                                                                                                                                                >  1:164334/1‑231 (MQ=255)
                                                                                           gcCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTc                                                                                                                                           >  1:1248100/1‑236 (MQ=255)
                                                                                              gtGGGCAGGCCGGTTTCGATCACGGCGCGGTCGAGCAGGGCATCGCCGAGGGCCTCGAGCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTgc                                                                                                                                                                                                        <  2:172422/171‑1 (MQ=255)
                                                                                                         ggTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGt                                                                                                                               >  2:497221/1‑234 (MQ=255)
                                                                                                          gTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTAc                                                                                                                            >  1:761549/1‑236 (MQ=255)
                                                                                                          gTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGAAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTAc                                                                                                                            >  1:1231598/1‑236 (MQ=255)
                                                                                                              cGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGc                                                                                                                                              >  2:1383473/1‑214 (MQ=255)
                                                                                                               gATCACGGCGCGGTCGATCAGGGCATCGCCGAGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCtgtg                                                                                                                                                                                                         <  1:1481078/153‑1 (MQ=255)
                                                                                                                  cACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGAAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCctc                                                                                                                     <  2:1231595/235‑1 (MQ=255)
                                                                                                                   aCGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGGATGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTg                                                                                                                   >  1:932422/1‑236 (MQ=255)
                                                                                                                     ggCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGt                                                                                                                               <  1:13810/222‑1 (MQ=255)
                                                                                                                              gATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagca                                                                                                        <  1:140671/236‑1 (MQ=255)
                                                                                                                              cttgaagccATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagca                                                                                                        >  1:402910/9‑236 (MQ=255)
                                                                                                                               aTCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagca                                                                                                        <  1:1309871/235‑1 (MQ=255)
                                                                                                                                        tCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGAt                                                                                                             <  1:1178210/221‑1 (MQ=255)
                                                                                                                                         cGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGCGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgc                                                                                                                                                                                                                                       <  2:438974/98‑1 (MQ=255)
                                                                                                                                         cGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGGTACCGGCCTCTGGCCGCTCAGCAGGTGGCCtgt                                                                                              >  1:252986/1‑235 (MQ=255)
                                                                                                                                           ccGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGAcgg                                                                                                                                                                        >  1:344502/1‑157 (MQ=255)
                                                                                                                                                    tCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGt                                                                                                                               <  2:223880/191‑1 (MQ=255)
                                                                                                                                                     cGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTcgc                                                                                                                                         <  1:502077/180‑1 (MQ=255)
                                                                                                                                                      gATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGc                                                                                >  1:1456013/1‑236 (MQ=255)
                                                                                                                                                             gtgGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGcc                                                                                                                >  2:369248/1‑197 (MQ=255)
                                                                                                                                                               ggCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGt                                                                                   >  2:1358282/1‑224 (MQ=255)
                                                                                                                                                               ggCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGt                                                                                   <  1:1358285/224‑1 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATTTTTCGCGGTAGGCATCCAATGCCGCCCAGGCCTgcgcgct                                                                                                                                                                                                                                  >  1:524396/1‑70 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgc                                                                                                                                                                                                                                       >  1:572025/1‑66 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgc                                                                                                                                                                                                                                       <  2:572023/66‑1 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgc                                                                                                                                                                                                                                       <  2:592166/66‑1 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgc                                                                                                                                                                                                                                       >  1:592168/1‑66 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGcc                                                                                                                >  1:171235/1‑185 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCtt                                                                          >  1:198074/1‑223 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAAt                                                               >  2:1218433/1‑234 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATg                                                              >  1:1284075/1‑235 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGc                                                             >  1:994815/1‑236 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGGCGCGGCTGCCGGGTACCGGCCGCTGGcc                                                                                                                >  1:165902/1‑185 (MQ=255)
                                                                                                                                                                         tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGGCGCGGCGGCCGGTTCCCGGCCTCTGGCCGAt                                                                                                             >  1:3267/1‑188 (MQ=255)
                                                                                                                                                                                  tctGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATg                                                    <  1:1218436/236‑1 (MQ=255)
                                                                                                                                                                                      cTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCg                                                                                                               <  2:767983/173‑1 (MQ=255)
                                                                                                                                                                                      cTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCg                                                                                                               >  1:767985/1‑173 (MQ=255)
                                                                                                                                                                                       tCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTg                                                                                                   >  1:630069/1‑184 (MQ=255)
                                                                                                                                                                                           tgttCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGAtt                                           <  1:526509/236‑1 (MQ=255)
                                                                                                                                                                                                 agtgATGTGTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGATGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATggg                                                  <  1:526317/220‑1 (MQ=255)
                                                                                                                                                                                                     aTGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGATTAc                                         <  2:949797/228‑1 (MQ=255)
                                                                                                                                                                                                            gcgGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCtgtg                                                                                             <  2:248939/169‑1 (MQ=255)
                                                                                                                                                                                                             cgGTAGGCATCGAAGGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGAtt                                           <  2:1009789/218‑1 (MQ=255)
                                                                                                                                                                                                                tAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGATTACGACCCCTGTAGGAGCGGGt                      <  1:369250/236‑1 (MQ=255)
                                                                                                                                                                                                                   gCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGATTACGACCCCTGTAGGAGCGGGtt                     <  1:103316/234‑1 (MQ=255)
                                                                                                                                                                                                                           cgtggccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATgggg                                                 >  1:647501/5‑198 (MQ=255)
                                                                                                                                                                                                                             ccgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGATTACGa                                       >  2:1342154/1‑206 (MQ=255)
                                                                                                                                                                                                                                cccAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGATTACGACCCCTGTAGGAGCGGGTTCACCCGCGAACACCg      >  1:926360/1‑236 (MQ=255)
                                                                                                                                                                                                                                 ccAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGGTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCgcg                                                                                                             >  2:373944/1‑130 (MQ=255)
                                                                                                                                                                                                                                   aGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTTGCCTCCTTGGAAATCTAATGCAGAt                                                         >  1:64509/1‑182 (MQ=255)
                                                                                                                                                                                                                                     gCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGATTACGACCCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGa  <  1:1383476/235‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGATTACGACCCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGa  <  2:527315/233‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       
GATCCGCACGTCCAGCCATTCGCCGGCCCGCACTGTACGGGCAAAGGTGCGTAGCTGGGTGCAGGTACGGCCACGCTCGCCCTGAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCTCCTTGGAAATCTAATGCAGATGCATGGGGGGGATTACGACCCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGA  >  AE015451/1435583‑1436046

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