Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,256,290 C→T 17.6% intergenic (+350/‑354) 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,2900CT17.6% 144.8 / 28.4 63intergenic (+350/‑354)asd/PP_1990Aspartate‑semialdehyde dehydrogenase/transposase
Reads supporting (aligned to +/- strand):  ref base C (30/21);  new base T (4/7);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-02

CGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACGATGGCGACACATACTCACCTTT  >  AE015451/2256090‑2256521
                                                                                                                                                                                                         |                                                                                                                                                                                                                                        
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                                                                                                                                                                                                                              <  2:114447/109‑1 (MQ=32)
                                                                                                        gCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGgatg                                                                                                                                                                                                                              >  1:114447/1‑109 (MQ=32)
                                                                                                        gCACGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATCACTACTCGGGGTGa                                                                                                              >  2:306873/1‑220 (MQ=255)
                                                                                                          aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAGGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGt                                                                                                                                                        <  1:774473/176‑1 (MQ=255)
                                                                                                          aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaag                                                                                            >  1:486586/1‑236 (MQ=255)
                                                                                                            ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCAc                                                                                                                                                                                           <  2:336945/140‑1 (MQ=255)
                                                                                                            ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCAc                                                                                                                                                                                           >  1:336946/1‑140 (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                                                                                                                                                        >  2:794676/1‑149 (MQ=255)
                                                                                                                                      aGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGt                                                                                                                                                        <  1:794677/149‑1 (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)
                                                                                                                                                                   ggAGCAAGCCCGGAGGGAGGGGACCCGGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGTTTTTTGGTTCCTTTTTGCCACAACAAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTa                                                                                                                        <  2:866855/153‑1 (MQ=18)
                                                                                                                                                                   ggAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtga                                                                                                                                   <  2:784371/141‑1 (MQ=255)
                                                                                                                                                                     aGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                                                                >  2:1029822/1‑79 (MQ=255)
                                                                                                                                                                     aGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                                                                <  1:1029823/79‑1 (MQ=11)
                                                                                                                                                                     aGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                                                                <  1:871210/79‑1 (MQ=11)
                                                                                                                                                                     aGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt                                                                                                                                                                                                >  2:871209/1‑79 (MQ=255)
                                                                                                                                                                            ccGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATc                                                                                                   >  2:563668/1‑164 (MQ=255)
                                                                                                                                                                                       ggACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACg                        <  1:1108575/228‑1 (MQ=255)
                                                                                                                                                                                           cccGGAGCGTAGCGCAGGGGCGGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCg                                       <  1:101747/209‑1 (MQ=255)
                                                                                                                                                                                             cGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaagcaagaagc                                                                      >  2:755356/1‑176 (MQ=255)
                                                                                                                                                                                                  cGTAGAGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACACCcaagaagcaagaagca                                                                                   <  2:815671/158‑1 (MQ=255)
                                                                                                                                                                                                    tAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGc                                        <  1:964113/199‑1 (MQ=255)
                                                                                                                                                                                                      gcgcAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACGATGGCGACACATACTCACCttt  >  1:938280/1‑235 (MQ=255)
                                                                                                                                                                                                         |                                                                                                                                                                                                                                        
CGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGCGAACAGTCCATTCACGATGGCGACACATACTCACCTTT  >  AE015451/2256090‑2256521

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: