Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,349,365 C→A 5.6% D20Y (GAC→TAC)  PP_1174 ← conserved protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,349,3650CA 5.6% 293.4 / 2.6 85D20Y (GAC→TAC) PP_1174conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base C (44/36);  new base A (2/3);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTCTTCCTTTTCCGCAACGCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTAGACCCGGCCGCAACCCTGCGCTATTTTCAGTTCTGGCCTCGCCCTC  >  AE015451/1349137‑1349591
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  
ctctTCCTTTTCCGCAACGCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGc                                                                                                                                                                                                                               >  1:651987/1‑234 (MQ=255)
           ccGCAACGCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCgcag                                                                                                                                                                                                                     <  2:383862/233‑1 (MQ=255)
             gCAACGCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGtt                                                                                                                                                                                                                <  2:560808/236‑1 (MQ=255)
             gCAACGCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGtt                                                                                                                                                                                                                <  1:440326/236‑1 (MQ=255)
             gCAACGCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGtt                                                                                                                                                                                                                >  1:164645/1‑236 (MQ=255)
             gCAACGCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGtt                                                                                                                                                                                                                <  1:889304/236‑1 (MQ=255)
              cAACGCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTTTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTg                                                                                                                                                                                                               >  1:225130/1‑236 (MQ=255)
                  gCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCAcc                                                                                                                                                                                                           <  1:641819/236‑1 (MQ=255)
                  gCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCAcc                                                                                                                                                                                                           <  1:143315/236‑1 (MQ=255)
                    ttCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGtt                                                                                                                                                                                                                <  1:256417/229‑1 (MQ=255)
                              aGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCAcc                                                                                                                                                                                                           <  1:228664/224‑1 (MQ=255)
                               gATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGtt                                                                                                                                                                                                                <  2:962087/218‑1 (MQ=255)
                                gtCGAGCGCTCCCTGGCGCCTCGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCgcag                                                                                                                                                                                                                     <  2:164643/211‑1 (MQ=255)
                                aTCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGtt                                                                                                                                                                                                                <  1:781907/217‑1 (MQ=255)
                                              ggCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTc                                                                                                                                                                                 <  2:651984/234‑1 (MQ=255)
                                                  ccACGGCCGCCAGCTAGCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGc                                                                                                                                                                                                                               <  2:155552/184‑1 (MQ=255)
                                                         cgccAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTg                                                                                                                                                                                                <  1:815347/208‑1 (MQ=255)
                                                              gCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACATTCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACt                                                                                                                                                                                                 <  1:648826/202‑1 (MQ=255)
                                                                    ttGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCAt                                                                                                                                                                           >  2:422028/1‑218 (MQ=255)
                                                                     tGGCCACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAg                                                                                                                                                                                                                          <  2:175488/170‑1 (MQ=255)
                                                                       gCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATg                                                                                                                                                      >  1:664532/1‑236 (MQ=255)
                                                                           ctccCCCGTGCTACGCGCTGCCCTGAATGGCCTGGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGtt                                                                                                                                                                                                                <  2:155290/174‑1 (MQ=255)
                                                                              ccccGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGtt                                                                                                                                                                                                                <  2:224487/170‑1 (MQ=255)
                                                                                ccGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATcc                                                                                                                                                                                            <  2:731778/189‑1 (MQ=255)
                                                                                 cGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAgc                                                                                                                                            <  1:992684/236‑1 (MQ=255)
                                                                                 cGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCACGCAGTGTtcgtc                                                                                                                                                               >  1:292892/1‑217 (MQ=255)
                                                                                      tACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGtt                                                                                                                                                                                                                <  2:241683/163‑1 (MQ=255)
                                                                                           gcTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGt                                                                                                                                   >  1:731282/1‑235 (MQ=255)
                                                                                                ccTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGc                                                                                                                                                                         >  1:836158/1‑192 (MQ=255)
                                                                                                  tGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTACGGTACTGATCCGGGTCGAGCGCGGGCATGCAGGGTTCGTCGAAGAAATGACCGCTCAgccggc                                                                                                                                        >  1:32951/1‑223 (MQ=255)
                                                                                                        gCCTTGTGCCGCGCAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTGGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCga                                                                                                                                                             <  1:825736/196‑1 (MQ=255)
                                                                                                            tgtgCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTg                                                                                                                              >  2:1070174/1‑223 (MQ=255)
                                                                                                            tgtgCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTACTTGCCATAGCGCAGTAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTc                                                                                                                                                                                 >  1:771059/1‑172 (MQ=255)
                                                                                                             gtgCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCgaagaa                                                                                                                                                         >  2:715169/1‑195 (MQ=255)
                                                                                                                 cgcgAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCgaa                                                                                                                                                            >  2:430929/1‑188 (MQ=255)
                                                                                                                 cgcgAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAAc                                                                                                                    >  1:851297/1‑228 (MQ=255)
                                                                                                                 cgcgAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCt                                                                                                              >  2:591520/1‑234 (MQ=255)
                                                                                                                 cgcgAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTg                                                                                                            >  1:764388/1‑236 (MQ=255)
                                                                                                                 cgcgAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCACCAGCGGCACCGTGCGGTa                                                                                                                                                                                                   >  2:204435/1‑149 (MQ=255)
                                                                                                                        gggCCGCACAGCCGCCCCGGCAAACAATCCATCGGTGCAGATCCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTATTTGCCATAGCGCACTAGCGGCACAGTGCGGTACTGATCCGCGTCGa                                                                                                                                                                                     <  2:771056/156‑1 (MQ=18)
                                                                                                                                      ccccGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAg                                                                                                                                    >  1:225891/1‑191 (MQ=255)
                                                                                                                                      ccccGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTATTTGCCATAg                                                                                                                                                                                                                          <  1:481306/105‑1 (MQ=255)
                                                                                                                                      ccccGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTATTTGCCATAg                                                                                                                                                                                                                          >  2:481303/1‑105 (MQ=255)
                                                                                                                                         cGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGa                                                                                                    >  1:514652/1‑220 (MQ=255)
                                                                                                                                          ggCAAACAATCCATCAGTGCAGAGCCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCagag                                                                                             <  1:265773/226‑1 (MQ=255)
                                                                                                                                                    ccATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAgccggccg                                                                                                                                      >  2:895772/1‑175 (MQ=255)
                                                                                                                                                     catcatTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCg                                                                                                                                                  >  2:364405/1‑162 (MQ=255)
                                                                                                                                                     catcatTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCg                                                                                                                                                  <  1:364408/162‑1 (MQ=255)
                                                                                                                                                      atcatTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCg                                                                                                                                                                                >  2:181661/1‑131 (MQ=255)
                                                                                                                                                      atcatTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGCGCGCAGAACCAGGCTTGGGATTCGAGcc                                                                                                 >  1:1032050/1‑210 (MQ=255)
                                                                                                                                                           tGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAAt                                                                  <  1:422031/236‑1 (MQ=255)
                                                                                                                                                                aaCCCCCCAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCAGGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCt                                                                                                                                                <  2:932345/153‑1 (MQ=255)
                                                                                                                                                                    cccAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTcaatcc                                                           <  2:664529/234‑1 (MQ=255)
                                                                                                                                                                     ccAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGc                                                                                                                                                 >  2:380850/1‑147 (MQ=255)
                                                                                                                                                                     ccAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAAc                                                                                                                    >  2:349393/1‑176 (MQ=255)
                                                                                                                                                                        aGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATCCGAg                                                                                                   >  2:179112/1‑190 (MQ=255)
                                                                                                                                                                              gTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTa                                                <  1:1070177/235‑1 (MQ=255)
                                                                                                                                                                               tGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCa                                                                                                                        >  2:397391/1‑162 (MQ=255)
                                                                                                                                                                               tGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCa                                                                                                                        <  1:397394/162‑1 (MQ=255)
                                                                                                                                                                                gCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTAGa                                              <  1:81782/235‑1 (MQ=255)
                                                                                                                                                                                       aCAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGt                                                                          >  2:499959/1‑200 (MQ=255)
                                                                                                                                                                                            gTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTAGACCCGGCCGCAAcc                                 >  1:745075/1‑236 (MQ=255)
                                                                                                                                                                                              gtaGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAg                                                                                                                       <  2:117465/147‑1 (MQ=255)
                                                                                                                                                                                              ctaGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGagag                                                                                  <  2:815599/184‑1 (MQ=255)
                                                                                                                                                                                              ataGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGtt                                                                         >  2:81780/1‑194 (MQ=255)
                                                                                                                                                                                                   ccccGGACTCGCTGACCATGGTGGTCTCCTGGTCTTTGCCATAg                                                                                                                                                                                                                          >  1:897885/1‑44 (MQ=255)
                                                                                                                                                                                                   ccccGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTAGACCCGGc                                        >  1:1074227/1‑222 (MQ=255)
                                                                                                                                                                                                   ccccGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTAGACCCGGCCGCAACCCTgcg                            >  2:938134/1‑234 (MQ=255)
                                                                                                                                                                                                        gACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAAtt                                                                 >  1:440227/1‑192 (MQ=255)
                                                                                                                                                                                                         aCTCGCTGTCCATGGCGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTcaatcca                                                          <  2:880751/198‑1 (MQ=255)
                                                                                                                                                                                                         aCTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCTGGGt                                                                                                                                                                                        >  1:555006/1‑72 (MQ=255)
                                                                                                                                                                                                          cTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGtt                                                                                                                                  >  2:772283/1‑125 (MQ=255)
                                                                                                                                                                                                          cTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTAGACCCGGCCGCAACCCTg                              >  2:964191/1‑225 (MQ=255)
                                                                                                                                                                                                          cTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTAGACCCGGCCGCAACCCTGCGCTATTTTCa                   >  1:665526/1‑236 (MQ=255)
                                                                                                                                                                                                          cTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGCTTACAAATTCAATCCAATCGCCTGt                                                 >  1:102848/1‑206 (MQ=255)
                                                                                                                                                                                                          cTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGATGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGATATGAc                                                                                                                                                    >  2:1086456/1‑107 (MQ=255)
                                                                                                                                                                                                          cTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCACCAGCTGCACCGTGCGGTa                                                                                                                                                                                                   >  2:55992/1‑60 (MQ=17)
                                                                                                                                                                                                            cGCTGGCCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCGCCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTAGACCCgg                                         <  1:127897/212‑1 (MQ=255)
                                                                                                                                                                                                               tGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAgg                                                                                            >  2:806931/1‑158 (MQ=255)
                                                                                                                                                                                                               tGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGgtgt                                                                                        <  1:981304/162‑1 (MQ=255)
                                                                                                                                                                                                                gACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGCTCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGCGGGCGTGGAGAGGACATg                                                                            >  2:59687/1‑173 (MQ=255)
                                                                                                                                                                                                                           tctcCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTAGACCCGGCCGCAACCCTGCGCTATTTTCAGTTCTGGCCTCGCCCTc  >  1:571395/1‑236 (MQ=255)
                                                                                                                                                                                                                            ctcCTGTTACTTGCCATAGCGCAGCAGTTGAAACGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGTTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCag                                                                                               <  1:51212/142‑1 (MQ=17)
                                                                                                                                                                                                                                   tCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGGGCCAGATAGTGTGGAggggcggttg                                                                           >  2:794158/1‑149 (MQ=14)
                                                                                                                                                                                                                                   tCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGCGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGGGCCAGATAGTGTGGAggggcggttg                                                                           <  1:794161/155‑7 (MQ=11)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  
CTCTTCCTTTTCCGCAACGCTTCCCGGCTAAGATCGAGCGCTCCCTGGCGCCACGGCCGCCAGCTATCTTGGCAACTCTCCCGTGCTACGCGCTTCCCTGAATGGCCTTGTGCCGCGAAAGGGCCGCACAGCCGCCCCGGCAAACAATCCATCATTGCAGAACCCCCAAGATCAGTGCTTTGTACAAGGTATAGGCCCCGGACTCGCTGACCATGGTGGTCTCCTGGTCCTTGCCATAGCGCAGCAGTTGCACCGTGCGGTACTGATCCGGGTCGAGCTCGGGCATGCAGTGTTCGTCGAAGAAATGACCGCTCAGCCGGCCGAGTTCGTGGGCGCAGAACCAGGCTTGGGATTCGAGCCAGAGGGTGTGGAGAGGACATGGTTTACAAATTCAATCCAATCGCCTGTAGACCCGGCCGCAACCCTGCGCTATTTTCAGTTCTGGCCTCGCCCTC  >  AE015451/1349137‑1349591

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