Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,692,221 T→A 5.3% E26V (GAA→GTA)  hdhA ← 7‑alpha‑hydroxysteroid dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,692,2210TA 5.3% 973.5 / 32.0 290E26V (GAA→GTA) hdhA7‑alpha‑hydroxysteroid dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (153/121);  new base A (9/7);  total (162/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

AGAGTTCCTGTTCGGAAGTAATATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAACCTTG  >  CP009273/1692072‑1692369
                                                                                                                                                     |                                                                                                                                                    
 gagTTCCTGTTCGGAAGTAATATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTc                                                                                                                                                     >  2:551194/1‑150 (MQ=255)
  agTTCCTGTTCGGAAGTAATATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCt                                                                                                                                                    <  1:4089802/150‑1 (MQ=255)
    ttCCTGTTCGGAAGTAATATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCttt                                                                                                                                                  >  2:1252171/1‑150 (MQ=255)
     tCCTGTTCGGAAGTAATATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTa                                                                                                                                                 >  1:4267326/1‑150 (MQ=255)
      ccTGTTCGGAAGTAATATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTAc                                                                                                                                                <  1:4564404/150‑1 (MQ=255)
                aGTAATATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGc                                                                                                                                      >  2:2195478/1‑150 (MQ=255)
                 gTAATATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCa                                                                                                                                     <  2:857447/150‑1 (MQ=255)
                   aaTATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCAcc                                                                                                                                   <  1:3910882/150‑1 (MQ=255)
                   aaTATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCAcc                                                                                                                                   <  1:2403391/150‑1 (MQ=255)
                    atatCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCAccc                                                                                                                                  <  1:4585760/150‑1 (MQ=255)
                      atCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGc                                                                                                                                >  1:1061961/1‑150 (MQ=255)
                      atCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGc                                                                                                                                >  2:1459148/1‑150 (MQ=255)
                        cacaACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCAc                                                                                                                              >  2:2288521/1‑150 (MQ=255)
                        cacaACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCAc                                                                                                                              >  2:268682/1‑150 (MQ=255)
                          caACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCt                                                                                                                            <  1:2803045/150‑1 (MQ=255)
                            aCGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCtgt                                                                                                                          >  2:2619495/1‑150 (MQ=255)
                             cGGCAGGCATTTGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATTGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTTCCAATACCTGCACCCGCACCtgtg                                                                                                                         <  2:654667/150‑1 (MQ=255)
                              ggcaggcaAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGtga                                                                                                                        <  1:3720129/150‑1 (MQ=255)
                              ggcaggcaAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGtga                                                                                                                        <  1:3340651/150‑1 (MQ=255)
                                caggcaAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGtgatg                                                                                                                      <  1:1246353/150‑1 (MQ=255)
                                caggcaAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGtgatg                                                                                                                      >  1:2066559/1‑150 (MQ=255)
                                 aggcaAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGtgatga                                                                                                                     >  2:3024042/1‑150 (MQ=255)
                                 aggcaAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGtgatga                                                                                                                     <  1:2202505/150‑1 (MQ=255)
                                     aaaTGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGgc                                                                                                                 <  1:3200946/150‑1 (MQ=255)
                                        tGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCa                                                                                                              <  1:1420354/150‑1 (MQ=255)
                                         gCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCAt                                                                                                             >  1:4045567/1‑150 (MQ=255)
                                             gACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTACTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTc                                                                                                         >  2:778472/1‑150 (MQ=255)
                                                 accCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTc                                                                                                     <  1:2955424/150‑1 (MQ=255)
                                                   ccAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGa                                                                                                   <  1:3048276/150‑1 (MQ=255)
                                                   ccAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGa                                                                                                   >  1:3154167/1‑150 (MQ=255)
                                                    cAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAg                                                                                                  <  1:3768676/150‑1 (MQ=255)
                                                     aGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGt                                                                                                 >  1:1009470/1‑150 (MQ=255)
                                                      gttgttGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGtc                                                                                                <  2:2126219/150‑1 (MQ=255)
                                                          ttgAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAg                                                                                            <  1:2108545/150‑1 (MQ=255)
                                                          ttgAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAg                                                                                            <  2:2462941/150‑1 (MQ=255)
                                                             aaTTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGtt                                                                                         >  1:4589065/1‑150 (MQ=255)
                                                               tttCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGt                                                                                       >  2:1783547/1‑150 (MQ=255)
                                                                ttCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTc                                                                                      >  2:2956497/1‑150 (MQ=255)
                                                                   gTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGa                                                                                   >  2:1415456/1‑150 (MQ=255)
                                                                     cTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAAt                                                                                 <  1:159638/150‑1 (MQ=255)
                                                                       acaacaTGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTa                                                                               <  1:2652542/150‑1 (MQ=255)
                                                                        caacaTGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTaa                                                                              >  1:1137258/1‑150 (MQ=255)
                                                                         aacaTGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTaaa                                                                             >  2:1489031/1‑150 (MQ=255)
                                                                           caTGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAaca                                                                           <  2:1592391/150‑1 (MQ=255)
                                                                           caTGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTATTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTACGTCGAGTCTCAGGTTTTCAGAATTAAaca                                                                           >  1:817193/1‑150 (MQ=255)
                                                                            aTGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGGGCCAGCTGTCGCGAATGTAATGGCGATTACTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAacac                                                                          <  1:765111/150‑1 (MQ=255)
                                                                            aTGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAacac                                                                          <  2:1562286/150‑1 (MQ=255)
                                                                              ggTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACAt                                                                        >  2:2328263/1‑150 (MQ=255)
                                                                              ggTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACAt                                                                        <  2:476146/150‑1 (MQ=255)
                                                                               gTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATa                                                                       >  2:1879461/1‑150 (MQ=255)
                                                                               gTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATa                                                                       <  2:88120/150‑1 (MQ=255)
                                                                                ttAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATaa                                                                      >  2:4022165/1‑150 (MQ=255)
                                                                                     tGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAAcctc                                                                 <  1:987486/150‑1 (MQ=255)
                                                                                     tGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTACTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAAcctc                                                                 <  1:1691164/150‑1 (MQ=255)
                                                                                       cGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCATTACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATCAAcctcct                                                               <  2:4476624/150‑1 (MQ=255)
                                                                                       cGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAAcctcct                                                               <  2:3105383/150‑1 (MQ=255)
                                                                                         tCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCtgt                                                             <  2:3419616/150‑1 (MQ=255)
                                                                                           ggCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTACTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGa                                                           >  2:954957/1‑150 (MQ=255)
                                                                                           ggCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTACTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGa                                                           >  2:3864762/1‑150 (MQ=255)
                                                                                             cGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAAc                                                         <  2:11910/150‑1 (MQ=255)
                                                                                             cGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAAc                                                         >  1:3659494/1‑150 (MQ=255)
                                                                                             cGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAAc                                                         >  2:1186795/1‑150 (MQ=255)
                                                                                              gTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACt                                                        <  1:1485797/150‑1 (MQ=255)
                                                                                              gTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACt                                                        >  2:3717789/1‑150 (MQ=255)
                                                                                                tAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTc                                                      <  1:663783/150‑1 (MQ=255)
                                                                                                 aaTATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCa                                                     >  1:834285/1‑150 (MQ=255)
                                                                                                  atatCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAg                                                    >  1:2466448/1‑150 (MQ=255)
                                                                                                    atCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTa                                                  <  1:1953018/150‑1 (MQ=255)
                                                                                                         tGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAAt                                             <  1:3412767/150‑1 (MQ=255)
                                                                                                            ccaccaCAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATaaa                                          >  2:1634054/1‑150 (MQ=255)
                                                                                                              accacAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAg                                        >  2:4608548/1‑150 (MQ=255)
                                                                                                                 acaGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTTTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCa                                     <  2:3280719/150‑1 (MQ=255)
                                                                                                                 acaGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCa                                     <  1:523331/150‑1 (MQ=255)
                                                                                                                  caGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCAc                                    >  2:134564/1‑150 (MQ=255)
                                                                                                                    gATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAGAAATAAAAGGCACCg                                  <  1:605735/150‑1 (MQ=255)
                                                                                                                     aTGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGc                                 <  1:2909373/150‑1 (MQ=255)
                                                                                                                        cgcCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCtg                              >  1:3279295/1‑150 (MQ=255)
                                                                                                                        cgcCAGCTGTCGCGAATGTAATGGCGATTTATTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCtg                              >  1:1372314/1‑150 (MQ=255)
                                                                                                                           cAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCtgtgt                           <  1:3562320/150‑1 (MQ=255)
                                                                                                                              cTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGcc                        >  2:2726340/1‑150 (MQ=255)
                                                                                                                              cTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGcc                        >  2:1099257/1‑150 (MQ=255)
                                                                                                                               tGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCca                       >  1:3988254/1‑150 (MQ=255)
                                                                                                                                      aaTGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACAttttt                <  1:3185262/150‑1 (MQ=255)
                                                                                                                                      aaTGTAACGGCGCTTTCATCACCAATACCAGCACCCGCACCCGTGATGACGGCGCACTTTCCGTCGAGCCTCAGGTTGTCAGAAATAAACCCATAAACCTCCTGTTAACTTCAGTACAAATACTAGCCATCGCCTGTGTGCCACAttttt                <  2:2913104/150‑1 (MQ=255)
                                                                                                                                         gTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTAt             >  2:2764387/1‑150 (MQ=255)
                                                                                                                                            aTGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCAt          >  2:1525623/1‑150 (MQ=255)
                                                                                                                                             tGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCAtt         >  1:4398023/1‑150 (MQ=255)
                                                                                                                                             tGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCAtt         >  2:4580533/1‑150 (MQ=255)
                                                                                                                                              ggCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTa        <  1:890292/150‑1 (MQ=255)
                                                                                                                                                cGATTACTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAAc      >  1:4187821/1‑150 (MQ=255)
                                                                                                                                                cGATTACTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAAc      >  1:3681937/1‑150 (MQ=255)
                                                                                                                                                 gATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAAcc     >  2:4303390/1‑150 (MQ=255)
                                                                                                                                                 gATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAAcc     >  2:4165493/1‑150 (MQ=255)
                                                                                                                                                  aTTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAACCt    <  1:2082813/150‑1 (MQ=255)
                                                                                                                                                  aTTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAACCt    >  1:174518/1‑150 (MQ=255)
                                                                                                                                                  aTTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAACCt    <  2:4212373/150‑1 (MQ=255)
                                                                                                                                                  aTTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAACCt    <  2:2150422/150‑1 (MQ=255)
                                                                                                                                                  aTTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAACCt    >  2:908751/1‑150 (MQ=255)
                                                                                                                                                    ttCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAACCTTg  <  2:2319992/150‑1 (MQ=255)
                                                                                                                                                    gtCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAACCTTg  <  2:1074032/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGAGTTCCTGTTCGGAAGTAATATCACAACGGCAGGCAAATGCCTGACCACCCAGTTGTTGAATTTCGTCTACAACATGGTTAGCTGCGTCGGCGTTAATATCACTGACCACCACAGATGCGCCAGCTGTCGCGAATGTAATGGCGATTTCTTTACCAATACCTGCACCCGCACCTGTGATGATGGCGCATTTTCCGTCGAGTCTCAGGTTGTCAGAATTAAACACATAAACCTCCTGTGAACTTCAGTAAAAATAAAAGGCACCGCCTGTGTGCCACATTTTTTATCATTAACCTTG  >  CP009273/1692072‑1692369

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