Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,414,820 T→C 12.6% L6L (CTT→CTC PP_2115 → conserved protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,414,8200TC12.6% 124.8 / 6.7 47L6L (CTT→CTCPP_2115conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base C (2/4);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGCGATCATTTTCGGAGCCCCCCATGACCGATTTCACCCCCGACGCCAT  >  AE015451/2414599‑2415047
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
cATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTCATCCACGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCg                                                                                                                                                                                                                        >  1:501746/1‑235 (MQ=255)
                  ggcatcggcggggcattgagctaaaatcaaggcaagatcccgatgctctcgctctcgctctcgctctcgctctcgctctcgcttcggcttttgattttgcttctaagcgcgcgatagcccaggcgccgcggattgcgacttcaggaggccgagcgcagggcttgcgaagggaggtgacgggcatccATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGcctc                                                                                                                                                                                                               <  2:434872‑M1/40‑1 (MQ=255)
                       cggcggggcattgagctaaaatcaaggcaagatcccgatgctctcgctctcgctctcgctctcgctctcgctctcgcttcggcttttgattttgcttctaagcgcgcgatagcccaggcgccgcggattgcgacttcaggaggccgagcgcagggcttgcgaagggaggtgacgggcatccATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaa                                                                                                                                                                                                  <  2:749465‑M1/53‑1 (MQ=255)
                            tttCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCCTCCATGCCCGTCACCTGCCTTCg                                                                                                                                                                                                                                   >  2:638704/1‑196 (MQ=18)
                               cGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaacatc                                                                                                                                                                                        >  1:473421/1‑236 (MQ=14)
                                                    tttGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaag                                                                                                                                                                   <  1:537819/236‑1 (MQ=31)
                                                          cgatgctctcgctctcgctctcgctctcgctctcgctctcgcttcggcttttgattttgcttctaagcgcgcgatagcccaggcgccgcggattgcgacttcaggaggccgagcgcagggcttgcgaagggaggtgacgggcatggATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTc                                                                                                                                                                                                       >  2:743018‑M1/147‑194 (MQ=255)
                                                          aaaCGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaa                                                                                                                                                                                                  <  2:476419/199‑1 (MQ=32)
                                                                                      cTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagctcaagatcaagatcaagatcaag                                                                                                                                     >  1:234115/1‑232 (MQ=21)
                                                                                        gcgcTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGTGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaag                                                                                                                                                       <  1:42044/212‑1 (MQ=14)
                                                                                          gcTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgag                                                                                                                               >  2:789194/1‑234 (MQ=37)
                                                                                              ccATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgag                                                                                                                               <  2:678549/230‑1 (MQ=38)
                                                                                                         cttttgattttgcttctaagcgcgcgatagcccaggcgccgcggattgcgacttcaggaggccgagcgcagggcttgcgaagggaggggacgggcatggATGCCCGTCACCTCCCTTCGCAAGCCCTGcgc                                                                                                                                                                                                                       >  2:867478‑M1/100‑130 (MQ=255)
                                                                                                                tgttgcttctaagcgcgcgatagttcaggcgacgcagattgcgacttcaggaggccgagcgtagggcttgcgaagggaggtgacgggcatggATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaag                                                                                                                         >  1:197625‑M1/93‑218 (MQ=255)
                                                                                                                      ttgtaagcgcgcggtagttcaggcgacgcagattgcgacttcagaaggccgagcgcagggcttgcgaagggaggtgacgggcatggATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatc                                                                                                                                                                                  <  2:383571‑M1/69‑1 (MQ=255)
                                                                                                                      ttctaagcgcgcgatagcccaggcgccgcggattgcgacttcaggaggccgagcgcagggcttgcgaagggaggtgacgggcctggATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagat                                                                                                                                                                             <  2:649037‑M1/74‑1 (MQ=255)
                                                                                                                                cccTCACTCCGGCCTTGCTCCCGGCAGGCCCGCGCCGAAGGCCCCCTCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaag                                                                                                                                                                                                 <  2:199114/130‑1 (MQ=17)
                                                                                                                                       tCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaag                                                                                                                                                       <  2:561202/165‑1 (MQ=31)
                                                                                                                                       tCCGGCCGTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaa                                                                                        <  1:93189/228‑1 (MQ=255)
                                                                                                                                                         cgacttcaggaggccgagcgtagggcttgcgaagggaggtgacgggcatggATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCt                                                               >  1:183839‑M1/52‑235 (MQ=255)
                                                                                                                                                         aGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCt                                                               <  1:817446/235‑1 (MQ=255)
                                                                                                                                                          gacttcaggaggtcgcgcgcagggcgtgcggagggaggtgacgggcatggACGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaag                                                                                                                                                                   >  1:431932‑M1/51‑134 (MQ=255)
                                                                                                                                                                 gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagaacaagat                                                                                                                                   >  2:25820/1‑159 (MQ=25)
                                                                                                                                                                   cGAAGGCCCTATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGt                                                                                             >  1:185292/1‑195 (MQ=255)
                                                                                                                                                                   cGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGatcaagatcaagatcgagatcaagatcaagagcgag                                                                                                             >  1:668533/1‑179 (MQ=255)
                                                                                                                                                                    gcaGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACcg                                                                        <  2:440258/213‑1 (MQ=255)
                                                                                                                                                                    gcaGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACcg                                                                        >  1:440261/3‑215 (MQ=255)
                                                                                                                                                                        gCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGc                                                <  1:317125/235‑1 (MQ=255)
                                                                                                                                                                                cTGCGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaa                                                                                        <  2:495294/187‑1 (MQ=255)
                                                                                                                                                                                   ttgcgaaggcaggtgacgggcatccATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaa                                                                                        <  1:645191‑M1/159‑1 (MQ=255)
                                                                                                                                                                                   ttgcgaaggcaggtgacgggcatccATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaa                                                                                        >  2:645188‑M1/26‑184 (MQ=255)
                                                                                                                                                                                       cTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaa                                                                                        <  2:227218/180‑1 (MQ=255)
                                                                                                                                                                                       cTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaa                                                                                        >  1:227220/1‑180 (MQ=255)
                                                                                                                                                                                       cTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGAcc                                                                         <  1:632623/195‑1 (MQ=255)
                                                                                                                                                                                        tCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATa                                                  <  1:115394/217‑1 (MQ=255)
                                                                                                                                                                                        tCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGCGATCATTTTCGGa                                   <  1:640768/232‑1 (MQ=255)
                                                                                                                                                                                           gcgcTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCTAGATCAAGACCAAGATCGAGATCAAGAGCGCTggg                                                                                              >  2:550425/1‑170 (MQ=255)
                                                                                                                                                                                             gcTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatcaagatccagat                                                                                                                                   >  2:844459/1‑131 (MQ=37)
                                                                                                                                                                                             gcTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCG‑‑‑‑‑‑AAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGCGATCATTTTCGGAGCCCCCCATGACCGa                    >  1:917185/1‑236 (MQ=255)
                                                                                                                                                                                               ttGGCGGGCATCCAGGCCAGTCCCCGGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaagatcaag                                                                                                                   <  2:173052/145‑1 (MQ=17)
                                                                                                                                                                                                        aTCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGcgcg                                                        >  2:849404/1‑195 (MQ=255)
                                                                                                                                                                                                            ATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTGCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGc                                                >  1:441279‑M1/1‑199 (MQ=255)
                                                                                                                                                                                                                  GTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGCGATCAttt                                        <  1:475805‑M1/201‑1 (MQ=255)
                                                                                                                                                                                                                    CACCTGCCTTCGCAAGGCCTGCGTTCGGCCTACTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaa                                                                                        >  2:50335‑M1/1‑151 (MQ=255)
                                                                                                                                                                                                                      ccTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAg                                                 >  2:385149/1‑188 (MQ=255)
                                                                                                                                                                                                                       cTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGCGATCATTTTCGGAGCCCCCCATGACCGATTTCACCCCCGACGCCAt  >  1:301870/1‑234 (MQ=255)
                                                                                                                                                                                                                        tGCCTTCGCAAGGCCTGCGTTCAGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaagatcaagatcgagatcaagatcaagatcgagat                                                                                                           <  1:224263/128‑1 (MQ=255)
                                                                                                                                                                                                                            ttCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGCGATCATTTTCGGAGCCCCCCATGACCGATTTCACCCCCg         <  1:912331/222‑1 (MQ=255)
                                                                                                                                                                                                                            ttCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGCGATCATTTTCGGAGCCCCCCATGACCGATTTCACCCCCg         >  2:912328/1‑222 (MQ=255)
                                                                                                                                                                                                                             tCGCAAGGCCTGCGTTCGGCATCCTTAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGCGATCATTTTCGGAg                                  <  2:69309/196‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
CATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCGCTATGCGCGACAATAGCGATCATTTTCGGAGCCCCCCATGACCGATTTCACCCCCGACGCCAT  >  AE015451/2414599‑2415047

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: