Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,173,457 T→G W169G (TGG→GGG)  melR → melibiose operon regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,173,4570TG99.1% 2209.0 / ‑6.3 580W169G (TGG→GGG) melRmelibiose operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base G (294/281);  total (295/285)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 608 total aligned reads displayed.

CATGGTGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCTGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGCACGTCA  >  DH1/4173309‑4173604
                                                                                                                                                    |                                                                                                                                                   
cATGGTGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTctggcggg                                                                                                                                                   >  2:1789500/1‑151 (MQ=255)
cATGGTGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTctggcggg                                                                                                                                                   >  3:935878/1‑151 (MQ=255)
cATGGTGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTctggcggg                                                                                                                                                   >  5:1471036/1‑151 (MQ=255)
     tGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAAcc                                                                                                                                              >  1:1198443/1‑151 (MQ=255)
     tGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAAcc                                                                                                                                              >  2:1610571/1‑151 (MQ=255)
      gATCAAATCACTGGCCACACCGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAAcc                                                                                                                                              <  2:1499499/150‑1 (MQ=255)
      gATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCg                                                                                                                                             <  6:778968/151‑1 (MQ=255)
      gATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCg                                                                                                                                             <  8:372336/151‑1 (MQ=255)
       aTCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGa                                                                                                                                            <  1:1314738/151‑1 (MQ=255)
         cAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGAt                                                                                                                                           <  1:1047536/150‑1 (MQ=255)
           aaTCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCt                                                                                                                                        >  2:306305/1‑151 (MQ=255)
            aTCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTg                                                                                                                                       <  1:242086/151‑1 (MQ=255)
            aTCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGAGTCTg                                                                                                                                       >  6:1077313/1‑151 (MQ=255)
               aCTGGCCACACAGCAACTTAGCCCGGTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTc                                                                                                                                    <  6:1923629/151‑1 (MQ=255)
                 tGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCa                                                                                                                                   >  3:347344/1‑150 (MQ=255)
                   gCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAAt                                                                                                                                 >  1:1883362/1‑150 (MQ=255)
                   gCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATa                                                                                                                                <  3:456639/151‑1 (MQ=255)
                    ccACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATaa                                                                                                                               <  3:1685427/151‑1 (MQ=255)
                     cacacaGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATaa                                                                                                                               >  3:1205286/1‑150 (MQ=255)
                       cacaGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAAc                                                                                                                            <  2:1096757/151‑1 (MQ=255)
                         caGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCt                                                                                                                          >  7:915385/1‑151 (MQ=255)
                          aGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTc                                                                                                                         <  4:823420/151‑1 (MQ=255)
                                 tAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGAGcac                                                                                                                   >  8:772949/1‑150 (MQ=255)
                                   gCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGcacaca                                                                                                                <  2:7338/151‑1 (MQ=255)
                                      cGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGcacacaca                                                                                                              >  3:1633814/1‑150 (MQ=255)
                                      cGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACaa                                                                                                             <  2:905150/151‑1 (MQ=255)
                                       gTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAACCCTCGCGCACACACaaa                                                                                                            <  4:123565/151‑1 (MQ=255)
                                         ttGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACaaaa                                                                                                           >  1:543418/1‑150 (MQ=255)
                                         ttGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACaaaa                                                                                                           <  1:1787326/150‑1 (MQ=255)
                                          tGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACaaaaa                                                                                                          <  1:964235/150‑1 (MQ=255)
                                            aaGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAg                                                                                                       <  8:649571/151‑1 (MQ=255)
                                             aGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGc                                                                                                      >  1:838180/1‑151 (MQ=255)
                                             aGTCCGCCGCTGGCAGCAGGAATTAAACAGTCCGAAAGATCAAATTCGCCAGATCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGc                                                                                                      <  8:1448059/151‑1 (MQ=255)
                                               tGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGt                                                                                                    >  7:809430/1‑151 (MQ=255)
                                                gcgcCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGtc                                                                                                   <  6:644903/151‑1 (MQ=255)
                                                gcgcCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGt                                                                                                    <  5:1825158/150‑1 (MQ=255)
                                                gcgcCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGt                                                                                                    <  8:2683/150‑1 (MQ=255)
                                                gcgcCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGt                                                                                                    <  5:435127/150‑1 (MQ=255)
                                                   ccgcTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGCGAACCGATTCTGGTCAATAACACCTCGCGCACACACAAAAACAGCGTCTcg                                                                                                >  8:50111/1‑151 (MQ=255)
                                                     gcTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGAGGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTcgcg                                                                                              >  4:1156003/1‑151 (MQ=255)
                                                       tGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGcc                                                                                            >  7:3928/1‑151 (MQ=255)
                                                          cagcagGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATg                                                                                         >  6:1790914/1‑151 (MQ=255)
                                                           agcagGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATgc                                                                                        >  4:1459571/1‑151 (MQ=255)
                                                             cagGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATgcgc                                                                                      >  6:1640469/1‑151 (MQ=255)
                                                              agGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAACCAGCGTCTCGCGCCATgcgc                                                                                      >  6:1416147/1‑150 (MQ=255)
                                                              agGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCa                                                                                     <  1:1555327/151‑1 (MQ=255)
                                                                 aaTTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAAtt                                                                                  <  5:839654/151‑1 (MQ=255)
                                                                 aaTTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAAtt                                                                                  >  8:181935/1‑151 (MQ=255)
                                                                  aTTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAAttt                                                                                 <  5:284673/151‑1 (MQ=255)
                                                                   ttAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAAtttt                                                                                >  3:821091/1‑151 (MQ=255)
                                                                    tAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGACTCGCGCCATGCGCAAtttt                                                                                >  2:1278887/1‑150 (MQ=255)
                                                                     aaaCAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTAt                                                                              >  7:623398/1‑151 (MQ=255)
                                                                       aCAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGt                                                                            >  6:1487019/1‑151 (MQ=255)
                                                                       aCAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGt                                                                            >  3:1434550/1‑151 (MQ=255)
                                                                        cAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGtt                                                                           >  8:1952817/1‑151 (MQ=255)
                                                                        cAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGtt                                                                           <  3:579830/151‑1 (MQ=255)
                                                                            ccGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGcc                                                                       >  5:1093911/1‑151 (MQ=255)
                                                                            ccGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGcc                                                                       <  4:1392047/151‑1 (MQ=255)
                                                                                   aGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCt                                                                >  7:1332099/1‑151 (MQ=255)
                                                                                   aGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTAGCGCCATGCGCAATTTTATGTTAGCCAGATGCt                                                                >  4:275267/1‑151 (MQ=255)
                                                                                     cAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTgg                                                              <  4:927239/151‑1 (MQ=255)
                                                                                     cAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTgg                                                              <  7:1656513/151‑1 (MQ=255)
                                                                                     cAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTgg                                                              >  8:109956/1‑151 (MQ=255)
                                                                                      aaaTTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTggg                                                             <  5:68166/151‑1 (MQ=255)
                                                                                         ttCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCtt                                                          >  4:1050438/1‑151 (MQ=255)
                                                                                          tCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCtt                                                          <  7:1462189/150‑1 (MQ=255)
                                                                                           cGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCGAGATGCTGGGCttt                                                         <  4:713466/150‑1 (MQ=255)
                                                                                                gCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGcc                                                   >  1:1392985/1‑151 (MQ=255)
                                                                                                gCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGcc                                                   <  5:1262905/151‑1 (MQ=255)
                                                                                                  tCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCg                                                  >  1:362117/1‑150 (MQ=255)
                                                                                                     ccATTGATGAAATTGGCCTGATTCTCAAGCGATTTAGCCTCTCTGGCTGGGAACCGATTCTTGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCTCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGaaaa                                              <  8:1669812/151‑1 (MQ=255)
                                                                                                       aTTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACt                                            >  7:317864/1‑151 (MQ=255)
                                                                                                       aTTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGGCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACt                                            <  7:259676/151‑1 (MQ=255)
                                                                                                       aTTGATGAAATTGGCCTGATGCTAAAGCGATTTAGCCTCTCTGCCGGCGAACCGATTCTGGTCAATTAAACCTCGCGCACACAGAAAAACAGCCTCTCGAGCAATGCGCAATTTTATCTTATCAAGATGGTGGGCTTTATTGCCGAAAACt                                            >  8:1556456/1‑151 (MQ=255)
                                                                                                         tgatgaAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTa                                           <  5:594183/150‑1 (MQ=255)
                                                                                                         tgatgaAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTAt                                          >  3:721186/1‑151 (MQ=255)
                                                                                                          gatgaAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTAt                                          <  6:522004/150‑1 (MQ=255)
                                                                                                           atgaAATTGGCCTGTTGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGa                                        <  4:181064/151‑1 (MQ=255)
                                                                                                           atgaAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATg                                         >  7:1594695/1‑150 (MQ=255)
                                                                                                             gaAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGAt                                       <  4:488886/150‑1 (MQ=255)
                                                                                                             gaAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATc                                      <  3:1006993/151‑1 (MQ=255)
                                                                                                             gaAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATc                                      <  6:897293/151‑1 (MQ=255)
                                                                                                                   ggCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGgcg                                 <  2:1929101/150‑1 (MQ=255)
                                                                                                                    gCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCt                               <  3:1113051/151‑1 (MQ=255)
                                                                                                                     ccTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTg                              <  3:130725/151‑1 (MQ=255)
                                                                                                                      cTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTg                              >  7:1812804/1‑150 (MQ=255)
                                                                                                                       tGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGTGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGAc                            >  8:1283435/1‑151 (MQ=255)
                                                                                                                        gATGCTCAAGCTATTTAGGCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAAGAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGGTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGTcc                           <  4:1643270/151‑1 (MQ=255)
                                                                                                                           gCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGGCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATc                        >  7:1175607/1‑151 (MQ=255)
                                                                                                                             tCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCa                       >  4:1930523/1‑150 (MQ=255)
                                                                                                                              cAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAAc                     >  1:1711095/1‑151 (MQ=255)
                                                                                                                                        aGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCt            <  4:349513/150‑1 (MQ=255)
                                                                                                                                        aGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCACCAACGATGTGGCt            <  6:1238761/150‑1 (MQ=255)
                                                                                                                                          ccTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAg         <  1:620255/151‑1 (MQ=255)
                                                                                                                                           ctctctGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGc        >  6:371403/1‑151 (MQ=255)
                                                                                                                                           ctctctGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGc        >  3:1829781/1‑151 (MQ=255)
                                                                                                                                            tctctGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGc        >  5:1617368/1‑150 (MQ=255)
                                                                                                                                              tctGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCt            >  3:1933031/1‑144 (MQ=255)
                                                                                                                                                 ggcgggGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGCACGt    <  2:462812/149‑1 (MQ=255)
                                                                                                                                                 ggcgggGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAACCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGCACg     <  2:547080/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATGGTGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCTGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGCACGTCA  >  DH1/4173309‑4173604

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: