Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,348,705 A→G 57.1% intergenic (‑70/+404) PP_3820 ← / ← galU Group II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,348,7050AG57.1% 5.5 / 82.6 63intergenic (‑70/+404)PP_3820/galUGroup II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (8/19);  new base G (15/21);  total (23/40)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

GCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCAC  >  AE015451/4348502‑4348930
                                                                                                                                                                                                           |                                                                                                                                                                                                                                 
gCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                                                                                                   <  1:680517/236‑1 (MQ=11)
  tGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTgg                                                                                                                                                                                                 <  1:55190/236‑1 (MQ=31)
       ccgccgGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGc                                                                                                                                                                                            <  1:378717/236‑1 (MQ=32)
                    tGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTa                                                                                                                                                                                         <  1:679717/226‑1 (MQ=33)
                     gCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAATTTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCGGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGGGGCAAGCGTACCGCTCTCCCa                                                                                                                                                                              >  1:231427/1‑236 (MQ=255)
                            cAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTa                                                                                                                                                                         >  2:460280/1‑234 (MQ=18)
                            cAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTCGTGGCGCGGTAAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCa                                                                                                                                                                                               >  1:710317/1‑212 (MQ=14)
                              ggCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTc                                                                                                                                                                     >  1:201915/1‑236 (MQ=17)
                                 aTCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcg                                                                                                                                                                    <  2:272317/234‑1 (MQ=33)
                                      ccTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAcc                                                                                                                                                             <  1:347202/236‑1 (MQ=33)
                                           gCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGCCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGGGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTgg                                                                                                                                                                                                 <  1:211632/195‑1 (MQ=14)
                                              tCTGGCCGGGATTCGGGATGACAGTTCCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACCCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                                                                                                        <  1:166384/201‑1 (MQ=17)
                                                         ttCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCg                                                                                                                                          <  1:744580/236‑1 (MQ=18)
                                                          tCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCTGCCGCTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCg                                                                                                                                          <  1:404601/235‑1 (MQ=12)
                                                          tCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAcc                                                                                                                                                             <  1:346681/216‑1 (MQ=17)
                                                          tCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCg                                                                                                                                          <  1:375567/235‑1 (MQ=33)
                                                            gggATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCgg                                                                                                                                         <  2:344147/234‑1 (MQ=33)
                                                                   cAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCt                                                                                                                                     <  1:62126/231‑1 (MQ=32)
                                                                         ccGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATATCCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCg                                                                                                                          >  1:419597/1‑236 (MQ=18)
                                                                           gCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGtt                                                                                                                        >  1:674413/1‑236 (MQ=33)
                                                                              gtgGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCg                                                                                                                     <  1:460281/236‑1 (MQ=18)
                                                                              gtgGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCg                                                                                                                     <  1:748968/236‑1 (MQ=18)
                                                                                 ggAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGcc                                                                                                                                               >  1:502219/1‑207 (MQ=33)
                                                                                 ggAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGcc                                                                                                                                               <  2:502218/207‑1 (MQ=33)
                                                                                 ggAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCt                                                                                                                                     <  1:546188/217‑1 (MQ=33)
                                                                                   aaaTTTTTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTc                                                                                                                                                                     <  1:344043/183‑1 (MQ=18)
                                                                                   aaaTTTCTGCATTACGGGTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCaa                                                                                                                <  1:323898/236‑1 (MQ=17)
                                                                                               tACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCtt                                                                                                                                                    <  2:357467/188‑1 (MQ=33)
                                                                                                         aaCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtg                                                                                           >  1:643594/1‑235 (MQ=33)
                                                                                                                       gCCTCCTGGCGGCATACCCTGTCCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCg                                                                                                                                          <  1:128998/174‑1 (MQ=16)
                                                                                                                             gggcggcAAAACCTGTTCGGCCCTTGGGGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGcc                                                                                                                                               <  1:90169/162‑1 (MQ=16)
                                                                                                                               gcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTg                                                                                                   >  1:385205/1‑205 (MQ=17)
                                                                                                                               gcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCAcc                                                                                    <  1:673027/220‑1 (MQ=18)
                                                                                                                                cggcATAACCTGTGCGGCCCTGGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgc                                                                                                                                                                   <  1:99553/140‑1 (MQ=255)
                                                                                                                                      aaCCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCtg                                                                         <  1:17774/224‑1 (MQ=17)
                                                                                                                                       cccTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCtg                                                                         <  1:197068/222‑1 (MQ=17)
                                                                                                                                              cGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACa                                                                                                                                                                                       >  2:676425/1‑106 (MQ=35)
                                                                                                                                               ggCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACCGCGTTCCGTGCAAGTAtt                                                                                                           >  2:497807/1‑181 (MQ=21)
                                                                                                                                               ggCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTTGCGAtt                                                                                             >  1:51428/1‑195 (MQ=18)
                                                                                                                                               ggCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTc                                                      >  2:354504/1‑234 (MQ=18)
                                                                                                                                                   cTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGACAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGt                                                  >  2:605098/1‑234 (MQ=18)
                                                                                                                                                        tggCGCGGTTACCCGCCACTTCCTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                                                                                                        <  1:84538/95‑1 (MQ=255)
                                                                                                                                                            gcgGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGACACCTCACCCCGCTGTGCCGCCACTATCCGCTTCCTGtt                                                 <  2:419596/226‑1 (MQ=18)
                                                                                                                                                                  aaCCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                                                                                                                                         <  2:71779/116‑1 (MQ=25)
                                                                                                                                                                   aCCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGAttt                                                                                            <  1:156748/176‑1 (MQ=18)
                                                                                                                                                                   aCCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCAtt                                  >  2:387091/1‑234 (MQ=18)
                                                                                                                                                                     cgccgcTTCCTACGGCTTCGGCGGACGTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAAGTCGCGCGACCGTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtg                                                                                           <  2:206188/175‑1 (MQ=11)
                                                                                                                                                                     cGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGc                              <  1:354505/236‑1 (MQ=18)
                                                                                                                                                                        cACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTc                           <  1:525988/236‑1 (MQ=18)
                                                                                                                                                                             actacGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGt                                              >  2:476125/1‑212 (MQ=17)
                                                                                                                                                                             actacGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGct                      >  1:728376/1‑236 (MQ=18)
                                                                                                                                                                             actacGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGCTCTCCCAGGGTAATTCGCGCGACCTTCCTg                                                                                                                                                       >  2:61732/1‑107 (MQ=25)
                                                                                                                                                                               tacGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGcttc                    <  1:72023/236‑1 (MQ=17)
                                                                                                                                                                                     ttCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAg              >  1:486649/1‑236 (MQ=18)
                                                                                                                                                                                      tCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCGTCAGa             >  1:162509/1‑236 (MQ=17)
                                                                                                                                                                                       cGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCCAGTATTGGGCTTTGGCGATTTGTGCCACCTCCCCCCGCtg                                                                         >  2:67974/1‑175 (MQ=18)
                                                                                                                                                                                        ggCTGACTTCTGCACGCTCGTCCCGTGGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                                                                                                   <  1:732632/52‑1 (MQ=255)
                                                                                                                                                                                          cTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCg                                                                                                                     >  2:285310/1‑128 (MQ=25)
                                                                                                                                                                                          cTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCg                                                                                                                     <  1:285311/128‑1 (MQ=37)
                                                                                                                                                                                                tCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATcc                                                           >  2:732156/1‑180 (MQ=21)
                                                                                                                                                                                                tCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATcc                                                           <  1:732157/180‑1 (MQ=21)
                                                                                                                                                                                                  tGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCAc  <  1:522296/235‑1 (MQ=32)
                                                                                                                                                                                                     aCGCTCGTCCCCTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCg                                                                                                                                          <  2:155669/96‑1 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                 
GCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCAC  >  AE015451/4348502‑4348930

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