Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,137,702 G→A S2F (TCC→TTC)  sugE ← quaternary ammonium compound‑resistance protein SugE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,137,7020GA95.2% 1268.5 / 34.2 375S2F (TCC→TTC) sugEquaternary ammonium compound‑resistance protein SugE
Reads supporting (aligned to +/- strand):  ref base G (10/8);  new base A (169/188);  total (179/196)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 375 total aligned reads displayed.

TTCATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGAAAAAACAA  >  DH1/4137555‑4137849
                                                                                                                                                   |                                                                                                                                                   
ttCATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAAc                                                                                                                                                   <  4:306108/150‑1 (MQ=255)
ttCATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAAc                                                                                                                                                   <  1:101372/150‑1 (MQ=255)
 acATAGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCAGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTCGGTATATATCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAAACCAGAACa                                                                                                                                                  >  4:1005789/2‑150 (MQ=255)
     ggCCCAGGCAAGTAGCGCCATACTGACAATCATGGCCGTCACAGTAATAACAGTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCa                                                                                                                                             <  6:1343022/150‑1 (MQ=255)
     cGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCa                                                                                                                                             <  4:197924/151‑1 (MQ=255)
      gCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAg                                                                                                                                            <  6:1264108/151‑1 (MQ=255)
      gCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATAACAg                                                                                                                                            <  8:362818/151‑1 (MQ=255)
      gCCCAGGAAAGTAGCGCCATACTGACAATCATCGCCGTCAGAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAg                                                                                                                                            <  7:723457/151‑1 (MQ=255)
        ccAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGc                                                                                                                                          >  1:1331310/1‑151 (MQ=255)
         cAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCCTACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAAAATATCAGGCt                                                                                                                                         >  3:137368/1‑151 (MQ=255)
         cAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCt                                                                                                                                         >  8:1241076/1‑151 (MQ=255)
          aGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTc                                                                                                                                        <  1:894634/151‑1 (MQ=255)
          aGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTc                                                                                                                                        <  1:276549/151‑1 (MQ=255)
          aGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTc                                                                                                                                        >  3:899323/1‑151 (MQ=255)
            gCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCg                                                                                                                                      >  5:1371955/1‑151 (MQ=255)
               aGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGgag                                                                                                                                   <  7:800219/151‑1 (MQ=255)
               aGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGgag                                                                                                                                   <  2:455163/151‑1 (MQ=255)
                  aGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAAt                                                                                                                                <  5:939321/151‑1 (MQ=255)
                     gcCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTgg                                                                                                                             <  3:748436/151‑1 (MQ=255)
                     gcCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTgg                                                                                                                             <  7:297245/151‑1 (MQ=255)
                       cATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTgggg                                                                                                                           >  5:472190/1‑151 (MQ=255)
                         tACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGcc                                                                                                                         >  1:1212486/1‑151 (MQ=255)
                          aCTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGcc                                                                                                                         <  1:1107199/150‑1 (MQ=255)
                               cAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTcccc                                                                                                                   >  1:540646/1‑151 (MQ=255)
                                aaTCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCg                                                                                                                  <  5:6454/151‑1 (MQ=255)
                                   catcGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTg                                                                                                               >  7:1228232/1‑151 (MQ=255)
                                    atcGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGa                                                                                                              <  1:45250/151‑1 (MQ=255)
                                       gCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGaaaa                                                                                                           <  2:1199134/151‑1 (MQ=255)
                                         cGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGaaa                                                                                                            >  4:1147214/1‑148 (MQ=255)
                                          gTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAag                                                                                                          >  1:314193/1‑149 (MQ=255)
                                              caGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGt                                                                                                    >  5:767138/1‑151 (MQ=255)
                                               aGTAATAACACTCGGCGTCAAACGAGTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTCGCGCCGACCCCGGTCAAAAGAGGCGtt                                                                                                   >  7:1025489/1‑151 (MQ=255)
                                               aGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGtt                                                                                                   >  5:1172608/1‑151 (MQ=255)
                                                 taataaCACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTg                                                                                                 >  2:948049/1‑151 (MQ=255)
                                                   ataaCACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGc                                                                                                >  8:244698/1‑150 (MQ=255)
                                                    taaCACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGc                                                                                                >  6:426564/1‑149 (MQ=255)
                                                     aaCACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAgg                                                                                             >  4:142716/1‑151 (MQ=255)
                                                     aaCACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAg                                                                                              <  6:922316/150‑1 (MQ=255)
                                                     aaCACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAg                                                                                              <  4:610902/150‑1 (MQ=255)
                                                          tCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGgtcgtc                                                                                        >  8:1086351/1‑151 (MQ=255)
                                                             gCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTg                                                                                     >  8:811159/1‑151 (MQ=255)
                                                              cGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGc                                                                                    >  8:57493/1‑151 (MQ=255)
                                                                 cAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCaa                                                                                  >  1:206772/1‑150 (MQ=255)
                                                                 cAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCaa                                                                                  >  3:529140/1‑150 (MQ=255)
                                                                 cAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCaa                                                                                  >  1:718265/1‑150 (MQ=255)
                                                                  aaaCGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCaa                                                                                  >  8:1360009/1‑149 (MQ=255)
                                                                  aaaCGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTACGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCaaac                                                                                >  6:1171716/1‑150 (MQ=255)
                                                                      gACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTg                                                                             >  6:792281/1‑150 (MQ=255)
                                                                         tAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATgg                                                                         >  4:1094867/1‑151 (MQ=255)
                                                                            aGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTg                                                                      <  2:80272/151‑1 (MQ=255)
                                                                             gCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGa                                                                     <  6:242867/151‑1 (MQ=255)
                                                                               cGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAAc                                                                   <  2:257723/151‑1 (MQ=255)
                                                                                gTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACg                                                                  <  1:967262/151‑1 (MQ=255)
                                                                                  gggTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGa                                                                <  5:1243162/151‑1 (MQ=255)
                                                                                   ggTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGat                                                               >  4:772541/1‑151 (MQ=255)
                                                                                   ggTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGat                                                               <  5:257153/151‑1 (MQ=255)
                                                                                    gTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCAGCAAAGCTGATGGTTGAACGGat                                                               >  5:1126396/1‑150 (MQ=255)
                                                                                     tataTTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGatat                                                             >  7:478456/1‑151 (MQ=255)
                                                                                      atatTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGatat                                                             <  6:1218737/150‑1 (MQ=255)
                                                                                      atatTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTCATGGTTGAACGGATAtt                                                            >  4:887750/1‑151 (MQ=255)
                                                                                       tatTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATAttt                                                           >  3:185224/1‑151 (MQ=255)
                                                                                       tatTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATAttt                                                           >  5:427854/1‑151 (MQ=255)
                                                                                        atTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATAtttt                                                          <  1:1272506/151‑1 (MQ=255)
                                                                                          ttCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGc                                                        >  3:688601/1‑151 (MQ=255)
                                                                                          ttCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGGGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTg                                                         <  4:692292/150‑1 (MQ=255)
                                                                                          ttCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTAGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTg                                                         <  7:695541/150‑1 (MQ=255)
                                                                                           tCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGc                                                        <  1:987484/150‑1 (MQ=255)
                                                                                            cAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCa                                                      <  3:475044/151‑1 (MQ=255)
                                                                                               gCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGt                                                   >  5:1085674/1‑151 (MQ=255)
                                                                                               gCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGt                                                   <  3:606535/151‑1 (MQ=255)
                                                                                                 cAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTgg                                                 <  3:1094833/151‑1 (MQ=255)
                                                                                                     ggCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAgg                                             <  2:335165/151‑1 (MQ=255)
                                                                                                      gCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAggg                                            <  2:355819/151‑1 (MQ=255)
                                                                                                            aCCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAAttt                                      <  7:140329/151‑1 (MQ=255)
                                                                                                            aCCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAAtt                                       <  7:855369/150‑1 (MQ=255)
                                                                                                             ccACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTc                                     <  5:1297843/151‑1 (MQ=255)
                                                                                                               aCTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCa                                    >  4:628306/1‑150 (MQ=255)
                                                                                                                 ttCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTTCCAGGTGGCAGGGGAATTTAAACt                                 >  8:764276/1‑151 (MQ=255)
                                                                                                                  tCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAAc                                  <  8:667909/149‑1 (MQ=255)
                                                                                                                  tCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAAc                                  <  5:1198533/149‑1 (MQ=255)
                                                                                                                    cagcagACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTAtt                              <  4:410013/151‑1 (MQ=255)
                                                                                                                        agACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATTTTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCt                          >  8:1001561/1‑151 (MQ=255)
                                                                                                                        agACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCt                          <  5:1010599/151‑1 (MQ=255)
                                                                                                                          aCCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCttt                        >  6:212784/1‑151 (MQ=255)
                                                                                                                           ccAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTtg                       >  6:76922/1‑151 (MQ=255)
                                                                                                                             aGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTtgtg                     >  8:1169059/1‑151 (MQ=255)
                                                                                                                              gCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGa                    <  6:693440/151‑1 (MQ=255)
                                                                                                                                aaTAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACt                  >  7:899294/1‑151 (MQ=255)
                                                                                                                                  tAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCACGGGAATTTCAACTATTACCTTTGTGACTAt                >  7:110731/1‑151 (MQ=255)
                                                                                                                                   aaCTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATc               <  5:58579/151‑1 (MQ=255)
                                                                                                                                   aaCTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATc               >  2:730242/1‑151 (MQ=255)
                                                                                                                                    aCTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCg              <  2:171391/151‑1 (MQ=255)
                                                                                                                                    aCTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCg              >  3:377888/1‑151 (MQ=255)
                                                                                                                                       aaGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAg           <  2:641440/151‑1 (MQ=255)
                                                                                                                                        aGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGc             >  3:1152970/1‑148 (MQ=255)
                                                                                                                                        aGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGc             >  4:1322401/1‑148 (MQ=255)
                                                                                                                                          aTAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGaaa        <  6:1073793/151‑1 (MQ=255)
                                                                                                                                          aTAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGaaa        <  8:907408/151‑1 (MQ=255)
                                                                                                                                               ccAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGAAAaaac    >  6:477099/1‑150 (MQ=255)
                                                                                                                                                cAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGAAAaaacaa  <  7:1324007/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTCATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGAAAAAACAA  >  DH1/4137555‑4137849

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: ‑