Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,457,128 T→A 5.2% M678L (ATG→TTG)  fadJ ← enoyl‑CoA hydratase/epimerase and isomerase/3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,457,1280TA 5.2% 458.3 / 8.5 128M678L (ATG→TTG) fadJenoyl‑CoA hydratase/epimerase and isomerase/3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (59/62);  new base A (4/3);  total (63/65)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 135 total aligned reads displayed.

ACCAGAGAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTC  >  NC_000913/2456861‑2457396
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                            
aCCAGAGAAAGAAAACGGTGGCCGGCATTTTGCCTTATCCCTTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCa                                                                                                                                                                                                                                                                       <  4:428845/275‑1 (MQ=255)
              aCGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACCTTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAATTCGTTTCATTATTGCAACCACTTCGACcgc                                                                                                                                                                                                                                                         >  4:392321/1‑275 (MQ=255)
                     ccGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCACTATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgaga                                                                                                                                                                                                                                                  >  2:117304/1‑275 (MQ=255)
                       ggCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgagaga                                                                                                                                                                                                                                                <  2:518023/275‑1 (MQ=255)
                            tgtgCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGa                                                                                                                                                                                                                                           <  2:102914/275‑1 (MQ=255)
                              tgCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGACCCACACTCCGTGGCAAGCAGTTGCATTATTGCAACCACTACGCCCGCGCCGAGAGAATCGat                                                                                                                                                                                                                                          >  4:470218/1‑274 (MQ=255)
                                          tCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGt                                                                                                                                                                                                                             >  2:545199/1‑275 (MQ=255)
                                          tCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGt                                                                                                                                                                                                                             <  2:299504/275‑1 (MQ=255)
                                                aaTGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCAATCTCGACCAAACGCTCGCAAGGGGTAATACGGGAACCATACTGCCTGCCAAGTCGTTGCCTTATTCCAACCACTTCGCCCGCCCACAAGGAATCGAAATAGCGAAACGGTCCAcag                                                                                                                                                                                                                       >  4:20751/1‑273 (MQ=255)
                                                                              agagacaggTCAATATGAGTTTACATTTTTACAAACAATCCGCCATTCAGCGCGGATTCATATAGCTTTGACCGTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCACAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATAcc                                                                                                                                                                                         <  4:223079/267‑1 (MQ=255)
                                                                                 gTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACAAAACGCTCGCAAGTGGTAAAACGGGAACCATACGGCGTGGCAAGTCGTTGCATAATTGCAACCACTTCGCACGCGCCGAGAAAATAGACTTACAGGAACGCACCACCAACAAAGGGGGAAAAACCAAAGGCaaaaaacacc                                                                                                                                                                                      >  4:514398/1‑268 (MQ=255)
                                                                                                 tttACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCaaa                                                                                                                                                                                             >  2:223080/1‑252 (MQ=255)
                                                                                                   tACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCAGGCAAATCCATTTCACGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCACAATAGCGAGCCCATAAccccaccact                                                                                                                                                                     >  4:268046/1‑268 (MQ=255)
                                                                                                     cATTTTGCCCAACACACCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCGTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTAGTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTcacgca                                                                                                                                                                   <  2:494105/274‑1 (MQ=255)
                                                                                                       ttttGCCAAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCacgcacgc                                                                                                                                                                <  2:353559/275‑1 (MQ=255)
                                                                                                        tgtGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGGCAGTTGCAGTTGTGTTCCAAAAACTTTCCCCACGCGCGCCCAACTCGACCAAACGCTCGCAAGGGGTAAAACGCGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGGGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCt                                                                                                                                                               <  2:471879/273‑1 (MQ=255)
                                                                                                                acacTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCa                                                                                                                                                                                                     <  1:3234/229‑1 (MQ=255)
                                                                                                                 cacTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGGCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGAGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCgcgc                                                                                                                                                                                                                                                       >  2:139185/1‑178 (MQ=255)
                                                                                                                        ccATTCAGCGCGGATTCATATAGCTTTGACCTTCGTATGGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCtca                                                                                                                                               <  2:268047/275‑1 (MQ=255)
                                                                                                                         cATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCtcat                                                                                                                                              <  2:459715/275‑1 (MQ=255)
                                                                                                                             cAGCGCGGATTCATATAGCTTTGACCTTCTTAATGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCGCCACGCGCGCCCATCTCGAACAAACGCTCGCAAGGGGTAAAACCGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCAGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAac                                                                                                                                          <  2:514399/275‑1 (MQ=255)
                                                                                                                                  cgGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAACACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCCCCGAGAGAATCCATATACCGGAACGGTCAACCCAGAAATGGCGGAAAACCAATGCCAAATCCCGCGCCACTACCCCCGTTCCGCCACCTACCGAAAACCTGCCAATCAacctca                                                                                                                                      >  4:377058/1‑270 (MQ=255)
                                                                                                                                     aTTCATATAGCTTTGACCTTCTGATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTa                                                                                                                                  <  2:392322/275‑1 (MQ=255)
                                                                                                                                      ttCTNNTAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAANTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCa                                                                                                                                                                                                     <  4:545198/207‑1 (MQ=255)
                                                                                                                                              gCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCattat                                                                                                                                                                                                                                                                           <  7:1040896/129‑1 (MQ=255)
                                                                                                                                                 ttGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCAATAATGCaa                                                                                                                                                                                                                                                                      <  6:631165/131‑1 (MQ=255)
                                                                                                                                                 ttGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCAATAATGCaa                                                                                                                                                                                                                                                                      >  5:631165/1‑131 (MQ=255)
                                                                                                                                                       ttcttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACtt                                                                                                                                                                                                                                                                >  6:769027/1‑131 (MQ=255)
                                                                                                                                                       ttcttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACtt                                                                                                                                                                                                                                                                <  5:769027/131‑1 (MQ=255)
                                                                                                                                                         cttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATAACCCCGTCACCCACGCTACGGATAACCCGCTCCACAACACCACATACTACTTCATTCCGCATaaa                                                                                                               >  4:263855/1‑274 (MQ=255)
                                                                                                                                                              tGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCa                                                                                                                                                                                                                                                                   <  5:21249/121‑1 (MQ=255)
                                                                                                                                                              tGCAGGTCAGTTGCAGTTGTGTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCa                                                                                                                                                                                                                                                                   >  6:21249/1‑121 (MQ=255)
                                                                                                                                                                  ggTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTc                                                                                                                                                                                                                                                               <  7:694296/121‑1 (MQ=255)
                                                                                                                                                                  ggTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTc                                                                                                                                                                                                                                                               >  8:694296/1‑121 (MQ=255)
                                                                                                                                                                       gttgcagttgTTTTCCAAAGACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGc                                                                                                                                                                                                                                                                        <  6:480774/107‑1 (MQ=255)
                                                                                                                                                                       gttgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGc                                                                                                                                                                                                                                                                        >  5:480774/1‑107 (MQ=255)
                                                                                                                                                                         tgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCattat                                                                                                                                                                                                                                                                           <  7:637315/102‑1 (MQ=255)
                                                                                                                                                                         tgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCattat                                                                                                                                                                                                                                                                           >  8:637315/1‑102 (MQ=255)
                                                                                                                                                                         tgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCg                                                                                                                                                                                                                                            <  6:1042190/133‑1 (MQ=255)
                                                                                                                                                                         tgcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGt                                                                                                                                                                                                                             >  8:1019655/1‑148 (MQ=255)
                                                                                                                                                                               tgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGa                                                                                                                                                                                                                                           >  5:460593/1‑128 (MQ=255)
                                                                                                                                                                               tgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGa                                                                                                                                                                                                                                           <  6:460593/128‑1 (MQ=255)
                                                                                                                                                                               tgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCg                                                                                                                                                                                                                       >  8:879829/1‑148 (MQ=255)
                                                                                                                                                                                 ttttCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCgag                                                                                                                                                                                                                     >  6:840297/1‑148 (MQ=255)
                                                                                                                                                                                       aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATgg                                                                                                                                                                                                               >  6:676236/1‑148 (MQ=255)
                                                                                                                                                                                       aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATgg                                                                                                                                                                                                               <  7:879829/148‑1 (MQ=255)
                                                                                                                                                                                         aaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGaaa                                                                                                                                                                                                                  <  7:317306/143‑1 (MQ=255)
                                                                                                                                                                                         aaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGaaa                                                                                                                                                                                                                  >  8:317306/1‑143 (MQ=255)
                                                                                                                                                                                            cTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGaa                                                                                                                                                                                                          <  7:1019655/148‑1 (MQ=255)
                                                                                                                                                                                             tttCCCCACGCGCGCACATCTCGACCAAACGCGCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAtt                                                                          <  1:100304/275‑1 (MQ=255)
                                                                                                                                                                                             gttCCCCACGCGAGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACAGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAt                                                                           <  1:73038/273‑1 (MQ=255)
                                                                                                                                                                                                 cccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCg                                                                                                                                                                                                                                            <  7:681626/109‑1 (MQ=255)
                                                                                                                                                                                                 cccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCg                                                                                                                                                                                                                                            >  8:681626/1‑109 (MQ=255)
                                                                                                                                                                                                     cgcgcgCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGc                                                                                                                                                                                                 <  5:676236/148‑1 (MQ=255)
                                                                                                                                                                                                     cgcgcgCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACCCAACGTACTGCTTCATTCAGCATCAACATCCCACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTg                                                                  >  2:427509/1‑275 (MQ=255)
                                                                                                                                                                                                           cccATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATa                                                                                                                                                                                           >  8:401974/1‑148 (MQ=255)
                                                                                                                                                                                                             cATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCgag                                                                                                                                                                                                                     <  8:721417/120‑1 (MQ=255)
                                                                                                                                                                                                             cATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCgag                                                                                                                                                                                                                     >  7:721417/1‑120 (MQ=255)
                                                                                                                                                                                                             cATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCCTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGACAGAATCGATATAGCGGAACGTTCCACCGAGAAATGGCGTAAAACCAATGCCAAATAcg                                                                                                                                                                                         >  6:216447/1‑147 (MQ=255)
                                                                                                                                                                                                                ctcGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGa                                                                                                                                                                                                                                  <  6:155333/104‑1 (MQ=255)
                                                                                                                                                                                                                ctcGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGa                                                                                                                                                                                                                                  >  5:155333/1‑104 (MQ=255)
                                                                                                                                                                                                                  cGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGcc                                                                                                                                                                                    >  6:1071141/1‑148 (MQ=255)
                                                                                                                                                                                                                    aCCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCAc                                                                                                                                                                                                                         >  8:178612/1‑109 (MQ=255)
                                                                                                                                                                                                                    aCCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCAc                                                                                                                                                                                                                         <  7:178612/109‑1 (MQ=255)
                                                                                                                                                                                                                      cAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAata                                                                                                                                                                                <  7:401974/148‑1 (MQ=255)
                                                                                                                                                                                                                       aaaCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCaa                                                                                                                                                                                              <  7:313235/133‑1 (MQ=255)
                                                                                                                                                                                                                         aCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCgag                                                                                                                                                                                                                     >  6:370905/1‑108 (MQ=255)
                                                                                                                                                                                                                         aCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATcc                                                                                                                                                                             <  5:1071141/148‑1 (MQ=255)
                                                                                                                                                                                                                                 aaGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCa                                                                                                                                                                                                     <  5:217423/116‑1 (MQ=255)
                                                                                                                                                                                                                                       tAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCt                                                                                                                                                               >  5:832010/1‑148 (MQ=255)
                                                                                                                                                                                                                                         aaaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTAc                                                                                                                                                             <  6:832010/148‑1 (MQ=255)
                                                                                                                                                                                                                                          aaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGc                                                                                                                                                                                                 <  6:680728/111‑1 (MQ=255)
                                                                                                                                                                                                                                          aaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGc                                                                                                                                                                                                 >  5:680728/1‑111 (MQ=255)
                                                                                                                                                                                                                                           aCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACgg                                                                                                                                                           >  6:536179/1‑148 (MQ=255)
                                                                                                                                                                                                                                            cGGGAACCATCCTGCGTGGCAAGTCGTTTCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGAc                           <  2:88088/275‑1 (MQ=255)
                                                                                                                                                                                                                                            cGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGa                                                                                                                                                          <  5:506759/148‑1 (MQ=255)
                                                                                                                                                                                                                                            cGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGa                                                                                                                                                          >  6:506759/1‑148 (MQ=255)
                                                                                                                                                                                                                                             gggAACCATACTGCGTGGCAAGTCGGTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGAGATAGAGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGAc                           <  4:171013/274‑1 (MQ=255)
                                                                                                                                                                                                                                              ggAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATa                                                                                                                                                        <  5:536179/148‑1 (MQ=255)
                                                                                                                                                                                                                                               gAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCgag                                                                                                                                                                                                                     >  3:79909/1‑86 (MQ=255)
                                                                                                                                                                                                                                               gAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTCAATGGCGGGACCgcact                         >  1:102041/1‑270 (MQ=255)
                                                                                                                                                                                                                                               gAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAAAACAACGCACTGCTTCATTCAGCATCAACATCCCACANCGTTCAGCAACCGTCGGTGCGGATATTCGCCCCTGCCTCTGTCGTCCAATCAGCGCGTAACTGGCggggaagtcct                         >  3:127758/1‑266 (MQ=255)
                                                                                                                                                                                                                                                  ccATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTcacgca                                                                                                                                                                   >  7:83545/1‑133 (MQ=255)
                                                                                                                                                                                                                                                  ccATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTcacgca                                                                                                                                                                   <  8:83545/133‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCt                                                                                                                                                    >  5:398444/1‑148 (MQ=255)
                                                                                                                                                                                                                                                   cATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTg                                                                                                                                                   <  5:216447/148‑1 (MQ=255)
                                                                                                                                                                                                                                                   cATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttcattc                                                                                                                       <  1:41304/176‑1 (MQ=255)
                                                                                                                                                                                                                                                    aTACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGACTCGATATAGCGGAACGGTCCACCGAGAAATGGCGGaaaa                                                                                                                                                                                                        <  8:663060/94‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttcagggTGGCAAGTGGTTGCAATATTGCGAACGCTGCGCCCGCGCCAAGAGATTCGATATGGCGGAACGGACCACGGAGAAAGGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGttttt                   <  4:269915/268‑1 (MQ=255)
                                                                                                                                                                                                                                                            tGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCtttt            >  2:349565/1‑274 (MQ=255)
                                                                                                                                                                                                                                                              gCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGaa                                                                                                                                                                                                                                 >  7:444495/1‑59 (MQ=255)
                                                                                                                                                                                                                                                              gCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAAcaca                                                                                                                                       <  2:335050/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tggTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATg                                                                                                                                                                                                                <  8:265186/69‑1 (MQ=38)
                                                                                                                                                                                                                                                                   tggTTGCAATAATGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGt                                                                                                                                   >  8:135299/3‑148 (MQ=255)
                                                                                                                                                                                                                                                                   tCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCt                                                                                                                                                 >  6:595322/1‑134 (MQ=255)
                                                                                                                                                                                                                                                                   tCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGAATAACCTGCt                                                                                                                                                 <  5:595322/134‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAAc                                                                                                                                     >  8:41214/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                    cGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTa                                                                                                                                  <  5:840297/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTc  <  1:10441/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ttGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACt                                                                                                                                >  8:656344/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                            
ACCAGAGAAAGAAAACGGTGGCCGGCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTC  >  NC_000913/2456861‑2457396

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