Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 3,935,805 C→A I255I (ATC→ATA)  mcrC → protein mcrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,935,8050CA98.6% 2419.8 / 12.3 646I255I (ATC→ATA) mcrCprotein mcrC
Reads supporting (aligned to +/- strand):  ref base C (5/4);  new base A (286/348);  total (292/354)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 646 total aligned reads displayed.

AGATGTCATTACTTTATCAAAAGTTTCTTTATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATCACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGCCTGAAAATGGC  >  DH1/3935657‑3935952
                                                                                                                                                    |                                                                                                                                                   
agaTGTCATTACTTTATCAAAAGTTTCTTTATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAAc                                                                                                                                                   <  3:879797/151‑1 (MQ=255)
   tGTCATTACTTTATCAAAAGTTTCTTTATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATATAAAATGGGATGCATCCAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATTGAAACTGACATAAcaa                                                                                                                                                 >  5:63193/1‑148 (MQ=255)
          aCTTTATCAAAAGTTTCTTTATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACNACCCGCNCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCtc                                                                                                                                          <  2:465613/150‑1 (MQ=255)
          aCTTTATCAAAAGTTTCTTTATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCtca                                                                                                                                         <  5:803480/151‑1 (MQ=255)
               aTCAAAAGTTTCTTTATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGa                                                                                                                                    <  6:111071/151‑1 (MQ=255)
                 gaaaaGTTTCTTTATGTATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGaa                                                                                                                                   <  1:204909/149‑1 (MQ=255)
                  caaaGTTTCTTTATGAATTTTGCCGTCGGGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGaaaa                                                                                                                                 <  8:1319439/150‑1 (MQ=255)
                  aaaaGTTTCTTTATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGaaaa                                                                                                                                 <  6:1847603/151‑1 (MQ=255)
                  aaaaGTTTCTTTATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGaaa                                                                                                                                  >  3:356186/1‑150 (MQ=255)
                  aaaaGTTTCTTAATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGAAGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGaaaa                                                                                                                                 <  8:1437459/151‑1 (MQ=255)
                   aaaGTTTCTTTATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGaaaa                                                                                                                                 >  6:1025694/1‑150 (MQ=255)
                           tttATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTAt                                                                                                                        >  6:1651322/1‑151 (MQ=255)
                            ttATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATCACCATTCGCTCATCAGAAAAAATACtt                                                                                                                          <  7:1511317/148‑1 (MQ=255)
                            ttATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACtt                                                                                                                          <  7:1117395/148‑1 (MQ=255)
                            ttATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACtt                                                                                                                          <  4:941461/148‑1 (MQ=255)
                            ttATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTa                                                                                                                         <  3:1592355/149‑1 (MQ=255)
                            ttATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATc                                                                                                                       >  1:371363/1‑151 (MQ=255)
                             tATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATc                                                                                                                       <  5:470657/150‑1 (MQ=255)
                             tATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCg                                                                                                                      >  2:855061/1‑151 (MQ=255)
                              aTGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCg                                                                                                                      <  7:44853/150‑1 (MQ=255)
                               tGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGtt                                                                                                                    >  1:1808738/1‑151 (MQ=255)
                                gAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGtt                                                                                                                    <  3:814431/150‑1 (MQ=255)
                                 aaTTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGa                                                                                                                  <  3:1887390/151‑1 (MQ=255)
                                  aTTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGAc                                                                                                                 <  4:1957614/151‑1 (MQ=255)
                                  aTTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGAc                                                                                                                 <  2:76773/151‑1 (MQ=255)
                                   ttttGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATCACCATTCGCTCATCAGAAAAAATACTTATCGTTGACg                                                                                                                >  6:1529582/1‑151 (MQ=255)
                                   ttttGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACg                                                                                                                >  1:1628556/1‑151 (MQ=255)
                                    tttGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACg                                                                                                                <  7:1677945/150‑1 (MQ=255)
                                      tGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGcc                                                                                                              <  1:489752/150‑1 (MQ=255)
                                        ccGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCaaa                                                                                                           <  4:606519/151‑1 (MQ=255)
                                           tcgtGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATAc                                                                                                        >  2:1923638/1‑151 (MQ=255)
                                            cgtGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGCCGCCAAATACt                                                                                                       >  1:896548/1‑151 (MQ=255)
                                             gtGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACta                                                                                                      <  6:1388657/151‑1 (MQ=255)
                                              tGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACtat                                                                                                     >  2:1017812/1‑151 (MQ=255)
                                              tGAATTAACCTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGACACTTAAATTGTTACATCGAATGGAAACTGACAAAACCAATCGCTCAACAGAAAAAATAATTATCGTAGACGCCAAAAACta                                                                                                      >  4:1045348/1‑150 (MQ=255)
                                               gAATTAACGTCTGCAAACACTAGCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACtat                                                                                                     <  8:1659759/150‑1 (MQ=255)
                                               gAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACtat                                                                                                     <  6:1301244/150‑1 (MQ=255)
                                                 aTTAACGTCTGCAAACACAACCCGCTCTTATTTATAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATCACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAag                                                                                                  <  1:1449727/151‑1 (MQ=255)
                                                   tAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAag                                                                                                  >  2:1600747/1‑149 (MQ=255)
                                                     aCGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCa                                                                                              <  5:21978/151‑1 (MQ=255)
                                                     aCGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCa                                                                                              <  1:100784/151‑1 (MQ=255)
                                                      cGTCTGCAAACACAAGCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATTGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCAt                                                                                             <  2:1311808/151‑1 (MQ=255)
                                                        tCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCAttt                                                                                           >  2:1195635/1‑151 (MQ=255)
                                                          tGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCAttttt                                                                                         >  5:1747491/1‑151 (MQ=255)
                                                           gCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCAtttt                                                                                          <  6:461804/149‑1 (MQ=255)
                                                           gCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCAttt                                                                                           <  6:313596/148‑1 (MQ=255)
                                                             aaaCACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCa                                                                                      <  4:74406/151‑1 (MQ=255)
                                                              aaCACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCac                                                                                     <  7:839396/151‑1 (MQ=255)
                                                              aaCACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCac                                                                                     <  7:544535/151‑1 (MQ=255)
                                                               acacAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCacg                                                                                    <  1:1352234/151‑1 (MQ=255)
                                                                 acaACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCacgac                                                                                  <  5:552679/151‑1 (MQ=255)
                                                                 acaACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCacgac                                                                                  <  7:956041/151‑1 (MQ=255)
                                                                  aaaCCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCacgacg                                                                                 <  4:318325/150‑1 (MQ=255)
                                                                     cccGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAAt                                                                              <  7:1328723/151‑1 (MQ=255)
                                                                      ccGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATg                                                                             <  8:683643/151‑1 (MQ=255)
                                                                        gCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATggg                                                                           >  5:433574/1‑151 (MQ=255)
                                                                        gCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATggg                                                                           <  8:1256947/151‑1 (MQ=255)
                                                                        gCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATggg                                                                           >  2:1871331/1‑151 (MQ=255)
                                                                           ctTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGa                                                                          <  5:890627/149‑1 (MQ=255)
                                                                              aTTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAg                                                                      <  7:658806/150‑1 (MQ=255)
                                                                                 tAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGaaa                                                                   <  7:1935860/150‑1 (MQ=255)
                                                                                  aaaaTGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGaaaa                                                                  >  4:1372156/1‑150 (MQ=255)
                                                                                    aaTGGGATTCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGTATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAAtt                                                               <  4:423415/151‑1 (MQ=255)
                                                                                     aTGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAAttt                                                              <  4:494125/151‑1 (MQ=255)
                                                                                       gggATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTc                                                             <  3:1157665/150‑1 (MQ=255)
                                                                                       gggATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTc                                                             <  3:1706306/150‑1 (MQ=255)
                                                                                          aTGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTc                                                         >  5:1880831/1‑151 (MQ=255)
                                                                                            gCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGc                                                       <  2:934665/151‑1 (MQ=255)
                                                                                            gCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGc                                                       >  5:285499/1‑151 (MQ=255)
                                                                                             cATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGc                                                       >  5:1407530/1‑150 (MQ=255)
                                                                                              aTCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCa                                                      >  8:478746/1‑150 (MQ=255)
                                                                                               tCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCaa                                                     >  1:834156/1‑150 (MQ=255)
                                                                                               tCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCaa                                                     >  5:1619260/1‑150 (MQ=255)
                                                                                                cGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCaaaa                                                   <  8:224561/151‑1 (MQ=255)
                                                                                                cGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCaaaa                                                   <  4:1532000/151‑1 (MQ=255)
                                                                                                cGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCaaa                                                    <  1:1173938/150‑1 (MQ=255)
                                                                                                cGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCaaa                                                    >  3:225546/1‑150 (MQ=255)
                                                                                                          ggATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCa                                          >  3:592479/1‑150 (MQ=255)
                                                                                                           gATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCaa                                         >  8:906017/1‑150 (MQ=255)
                                                                                                              cAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACtga                                     <  3:1509939/151‑1 (MQ=255)
                                                                                                                  cACTTAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACtgatga                                  >  7:331418/1‑150 (MQ=255)
                                                                                                                     ttAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAAtt                               >  3:1356759/1‑150 (MQ=255)
                                                                                                                      tAATTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTa                              >  6:972807/1‑150 (MQ=255)
                                                                                                                       aaTTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACt                            <  8:1842871/151‑1 (MQ=255)
                                                                                                                        aTTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACtt                           <  7:1625039/151‑1 (MQ=255)
                                                                                                                        aTTTGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACtt                           >  2:582325/1‑151 (MQ=255)
                                                                                                                          ttGTTACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTAt                         >  6:1200219/1‑151 (MQ=255)
                                                                                                                              tACCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGt                      <  6:1979937/150‑1 (MQ=255)
                                                                                                                               aCCTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCg                    >  2:1783520/1‑151 (MQ=255)
                                                                                                                                ccTCGAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGt                   <  5:1284565/151‑1 (MQ=255)
                                                                                                                                    gAATGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTaa                >  6:882895/1‑150 (MQ=255)
                                                                                                                                      aTGGAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGc             >  4:1797215/1‑151 (MQ=255)
                                                                                                                                        ggAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGCCt           >  5:1479108/1‑151 (MQ=255)
                                                                                                                                        ggAAACTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGCCt           >  8:785781/1‑151 (MQ=255)
                                                                                                                                         atgttccaCATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGCCTg          >  7:1752917/8‑151 (MQ=255)
                                                                                                                                          aaaCTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGCCTGa         <  8:38998/151‑1 (MQ=255)
                                                                                                                                          aaaCTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGCCTGa         <  1:686337/151‑1 (MQ=255)
                                                                                                                                            aCTGACATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGCCTGaaa       <  2:29546/151‑1 (MQ=255)
                                                                                                                                                aCATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCCAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGCCTGAAAATgg   <  6:664715/151‑1 (MQ=255)
                                                                                                                                                aCATAACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGCCTGAAAATgg   >  8:611918/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGATGTCATTACTTTATCAAAAGTTTCTTTATGAATTTTGCCGTCGTGAATTAACGTCTGCAAACACAACCCGCTCTTATTTAAAATGGGATGCATCGAGTATATCGGATCAGTCACTTAATTTGTTACCTCGAATGGAAACTGACATCACCATTCGCTCATCAGAAAAAATACTTATCGTTGACGCCAAATACTATAAGAGCATTTTTTCACGACGAATGGGAACAGAAAAATTTCATTCGCAAAATCTTTATCAACTGATGAATTACTTATGGTCGTTAAAGCCTGAAAATGGC  >  DH1/3935657‑3935952

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