Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,256,267 C→T 15.1% intergenic (+327/‑377) asd → / → PP_1990 Aspartate‑semialdehyde dehydrogenase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,256,2670CT15.1% 150.8 / 20.6 62intergenic (+327/‑377)asd/PP_1990Aspartate‑semialdehyde dehydrogenase/transposase
Reads supporting (aligned to +/- strand):  ref base C (32/20);  new base T (4/5);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-02

ATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGC  >  AE015451/2256057‑2256483
                                                                                                                                                                                                                   |                                                                                                                                                                                                                         
aTGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCCCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGa                                                                                                                                                                                                             >  2:958448/1‑225 (MQ=255)
         tGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCggc                                                                                                                                                                                                             >  1:662291/1‑215 (MQ=255)
                tctGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAg                                                                                                                                                                                                       >  1:520428/1‑215 (MQ=255)
                                 cgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGt                                                                                                                                                                         >  2:1108574/1‑228 (MQ=40)
                                  gttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGgatga                                                                                                                                                                                       <  2:1033504/213‑1 (MQ=40)
                                              tttgcttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCa                                                                                                                                                      >  2:1119278/1‑234 (MQ=255)
                                                         ttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGttttt                                                                                                                                                                     >  2:964112/1‑208 (MQ=35)
                                                                          ttCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAg                                                                                                                                                                                                  <  2:553117/162‑1 (MQ=33)
                                                                                    gCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATgg                                                                                                                                                                                    >  1:484770/1‑166 (MQ=33)
                                                                                     cAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAATAAGGGACCCGCAGTa                                                                                                                               >  1:611098/1‑217 (MQ=255)
                                                                                       gATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATgc                                                                                                          >  1:122949/1‑236 (MQ=255)
                                                                                                 tCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAg                                                                                                                             >  1:440252/1‑207 (MQ=255)
                                                                                                    ggaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGCAGCGCAGGGGACGgag                                                                                                                                                                                         >  2:193153/1‑144 (MQ=21)
                                                                                                     ggggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGCCccgccg                                                                                                                                 <  2:416486/197‑1 (MQ=255)
                                                                                                     gaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtga                                                                                             <  1:729803/235‑1 (MQ=255)
                                                                                                             gcgcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCCGGGGCCGGATGATgg                                                                                                                                                                                    >  2:211131/1‑141 (MQ=21)
                                                                                                               gcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATg                                                                                                           >  2:655462/1‑211 (MQ=255)
                                                                                                               gcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGAACCGACGTAAGGGCGGAAAgg                                                                                                                   >  1:40586/1‑203 (MQ=255)
                                                                                                                     gtggCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGGCCGGAGCGAGGGGACCCCGGAGCGCAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCAc                                                                                                                                                     <  1:511695/161‑1 (MQ=255)
                                                                                                                                       gggCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTaa                                                                                                                              >  1:760139/1‑169 (MQ=255)
                                                                                                                                         gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGgatg                                                                                                                                                                                        >  1:114447/1‑109 (MQ=32)
                                                                                                                                         gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGgatg                                                                                                                                                                                        <  2:114447/109‑1 (MQ=32)
                                                                                                                                         gCACGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATCACTACTCGGGGTGa                                                                        >  2:306873/1‑220 (MQ=255)
                                                                                                                                           aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAGGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGt                                                                                                                  <  1:774473/176‑1 (MQ=255)
                                                                                                                                           aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaag                                                      >  1:486586/1‑236 (MQ=255)
                                                                                                                                             ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCAc                                                                                                                                                     >  1:336946/1‑140 (MQ=255)
                                                                                                                                             ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCAc                                                                                                                                                     <  2:336945/140‑1 (MQ=255)
                                                                                                                                             ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGgt                                                                            >  1:1069021/1‑212 (MQ=255)
                                                                                                                                               ctGCCCGTCAAGCGCTGGGCGCCAGGGTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCCAGCCCGGAGTGAGGGGAACCCGGAGTGTAGCGTAGGGGCCGGTTGATGGGAGCCGAGGGTTTTTGGTTCCTTTTTGCCAc                                                                                                                                                     <  1:386632/137‑1 (MQ=11)
                                                                                                                                                tGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAggg                                                                                                                                         <  2:1041889/148‑1 (MQ=255)
                                                                                                                                                tGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAggg                                                                                                                                         >  1:1041890/1‑148 (MQ=255)
                                                                                                                                                   ccGTCAAGCGCTGGGCCCCAGGATGGGGCCGGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtga                                                                                             <  1:1115390/189‑1 (MQ=255)
                                                                                                                                                    cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGGGCAGTCCGGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTACTTTTTGCCACGACAAAAAg                                                                                                                                           <  2:596846/142‑1 (MQ=17)
                                                                                                                                                      tCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGGCCGGAGCGAGGGGACCCCGGAGCGCAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACg                                                                                                                                                    <  2:20909/132‑1 (MQ=255)
                                                                                                                                                        aaGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGAAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAg                                                                                                                                           <  2:760138/139‑1 (MQ=14)
                                                                                                                                                        aaGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCca                                                                 >  2:940185/1‑212 (MQ=255)
                                                                                                                                                           ggcTGGGTCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCt                                                                                                   <  1:801438/174‑1 (MQ=255)
                                                                                                                                                                 ggCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCt                                                                                   <  2:980075/186‑1 (MQ=255)
                                                                                                                                                                   ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTg                                                                                                  <  2:452476/169‑1 (MQ=255)
                                                                                                                                                                    cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTc                                                                               >  1:701040/1‑187 (MQ=255)
                                                                                                                                                                    cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCca                                                           >  2:795901/1‑207 (MQ=255)
                                                                                                                                                                       aGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGt                                                                                                                  <  1:794677/149‑1 (MQ=255)
                                                                                                                                                                       aGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGt                                                                                                                  >  2:794676/1‑149 (MQ=255)
                                                                                                                                                                                     ggcgcGGTCCTCCGGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAgg                                                                                                                                                                                                 <  2:320186/56‑1 (MQ=255)
                                                                                                                                                                                     aGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGCAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCAc                                                                                                                                                     <  2:461160/101‑1 (MQ=11)
                                                                                                                                                                                     aGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGCAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCAc                                                                                                                                                     >  1:461161/1‑101 (MQ=11)
                                                                                                                                                                                     aGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCg                                                                                                                                                                              >  2:208007/1‑76 (MQ=40)
                                                                                                                                                                                      gcgcGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGCATACCTACTCGGCGTGAATGCCCACATCcaagaagcaagaagc                                              >  1:351929/1‑202 (MQ=255)
                                                                                                                                                                                      gcgcGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCt               >  2:603895/1‑233 (MQ=255)
                                                                                                                                                                                      gcgcGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGACGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAAt                                                                      >  1:619913/1‑178 (MQ=255)
                                                                                                                                                                                            tcctccCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGAc                                                                                                                                                                            >  2:22368/1‑71 (MQ=255)
                                                                                                                                                                                            tcctccCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGgt                                                                            >  2:157524/1‑166 (MQ=255)
                                                                                                                                                                                            tcctccCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGAc       >  1:256560/1‑235 (MQ=255)
                                                                                                                                                                                                cccGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGc  <  1:940186/236‑1 (MQ=255)
                                                                                                                                                                                                  cGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAgc                                                                                                                                                                                                      <  1:499644/39‑1 (MQ=40)
                                                                                                                                                                                                    ggAGCAAGCCCGGAGGGAGGGGACCCGGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGTTTTTTGGTTCCTTTTTGCCACAACAAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTa                                                                                  <  2:866855/153‑1 (MQ=18)
                                                                                                                                                                                                    ggAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtga                                                                                             <  2:784371/141‑1 (MQ=255)
                                                                                                                                                                                                      aGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                          >  2:871209/1‑79 (MQ=255)
                                                                                                                                                                                                      aGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                          <  1:871210/79‑1 (MQ=11)
                                                                                                                                                                                                      aGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                          <  1:1029823/79‑1 (MQ=11)
                                                                                                                                                                                                      aGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                          >  2:1029822/1‑79 (MQ=255)
                                                                                                                                                                                                             ccGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATc                                                             >  2:563668/1‑164 (MQ=255)
                                                                                                                                                                                                                   |                                                                                                                                                                                                                         
ATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGC  >  AE015451/2256057‑2256483

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