Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,980,087 T→G 5.7% E66D (GAA→GAC PP_2605 ← conserved protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,980,0870TG 5.7% 233.3 / 3.0 66E66D (GAA→GACPP_2605conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base T (41/21);  new base G (2/2);  total (43/23)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

GCCCTAGCCGTCGCCCGTCGTAAGGTAATGCACACAGCGAGCGACCGCCGCGTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTCACATTCAAGGGGGGGCATGGAAGGGGTCCGCTCGGATGTTGGAAGAAGTGGC  >  AE015451/2979855‑2980318
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       
gCCCTAGCCGTCGCCCGTCGTAAGGTAATGCACACAGCGAGCGACCGCCGCGTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTcg                                                                                                                                                                                                                                      <  1:1178876/236‑1 (MQ=255)
   cTAGCCGTCGCCCGTCGTAAGGTAATGCACACAGCGAGCGACCGCCGCGTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTcgac                                                                                                                                                                                                                                    >  1:426309/1‑235 (MQ=255)
        cGTCGCCCGTCGTAAGGTAATGCACACAGCGAGCGACCGCCGCGTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATg                                                                                                                                                                                                                                >  2:795852/1‑234 (MQ=255)
        cGTCGCCCGTCGTAAGGTAATGCACACAGCGAGCGACCGCCGCGTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGc                                                                                                                                                                                                                               <  1:422488/235‑1 (MQ=255)
                               acacaGCGAGCGACCGCCGCGTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGACCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGc                                                                                                                                                                                                                        >  1:35521/1‑219 (MQ=255)
                                          gaCCGCCGCGTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCcac                                                                                                                                                                                            >  1:801374/1‑236 (MQ=255)
                                             cgccgcGTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCcaca                                                                                                                                                                                           <  1:283794/234‑1 (MQ=255)
                                               ccgcGTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTgc                                                                                                                                                                                                       <  2:1140841/220‑1 (MQ=255)
                                                cgcgTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAaccac                                                                                                                                                                                      >  1:1113826/1‑236 (MQ=255)
                                                        cTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTcg                                                                                                                                                                                                                                      <  2:819637/180‑1 (MQ=255)
                                                           gatggcagtgcACGCATGTCAGGCATACGGCTACCGCACTGGAAGTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGAAGCCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGGATCca                                                                                                                                                                                             >  1:665612/11‑218 (MQ=255)
                                                           ccATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCa                                                                                                                                                                           <  2:976201/236‑1 (MQ=255)
                                                             aTCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCa                                                                                                                                                                           >  1:910356/1‑234 (MQ=255)
                                                                        gCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTgc                                                                                                                                                                                                       >  1:496520/1‑195 (MQ=255)
                                                                        gCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTgcc                                                                                                                                                               >  1:996015/1‑235 (MQ=255)
                                                                               gggCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGc                                                                                                                                                                                                                        >  2:663809/1‑171 (MQ=255)
                                                                               gggCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGc                                                                                                                                                                                                                        <  1:663812/171‑1 (MQ=255)
                                                                                       ggCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCGTTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGgc                                                                                                                                                                        >  1:776833/1‑211 (MQ=255)
                                                                                         cTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGCCCAGGCGCTGGTGCCGccc                                                                                                                                                           >  1:51032/1‑222 (MQ=255)
                                                                                              gCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCcg                                                                                                                                        >  1:478883/1‑236 (MQ=255)
                                                                                               cACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTAt                                                                                                                                                                                                >  2:526791/1‑179 (MQ=255)
                                                                                                aCTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCcgc                                                                                                                                       <  1:795855/235‑1 (MQ=255)
                                                                                                  tGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGcgac                                                                                                                                    <  1:483673/236‑1 (MQ=255)
                                                                                                            cGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGt                                                                                                                           >  1:1128739/1‑235 (MQ=255)
                                                                                                             ggCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAgctggctg                                                                                                                                                  >  1:1115213/1‑211 (MQ=255)
                                                                                                             ggCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCAc                                                                                                                              <  1:1032725/231‑1 (MQ=255)
                                                                                                             ggCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCAc                                                                                                                              >  2:1032722/1‑231 (MQ=255)
                                                                                                             ggCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCa                                                                                                                         >  1:429991/1‑236 (MQ=255)
                                                                                                                  tgtgCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAgct                                                                                                                                                       >  1:800622/1‑201 (MQ=255)
                                                                                                                  tgtgCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGctggctgg                                                                                                                                                 >  1:824488/1‑207 (MQ=255)
                                                                                                                  tgtgCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTggctggct                                                                                                                                               >  1:516844/1‑209 (MQ=255)
                                                                                                                  tatgctgcTCGAAGAACGCTACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCa                                                                                                                    >  1:761188/3‑236 (MQ=255)
                                                                                                                            gaagaaCGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGcgcg                                                                                                                                                                                                    >  2:72453/1‑146 (MQ=255)
                                                                                                                            gaagaaCGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATcc                                                                                                                                                                                              >  1:35079/1‑152 (MQ=255)
                                                                                                                            gaagaaCGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCtggctggcgt                                                                                                                                              >  2:811529/1‑198 (MQ=255)
                                                                                                                            gaagaaCGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGGATCCACCACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCttt                                                                                                                                             >  1:225242/1‑201 (MQ=255)
                                                                                                                                 aCGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGt                                                                                                      >  1:495423/1‑235 (MQ=255)
                                                                                                                                      accaTTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCCCCCCCCCGACCAGGCGCTGGTgc                                                                                                                                                                >  1:906073/1‑172 (MQ=255)
                                                                                                                                                 gccATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTgc                                                                                                                                                                                                       <  1:217798/121‑1 (MQ=255)
                                                                                                                                                   cATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAgccaggcc                                                                                    <  1:927172/235‑1 (MQ=255)
                                                                                                                                                        cgcAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGcc                                                                                                                     >  1:67881/1‑197 (MQ=255)
                                                                                                                                                        cgcAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAgg                                                                               >  1:864339/1‑235 (MQ=255)
                                                                                                                                                        cgcAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGAGGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGgc                                                                              >  1:714455/1‑236 (MQ=255)
                                                                                                                                                         gcAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCAGGATGCGCGGCGGtt                                                                                                     >  1:794302/1‑212 (MQ=255)
                                                                                                                                                            ggAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACggcggc                                                                                             <  2:877699/217‑1 (MQ=255)
                                                                                                                                                             gAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTgc                                                                                                                                                                                                       <  2:283882/110‑1 (MQ=255)
                                                                                                                                                             gAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCg                                                                          <  1:931823/235‑1 (MQ=255)
                                                                                                                                                                ggCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGc                                                                      >  1:511660/1‑236 (MQ=255)
                                                                                                                                                                     gTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGAACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGtt                                                                                                                                                                                                         >  2:321329/1‑100 (MQ=255)
                                                                                                                                                                              ggTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTAc                                                         <  1:939487/235‑1 (MQ=255)
                                                                                                                                                                                            gtgtGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAgcccggcccgga                                                                                >  1:175607/1‑194 (MQ=255)
                                                                                                                                                                                                gggTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTCACATTCAAggggggg                                       >  1:1110971/1‑235 (MQ=255)
                                                                                                                                                                                                 ggTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTCACATTCAAgggggg                                        <  2:761185/233‑1 (MQ=255)
                                                                                                                                                                                                  gTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTCACATTCAAGGGGGGGCAt                                    >  1:604130/1‑236 (MQ=255)
                                                                                                                                                                                                        gcacgcaGGTAGCGCTCCGGGGCATCGTCGGCGTCGACGATGCCGTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCcg                                                                                                                                        <  2:35521/130‑1 (MQ=255)
                                                                                                                                                                                                         cacgcaGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGAACGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTACGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTa                                                                                                    >  2:632046/1‑165 (MQ=255)
                                                                                                                                                                                                          acgcaGATAGCGCTCCGGGGCATCGTCGAAGCCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCtt                                                                                                                                              >  1:1033854/1‑122 (MQ=255)
                                                                                                                                                                                                          acgcaGATAGCGCTCCGGGGCATCGTCGAAGCCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCtt                                                                                                                                              <  2:1033851/122‑1 (MQ=255)
                                                                                                                                                                                                              aGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTcaca                                                   >  1:886489/1‑209 (MQ=255)
                                                                                                                                                                                                              aGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTCACATTCAAGGGGGGGCATGGAAGGGGTCCg                        >  1:642715/1‑236 (MQ=255)
                                                                                                                                                                                                                  aGCGCTCCGGGGAATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTa                                                          <  1:572832/198‑1 (MQ=255)
                                                                                                                                                                                                                           gggCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTCACATTCAAg                                             <  2:1176943/202‑1 (MQ=255)
                                                                                                                                                                                                                            ggCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTCACATTCAAggggggg                                       >  1:880566/1‑207 (MQ=255)
                                                                                                                                                                                                                                   tcgGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTCACATTCAAGGGGGGGCATGGAAGGGGTcc                         <  2:429988/214‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTCACATTCAAGGGGGGGCATGGAAGGGGTCCGCTCGGATGTTGGAAGAAGTGGc  <  1:72453/235‑1 (MQ=255)
                                                                                                                                                                                                                                      gCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTAcggcg                                                                                               >  1:124800/1‑141 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       
GCCCTAGCCGTCGCCCGTCGTAAGGTAATGCACACAGCGAGCGACCGCCGCGTGCACTGCCATCGCCTTCACGCATGTCGGGCATACGGCTACCGCACTGGAACTGCACGGCTGTGTGCTGCTCGAAGAACGCCACCATTTCCGGACCATGGCGCAGGAAGGCGCGTTGACGCGGGTCCAGCTCGCGGGTGTGGGTTTGCGCACGCAGATAGCGCTCCGGGGCATCGTCGGCTTCGACGATGCCTTCGGCGATGGCCAGCGGGTTGCGCGGTATCCACAACCACCCCCCCGACCAGGCGCTGGTGCCGCCCAGCTGGCTGGCTTTTTCCGCGACGACCACTGTCAGGCCATGATGCGCGGCGGTTACGGCGGCAGCCAGGCCAGAGGCGCCGGAGCCGATGACGAGTACGTCACATTCAAGGGGGGGCATGGAAGGGGTCCGCTCGGATGTTGGAAGAAGTGGC  >  AE015451/2979855‑2980318

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: