Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,087,428 G→A 7.9% intergenic (+241/+677) PP_1865 → / ← aroH transposase/Phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,087,4280GA 7.9% 199.3 / 3.2 61intergenic (+241/+677)PP_1865/aroHtransposase/Phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base A (3/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTC  >  AE015451/2087208‑2087658
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                      
ggTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTa                                                                                                                                                                                                                         <  1:704740/236‑1 (MQ=255)
   cccTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGa                                                                                                                                                                                                                       <  1:1149858/235‑1 (MQ=255)
        tttGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCTCGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCTCGCTTAGAAGc                                                                                                                                                                                                                    <  2:1185753/233‑1 (MQ=255)
                                   ccGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCGCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCgg                                                                                                                                                                                         <  2:539664/233‑1 (MQ=255)
                                     gTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGa                                                                                                                                                                                      <  2:1075731/234‑1 (MQ=255)
                                          cTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcggccg                                                                                                                                                                               >  1:110000/1‑236 (MQ=255)
                                                            tACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCg                                                                                                                                                             <  1:973706/236‑1 (MQ=255)
                                                               gCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCgg                                                                                                                                                            >  2:847790/1‑234 (MQ=255)
                                                               gCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCg                                                                                                                                                          <  1:1069850/236‑1 (MQ=255)
                                                                 tGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCAccgcc                                                                                                                                                                       >  2:1014113/1‑221 (MQ=255)
                                                                  gcgcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAg                                                                                                                                                        <  1:945196/235‑1 (MQ=255)
                                                                        cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGc                                                                                                                                                 <  1:847795/236‑1 (MQ=255)
                                                                             tCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGt                                                                                                                                            <  1:1001930/236‑1 (MQ=255)
                                                                               cccTCGCGCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGTGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTTGCCTGAACTACCGCGCACTTAGAAGCACGATCaaaagcgaaaagcgctta                                                                                                                                                                                            <  1:76803/186‑5 (MQ=17)
                                                                                                  cccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGCCCTCCTGAAGTCGCAATCTGCGTTGCCTGAACTCCcgcgcc                                                                                                                                                                                                                             >  2:115054/1‑133 (MQ=255)
                                                                                                   ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTTGCCTGAACTACCGCGCACTTAGAAGCACGATCAAAAGCgaa                                                                                                                                                                                                     >  2:1147853/1‑156 (MQ=14)
                                                                                                             gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGc                                                                                                            <  1:1105551/236‑1 (MQ=255)
                                                                                                                cTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCt                                                                                                                                                   <  1:47360/194‑1 (MQ=255)
                                                                                                                            tccatGGGACCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTa                                                                                             >  1:709336/1‑236 (MQ=255)
                                                                                                                              catGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGGAGTCGCCATCTGCGTCGCCTGGACTACCgcgcgcct                                                                                                                                                                                                                          >  1:500950/1‑107 (MQ=255)
                                                                                                                                  gggCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCTGGCATGGCTGCCCTCGGCCt                                                                                                   <  2:789268/224‑1 (MQ=255)
                                                                                                                                    gCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCCAATTCGGTGATCGGCCGGCACCGCCGATg                                                                                                                                                                   >  1:544638/1‑158 (MQ=255)
                                                                                                                                     ccTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCt                                                                                                              >  1:1228635/1‑210 (MQ=255)
                                                                                                                                         aGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCt                                                                                                   <  2:1146402/217‑1 (MQ=255)
                                                                                                                                           cgcTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCCCTTAGAAGCAAGATCAAAAGCGAGAAACGAATTCGGTGa                                                                                                                                                                                      >  1:1207569/1‑132 (MQ=255)
                                                                                                                                                  cGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGt                                                                                                                                                                                        >  1:536491/1‑123 (MQ=255)
                                                                                                                                                     gCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTTGCCTGAACTAcc                                                                                                                                                                                                                                  <  2:153610/78‑1 (MQ=255)
                                                                                                                                                     gCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTAcc                                                                                                                                                                                                                                  >  1:153613/1‑78 (MQ=255)
                                                                                                                                                           aTGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGttgctcttgcgctt                                                                >  2:487161/1‑234 (MQ=255)
                                                                                                                                                               ccGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCCTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGt                                                                                                                                            <  1:748559/154‑1 (MQ=255)
                                                                                                                                                                     tctcCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCt                                                                                                                                                   <  2:543417/141‑1 (MQ=255)
                                                                                                                                                                      ctccctccGCAGTGCCTCCCCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCa                                                                                                                   <  2:83107/172‑1 (MQ=255)
                                                                                                                                                                      ctccctccGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTc                                                                                                                                           >  1:948725/1‑148 (MQ=255)
                                                                                                                                                                      ctccctccGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTc                                                                                                                                           <  2:948720/148‑1 (MQ=255)
                                                                                                                                                                      ctccctccGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGgtgt                                                                                              >  1:739057/1‑193 (MQ=255)
                                                                                                                                                                      ctccctccGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCCTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGAGCGGCCGGCACCGCCGATGCAGTCGGCGAGGCTCTGCCCg                                                                                                                                              >  1:1113032/1‑145 (MQ=255)
                                                                                                                                                                         cctccGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTgg                                                   <  2:709331/233‑1 (MQ=255)
                                                                                                                                                                            ccGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGccc                                                                                                          >  2:688515/1‑175 (MQ=255)
                                                                                                                                                                              gCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGgaa                                                                                      <  2:109998/193‑1 (MQ=255)
                                                                                                                                                                              gCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTg                                             >  1:1250687/1‑234 (MQ=255)
                                                                                                                                                                                aaTGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTgg                                                                                                 >  2:1276504/1‑180 (MQ=255)
                                                                                                                                                                                aaTGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGGCGCCTGAACTACCGCGCGCTTAGAAGCa                                                                                                                                                                                                                   >  2:230580/1‑66 (MQ=255)
                                                                                                                                                                                 aTGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGatat                                                                                                                           <  2:1228630/153‑1 (MQ=255)
                                                                                                                                                                                          cTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCCTCTGGATCAGGATAGCGCCGGCATGGCTGCCCTCCGc                                                                                                     >  2:178119/1‑166 (MQ=255)
                                                                                                                                                                                          cTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGttgctcttgctcttg                                                               >  2:990924/1‑204 (MQ=255)
                                                                                                                                                                                          cTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGttgctcttgctcttg                                                               <  1:990929/204‑1 (MQ=255)
                                                                                                                                                                                                 tcctGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTc                                         >  2:421261/1‑219 (MQ=255)
                                                                                                                                                                                                 tcctGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCGCTTGCTCTTGCTTTgg                                                   >  2:558547/1‑209 (MQ=255)
                                                                                                                                                                                                     gAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCt                                                                                                                                                                                                                           >  2:1123740/1‑37 (MQ=255)
                                                                                                                                                                                                     gAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCt                                                                                                                                                                                                                           <  1:1123745/37‑1 (MQ=255)
                                                                                                                                                                                                            cAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGgtgt                                                                                              >  2:445897/1‑155 (MQ=255)
                                                                                                                                                                                                            cAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGgtgt                                                                                              <  1:445902/155‑1 (MQ=255)
                                                                                                                                                                                                            cAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGCGAt                                                                                                                                                                                     >  1:130619/1‑68 (MQ=255)
                                                                                                                                                                                                                cTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGtt                                                                            >  2:617930/1‑169 (MQ=255)
                                                                                                                                                                                                                   cGTCGCCTGGACTCCCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAgc                                    <  1:565945/206‑1 (MQ=255)
                                                                                                                                                                                                                    gTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTaa                                      <  2:983861/203‑1 (MQ=255)
                                                                                                                                                                                                                      cGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcggcc                                                                                                                                                                                >  1:1078647/1‑63 (MQ=255)
                                                                                                                                                                                                                      cGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCTCGGTCATCTGGTTAAggtc                                                                                                                             >  1:332702/1‑112 (MQ=255)
                                                                                                                                                                                                                      cGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGGGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGgtgt                                                                                              >  1:172299/1‑145 (MQ=255)
                                                                                                                                                                                                                       gCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCAGCTGGATAAGGATATCGCCGGCa                                                                                                                   <  2:929452/123‑1 (MQ=255)
                                                                                                                                                                                                                        ccTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTc  <  1:1014118/235‑1 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                      
GGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTC  >  AE015451/2087208‑2087658

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: