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,7020GA96.7% 1077.6 / 14.1 306S2F (TCC→TTC) sugEquaternary ammonium compound‑resistance protein SugE
Reads supporting (aligned to +/- strand):  ref base G (7/3);  new base A (127/166);  total (136/170)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

TTCATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGAAAAAACAAACG  >  DH1/4137555‑4137852
                                                                                                                                                   |                                                                                                                                                      
 tCATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACa                                                                                                                                                     <  1:837176/150‑1 (MQ=255)
 tCATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACa                                                                                                                                                     <  7:273580/150‑1 (MQ=255)
  cATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACata                                                                                                                                                   <  4:110229/151‑1 (MQ=255)
  cATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACat                                                                                                                                                    <  7:907172/150‑1 (MQ=255)
  cATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACat                                                                                                                                                    <  1:269789/150‑1 (MQ=255)
  cATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACat                                                                                                                                                    <  1:509857/150‑1 (MQ=255)
  cATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCCCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACata                                                                                                                                                   <  8:1051101/151‑1 (MQ=255)
   aTCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACatat                                                                                                                                                  <  8:349101/151‑1 (MQ=255)
   aTCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACatat                                                                                                                                                  >  6:305882/1‑151 (MQ=255)
   aTCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACatat                                                                                                                                                  >  5:945648/1‑151 (MQ=255)
      gCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATc                                                                                                                                                 >  5:44718/1‑149 (MQ=255)
      ccccAGCCAAGTAGCGCCATACAGACAATCATCGCCGTCAAAGTAATAAAACTCTGCTTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGCCCCATACCACTTCCAGCAGACCAGAAATAACTAAGATAATCCAGAACAAATCac                                                                                                                                               >  6:29019/2‑150 (MQ=255)
       cccAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAgg                                                                                                                                              <  6:574385/151‑1 (MQ=255)
          aGGCATGTAGCGCAATACTTACAATCATCGCCGTCAGAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTc                                                                                                                                           <  1:978568/151‑1 (MQ=255)
          aGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTc                                                                                                                                           <  5:707544/151‑1 (MQ=255)
           tgCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTcc                                                                                                                                          >  4:608107/2‑151 (MQ=255)
              aaGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCgg                                                                                                                                        >  3:154751/1‑150 (MQ=255)
              aaGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCgg                                                                                                                                        >  6:442470/1‑150 (MQ=255)
              aaGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCgg                                                                                                                                        >  1:374427/1‑150 (MQ=255)
               aGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCgg                                                                                                                                        >  8:586947/1‑149 (MQ=255)
               aGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGgag                                                                                                                                      <  1:942609/151‑1 (MQ=255)
                gTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGgaga                                                                                                                                     <  2:119299/151‑1 (MQ=255)
                gTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGgaga                                                                                                                                     <  3:566146/151‑1 (MQ=255)
                 tAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGgaga                                                                                                                                     >  8:828053/1‑150 (MQ=255)
                  aGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGaa                                                                                                                                    <  4:864303/150‑1 (MQ=255)
                       cATAGTGAGAATCATCGCCGTCACAGTAATTACACTCGTCGTCATAGGTCTATAGCCGTGGGTATATTTCAGGGCAACGGCGCATTCGACTTCGAGCAGACGAGCAATAAGTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTgggg                                                                                                                              <  8:219553/151‑1 (MQ=255)
                       cATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTgggg                                                                                                                              <  5:180509/151‑1 (MQ=255)
                          aCTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACATCCAGCAGACCAGCAATAACTAAGATAATCAAGAACATATCAGGCTCCGGAGAATTGGGGCCg                                                                                                                           >  6:823878/1‑151 (MQ=255)
                              aCAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTccc                                                                                                                       >  3:198096/1‑151 (MQ=255)
                               cAATCATCGCCGTCACAGTAATAACACTGGGCGTCAAACGACTAAAGCCGTTGGTATATTTCAGGCCAACCGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAAGATATCAGTCTCCGGAGAATTGGGGCCGTcccc                                                                                                                      <  2:747473/151‑1 (MQ=255)
                               cAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTcccc                                                                                                                      <  3:528106/151‑1 (MQ=255)
                                aaTCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCg                                                                                                                     >  3:896590/1‑151 (MQ=255)
                                  tcatcGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGt                                                                                                                   >  1:11809/1‑151 (MQ=255)
                                      cGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTTaa                                                                                                                <  8:489354/150‑1 (MQ=255)
                                         cGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGaaa                                                                                                               >  4:272534/1‑148 (MQ=255)
                                            cacaGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGc                                                                                                         <  4:385531/151‑1 (MQ=255)
                                             acaGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCg                                                                                                        >  7:746969/1‑151 (MQ=255)
                                                gTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGttt                                                                                                     <  7:617668/151‑1 (MQ=255)
                                                gTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGttt                                                                                                     <  3:299354/151‑1 (MQ=255)
                                                 taataaCACTCGGCGTCAAACGACTAAAGCCGTTGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCATCAATAACTAAGATAATCCAGAACATATCAGGCTACGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGttt                                                                                                     <  1:471424/150‑1 (MQ=255)
                                                 taataaCACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTg                                                                                                    >  8:539184/1‑151 (MQ=255)
                                                  aataaCACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGc                                                                                                   >  2:815001/1‑151 (MQ=255)
                                                    taaCACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGc                                                                                                   >  6:226360/1‑149 (MQ=255)
                                                     aaCACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAg                                                                                                 <  7:487260/150‑1 (MQ=255)
                                                       cacTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGgtg                                                                                              >  1:578749/1‑150 (MQ=255)
                                                       cacTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGgtc                                                                                              <  3:1029430/151‑1 (MQ=255)
                                                       cacTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGgtc                                                                                              <  3:420380/151‑1 (MQ=255)
                                                        acTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGgtcg                                                                                             <  1:787401/151‑1 (MQ=255)
                                                          tCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGgtcgtc                                                                                           >  7:289530/1‑151 (MQ=255)
                                                          tCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGgtcgtc                                                                                           >  5:522429/1‑151 (MQ=255)
                                                              cGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGc                                                                                       <  4:202792/151‑1 (MQ=255)
                                                                 cAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCaaa                                                                                    <  3:942484/151‑1 (MQ=255)
                                                                 cAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCaa                                                                                     >  6:1039863/1‑150 (MQ=255)
                                                                      gACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGa                                                                               <  2:432435/151‑1 (MQ=255)
                                                                       aCTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGAt                                                                              >  6:714431/1‑151 (MQ=255)
                                                                         tAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATgg                                                                            >  8:118118/1‑151 (MQ=255)
                                                                          aaaGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATgg                                                                            <  2:318593/150‑1 (MQ=255)
                                                                          aaaGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATg                                                                             <  4:450460/149‑1 (MQ=255)
                                                                          aaaGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGt                                                                           >  7:663849/1‑151 (MQ=255)
                                                                           aaGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGtt                                                                          >  1:714316/1‑151 (MQ=255)
                                                                             gCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGa                                                                        <  3:469311/151‑1 (MQ=255)
                                                                              ccGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGaa                                                                       <  2:38852/151‑1 (MQ=255)
                                                                              ccGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGaa                                                                       <  3:775700/151‑1 (MQ=255)
                                                                                gTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACg                                                                     <  2:27966/151‑1 (MQ=255)
                                                                                gTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACg                                                                     <  5:347049/151‑1 (MQ=255)
                                                                                 tGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACgg                                                                    >  4:761028/1‑151 (MQ=255)
                                                                                      atatTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGatat                                                                <  8:887693/150‑1 (MQ=255)
                                                                                      atatTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATAtt                                                               <  1:63355/151‑1 (MQ=255)
                                                                                       tatTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATAttt                                                              >  7:214858/1‑151 (MQ=255)
                                                                                         tttCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGCGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTg                                                            >  2:970792/1‑151 (MQ=255)
                                                                                           tCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGcc                                                          >  1:19380/1‑151 (MQ=255)
                                                                                             aGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAg                                                        <  3:352714/151‑1 (MQ=255)
                                                                                                 cAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTgg                                                    <  1:303252/151‑1 (MQ=255)
                                                                                                 cAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTgg                                                    <  8:776839/151‑1 (MQ=255)
                                                                                                            aCCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAAtt                                          <  4:324714/150‑1 (MQ=255)
                                                                                                             ccACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATGAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTc                                        >  5:829742/1‑151 (MQ=255)
                                                                                                              cACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCa                                       <  5:207513/151‑1 (MQ=255)
                                                                                                               aCTTCCAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCa                                       >  7:991291/1‑150 (MQ=255)
                                                                                                                   ccAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGGTGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTa                                   <  7:1071760/150‑1 (MQ=255)
                                                                                                                   ccAGCAGACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTAt                                  <  6:624448/151‑1 (MQ=255)
                                                                                                                        agACCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCt                             <  6:1006558/151‑1 (MQ=255)
                                                                                                                          aacatccaTAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGGCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCttt                           <  4:896018/144‑1 (MQ=255)
                                                                                                                          aCCAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCttt                           <  8:703745/151‑1 (MQ=255)
                                                                                                                           ccAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTtg                          <  5:665063/151‑1 (MQ=255)
                                                                                                                           ccAGCAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTtg                          >  6:976703/1‑151 (MQ=255)
                                                                                                                              ccAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGa                       <  3:729349/150‑1 (MQ=255)
                                                                                                                               cAATAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGAc                      >  2:878036/1‑151 (MQ=255)
                                                                                                                                aaTAACTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACt                     <  6:463127/151‑1 (MQ=255)
                                                                                                                                   aaCTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTAtg                  >  6:899787/1‑150 (MQ=255)
                                                                                                                                   aaCTAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATc                  >  3:662742/1‑151 (MQ=255)
                                                                                                                                      tAAGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGc                >  4:942802/1‑150 (MQ=255)
                                                                                                                                       aaGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGc                <  5:977827/149‑1 (MQ=255)
                                                                                                                                       aaGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCa               <  8:371769/150‑1 (MQ=255)
                                                                                                                                       aaGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAg              <  5:1040797/151‑1 (MQ=255)
                                                                                                                                        aGATAATCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGc                >  2:607859/1‑148 (MQ=255)
                                                                                                                                            aaTCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGaaaa          >  6:1079297/1‑150 (MQ=255)
                                                                                                                                             aTCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTGAACTATTACCTTTGTGACTATCGCAGaaaaa         >  6:382055/1‑150 (MQ=255)
                                                                                                                                             aTCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGaaaaaa        <  2:178512/151‑1 (MQ=255)
                                                                                                                                              tCCAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGAAAaaac       >  5:312570/1‑151 (MQ=255)
                                                                                                                                                cAGAACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGAAAaaaca      >  1:397449/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TTCATCGCCCAGGCAAGTAGCGCCATACTGACAATCATCGCCGTCACAGTAATAACACTCGGCGTCAAACGACTAAAGCCGTGGGTATATTTCAGGCCAACGGCCCATACCACTTCCAGCAGACCAGCAATAACTAAGATAATCCAGGACATATCAGGCTCCGGAGAATTGGGGCCGTCCCCGGTGAAAAGAGGCGTTTGCAGGTCGTCCTGCAAAGCTGATGGTTGAACGGATATTTTGCCAGGTGGCAGGGGAATTTCAACTATTACCTTTGTGACTATCGCAGAAAAAACAAACG  >  DH1/4137555‑4137852

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: