Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,830 A→G 6.5% intergenic (+22/+18) tonB → / ← yciA membrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,8300AG 6.5% 610.8 / 17.4 187intergenic (+22/+18)tonB/yciAmembrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base A (84/88);  new base G (8/5);  total (93/94)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01

Only 100 of 188 total aligned reads displayed.

CAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACG  >  NC_000913/1311566‑1312088
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                  
     gcatgatgctgCACGAGGACAGGCAGTGCGCATTGAAGTGCAGGTTAAAGTTAAATTGGACGTTACGCCGGATGGTCGAGGTGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                                                                                                                                                                                     <  4:210361/265‑1 (MQ=255)
                                           ggCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGCAATTCAGTAAGCAGACAGTCAACAGCCTCCGACCGGAGACTTTTGACTCTTACTTACAAGGTCAGGCCCTAGGTTTTCCttgcgg                                                                                                                                                                                                               >  4:263983/1‑271 (MQ=255)
                                                gttaaagttaaaTTTGACGTTACGCCGGACGGTCGCGTGGATAACGTACAAATCCTTTCAGCCAAGCCTGCGAACATGCTTGAGCGTGGGTGAAAAATGCGATGCGTAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGGGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTTCGACCGGAGGCTTTTGACTATTACTTAACAGGTAAGGCGCGAGGTTTGCCTTCAGGATCAAc                                                                                                                                                                                                         >  1:87667/1‑275 (MQ=255)
                                                     agttaaaTTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTc                                                                                                                                                                                                                  >  4:529725/1‑262 (MQ=255)
                                                       ttaaaTTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGAc                                                                                                                                                                                                   <  4:486235/275‑1 (MQ=255)
                                                         aaaTTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGa                                                                                                                                                                                                              <  2:529725/262‑1 (MQ=255)
                                                                 cGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTGGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat                                                                                                                                                                                         <  2:133668/275‑1 (MQ=255)
                                                                     actccGGATGGTCGCGTGGATAACGTACAAATCCTCTCACCCAAGCCTGCGAACCTGTCTGAGCGTGAGGTGAACAATGCGATGCGCAGATTGCCTTATGGGCGGGTTAAGCCAGGCCGTGGGTTTCTTGTGACTTTCCTGTTTTAAATTAAAGGCACCCCCGCAATTCCGTAACCAGAAAACCCAAAGCCTCCGGCCGGCGGCTTATGTCTATTACTCAACCGGCAAGGCGCCAGGTTTTCCATCAGGACCAACCTCGACATCCTTAaataatc                                                                                                                                                                                     >  4:238326/4‑274 (MQ=255)
                                                                                                  aaTCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGCTCAACTGCCACATACTTAAATAATGCTTCCGTCGCTTTATAGCGTTGCCCaa                                                                                                                                                          >  4:151058/1‑273 (MQ=255)
                                                                                                        ctcAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTc                                                                                                                                                   >  2:324865/1‑274 (MQ=255)
                                                                                                               aaGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGGTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTa                                                                                                                                           <  2:460716/275‑1 (MQ=255)
                                                                                                                 gCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTa                                                                                                                                           >  2:97361/1‑273 (MQ=255)
                                                                                                                 gCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATAACCTGTTTAAAATTGACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAATCCTCCGACCGGAGGCTCTTGACGATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGCCGCtg                                                                                                                                           >  2:253855/1‑272 (MQ=255)
                                                                                                                     gCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATCGGTTCAGACGATACttttt                                                                                                                                     >  2:307548/1‑275 (MQ=255)
                                                                                                                        aaCATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCgg                                                                                                                                                                                                                                                                <  7:115321/149‑1 (MQ=255)
                                                                                                                         aCATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                                                                                                                                                                               <  7:820717/149‑1 (MQ=255)
                                                                                                                            tGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCca                                                                                                                              <  2:484725/275‑1 (MQ=255)
                                                                                                                              tttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtt                                                                                                                                                                                                                                                          >  8:720592/1‑149 (MQ=255)
                                                                                                                              tttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtt                                                                                                                                                                                                                                                          <  8:238758/149‑1 (MQ=255)
                                                                                                                              tttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtt                                                                                                                                                                                                                                                          >  7:248574/1‑149 (MQ=255)
                                                                                                                               ttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCttt                                                                                                                                                                                                                                                         <  8:343566/149‑1 (MQ=255)
                                                                                                                               ttGAGCGTGAGGTGAAAAATACGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCTTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCttt                                                                                                                                                                                                                                                         <  7:888186/149‑1 (MQ=255)
                                                                                                                                tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                                                                                                                                                                                        <  6:193952/149‑1 (MQ=255)
                                                                                                                                tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                                                                                                                                        <  6:605746/149‑1 (MQ=255)
                                                                                                                                 gAGCGTGAGGTGAAAAATGCGATGCGAAGATGTCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACtt                                                                                                                         <  1:164906/275‑1 (MQ=255)
                                                                                                                                   gCGTGAGGTGAAAAAGTCGATGCGCAGAGGGCGTTATGAGCCGGGTAAGCCCTGCAGTGGGATTGTGGTTAATATCCTGTTTAAAAATAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTc                                                                                                                        <  4:307548/274‑1 (MQ=255)
                                                                                                                                    cGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGACt                                                                                                                                                                                                                                                    <  7:720592/149‑1 (MQ=255)
                                                                                                                                        aGGGGAAAAATGCGATGTGTAGATGGCGTTATGAGCCAGGTAAGCCAGGCAGTGGGATTGTGGTGAATAACCTGTTTAAAATTAACGGCAGCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGCCCGGAGGCTTTTGGCTATTa                                                                                                                                                                                                                                                <  8:875472/149‑1 (MQ=255)
                                                                                                                                         ggTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGCCTATTAc                                                                                                                                                                                                                                               <  8:542301/149‑1 (MQ=255)
                                                                                                                                          gTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACt                                                                                                                                                                                                                                              >  5:750765/1‑149 (MQ=255)
                                                                                                                                            gAAAAATGCGATACGCAGATGGCGTTATGAGCCGGGTAAGCAAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCCCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAAGGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATg                                                                                                              <  2:190744/275‑1 (MQ=255)
                                                                                                                                             aaaaaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCaa                                                                                                                                                                                                                                           >  7:529061/1‑149 (MQ=255)
                                                                                                                                             aaaaaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATg                                                                                                              <  4:97361/274‑1 (MQ=255)
                                                                                                                                                aaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAg                                                                                                                                                                                                                                        <  5:467512/149‑1 (MQ=255)
                                                                                                                                                  tgcgatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAAcc                                                                                                                                                                                                        >  4:287591/1‑179 (MQ=255)
                                                                                                                                                    cgatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGAcc                                                                                                        >  2:68221/1‑273 (MQ=255)
                                                                                                                                                     gatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCATCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAg                                                                                                                                                                                                                                   >  8:641775/1‑149 (MQ=255)
                                                                                                                                                     gatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAg                                                                                                                                                                                                                                   >  6:375985/1‑149 (MQ=255)
                                                                                                                                                          gcAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTc                                                                                                                                                                                                                                             >  8:423542/1‑134 (MQ=255)
                                                                                                                                                                 gCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                                                                                                                                                                                                       >  5:936134/1‑149 (MQ=255)
                                                                                                                                                                   gTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                                                                                                                                                                                     <  6:912360/149‑1 (MQ=255)
                                                                                                                                                                   gTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                                                                                                                                                                                     >  5:801872/1‑149 (MQ=255)
                                                                                                                                                                      aTGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAACCCTCCGGTCGGAGGCtttt                                                                                                                                                                                                                                                        <  8:747504/111‑1 (MQ=255)
                                                                                                                                                                        gAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAg                                                                                                                                                                                                                <  5:586553/149‑1 (MQ=255)
                                                                                                                                                                         aGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAgg                                                                                                                                                                                                               <  5:507144/149‑1 (MQ=255)
                                                                                                                                                                                  aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTc                                                                                                                        <  1:3347/227‑1 (MQ=255)
                                                                                                                                                                                   gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                                                                                                                                                                                            >  7:775485/1‑126 (MQ=255)
                                                                                                                                                                                   gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAAAAGGTAAGGCGCGAGGTTTTACTTCAGGATCAACcgcg                                                                                                                                                                                                     >  7:280963/1‑149 (MQ=255)
                                                                                                                                                                                        gCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATagcagc                                                                   >  2:517323/1‑274 (MQ=255)
                                                                                                                                                                                         cAGTGGGATTGTGGTGAATATCCTGTTGAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAgcagca                                                                  <  2:238326/274‑1 (MQ=255)
                                                                                                                                                                                           gTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtt                                                                                                                                                                                                                         >  5:159173/1‑121 (MQ=255)
                                                                                                                                                                                                 ttGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataa                                                                                                                                                                                       <  6:673680/149‑1 (MQ=255)
                                                                                                                                                                                                   gtggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATg                                                                                                                                                                                     <  6:750765/149‑1 (MQ=255)
                                                                                                                                                                                                      gtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACa                                                                                                                                                                                                  <  4:700267/133‑1 (MQ=255)
                                                                                                                                                                                                      gtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCt                                                                                        >  2:399580/1‑239 (MQ=255)
                                                                                                                                                                                                      gtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGa                                                    >  2:581361/1‑275 (MQ=255)
                                                                                                                                                                                                       tgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc                                                                                                                                                                                 >  7:583480/1‑149 (MQ=255)
                                                                                                                                                                                                         aaTATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCAt                                                                         >  2:613031/1‑251 (MQ=255)
                                                                                                                                                                                                          atatCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGac                                                                                    <  4:399580/239‑1 (MQ=255)
                                                                                                                                                                                                              ccTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                                                                                                                                                                                        >  4:661277/1‑103 (MQ=255)
                                                                                                                                                                                                              ccTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAAc                                             >  2:560233/1‑274 (MQ=255)
                                                                                                                                                                                                               cTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCtt                                                                                                                                                                         >  8:136475/1‑149 (MQ=255)
                                                                                                                                                                                                                        aTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGtt                                                                                                                                                                >  6:840105/1‑149 (MQ=255)
                                                                                                                                                                                                                        aTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGtt                                                                                                                                                                <  5:230381/149‑1 (MQ=255)
                                                                                                                                                                                                                         ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTg                                                                                                                                                               >  6:406165/1‑149 (MQ=255)
                                                                                                                                                                                                                         ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTg                                                                                                                                                               >  5:555130/1‑149 (MQ=255)
                                                                                                                                                                                                                         ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCGTCTGTCGCTTTATAGCGTTg                                                                                                                                                               >  6:898071/1‑149 (MQ=255)
                                                                                                                                                                                                                           aaCGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGAc                                                                                                         <  4:597897/201‑1 (MQ=255)
                                                                                                                                                                                                                              ggCACCACCGAAATTCAGTAAGAAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCa                            <  4:47783/275‑1 (MQ=255)
                                                                                                                                                                                                                                caccacCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGgcg                                                                                                                                                                                                                                <  4:439558/77‑1 (MQ=255)
                                                                                                                                                                                                                                    accGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGtt                                                                                                                                                    >  8:193663/1‑149 (MQ=255)
                                                                                                                                                                                                                                     ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACtt                                                                                                                                        >  3:192142/1‑160 (MQ=255)
                                                                                                                                                                                                                                     ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCGCACTTc                                                                                                                        <  2:282962/176‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaTTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGa                                                                                                                                                <  7:193663/149‑1 (MQ=255)
                                                                                                                                                                                                                                         aaTTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACtttttt                                                                                                                                    <  1:192142/160‑1 (MQ=255)
                                                                                                                                                                                                                                           ttCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTCTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGc                                                                                                                                             >  6:324077/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGca             <  2:263983/274‑1 (MQ=255)
                                                                                                                                                                                                                                               gTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACt                                                                                                                                         >  7:160330/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   gCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTc                                                                                                                                                                                                                      <  8:415476/68‑1 (MQ=255)
                                                                                                                                                                                                                                                   gCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTc                                                                                                                                                                                                                      >  7:415476/1‑68 (MQ=255)
                                                                                                                                                                                                                                                    cAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAATAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACtttttt                                                                                                                                    <  5:323752/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     aGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTc                                                                                                                                   >  6:649362/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAg                                                                                                                                                                    >  7:478874/1‑115 (MQ=255)
                                                                                                                                                                                                                                                      gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAg                                                                                                                                                                    <  8:478874/115‑1 (MQ=255)
                                                                                                                                                                                                                                                      gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCa                                                                                                                                  <  5:876654/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCAc                                                                                                                                 >  5:128231/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        taGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCAcc                                                                                                                                >  5:519675/2‑149 (MQ=255)
                                                                                                                                                                                                                                                        taGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTCTTTCTcc                                                                                                                                >  7:858175/2‑149 (MQ=255)
                                                                                                                                                                                                                                                          gTCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCca                                                                                                                              >  6:107083/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           tCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCcac                                                                                                                             >  6:550522/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             aaaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGTCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCacac                                                                                                                           >  5:163913/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             aaaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTc                                                                                                                                   <  7:938597/141‑1 (MQ=255)
                                                                                                                                                                                                                                                             aaaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCacac                                                                                                                           >  5:377154/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              aaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCt                                                                                                                                                                                >  3:159880/1‑95 (MQ=255)
                                                                                                                                                                                                                                                               aaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACtt                                                                                                                         >  8:797533/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                aGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTc                                                                                                                        >  5:609626/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  ccTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg                                                                                                                                                                            <  1:159880/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCaa                                                                                                                      <  7:337789/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGcc                                                                                                                                                             >  8:359880/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                    tCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGcc                                                                                                                                                             <  7:359880/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                  
CAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAATGCTGACCGATGTCGTCCCTTTCTGGACACAGCGTGCATAGCAGCACACCACATCGCCGACCGCAACCGGCCGTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACG  >  NC_000913/1311566‑1312088

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