Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,874,049 G→C 5.1% C142W (TGC→TGG) ‡ dgoK ← 2‑dehydro‑3‑deoxygalactonokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,874,0490GC 5.1% 544.4 / 13.9 157C142W (TGC→TGG) ‡dgoK2‑dehydro‑3‑deoxygalactonokinase
Reads supporting (aligned to +/- strand):  ref base G (74/75);  new base C (4/4);  total (78/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 158 total aligned reads displayed.

AAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACGT  >  NC_000913/3873902‑3874198
                                                                                                                                                   |                                                                                                                                                     
aaGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCa                                                                                                                                                     <  2:2193226/150‑1 (MQ=255)
    ccAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATgg                                                                                                                                                 >  1:1089403/1‑150 (MQ=255)
     cAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATggg                                                                                                                                                <  1:2044714/150‑1 (MQ=255)
       gCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGtt                                                                                                                                              >  1:205044/1‑150 (MQ=255)
          gtgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCg                                                                                                                                           >  2:1975856/1‑150 (MQ=255)
           tgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgcagg                                                                                                                                          <  1:578215/150‑6 (MQ=255)
           tgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgcagg                                                                                                                                          <  1:2294108/150‑6 (MQ=255)
            gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCggg                                                                                                                                         <  1:1927922/150‑1 (MQ=255)
            gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCggg                                                                                                                                         <  2:2329866/150‑1 (MQ=255)
            gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgcaggc                                                                                                                                         <  2:852649/150‑7 (MQ=255)
             aaGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGGGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGc                                                                                                                                        <  2:570837/150‑1 (MQ=255)
                gCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATg                                                                                                                                     >  2:294978/1‑150 (MQ=255)
                 cATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGa                                                                                                                                    >  1:21824/1‑150 (MQ=255)
                 cATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGa                                                                                                                                    >  1:1984914/1‑150 (MQ=255)
                   tCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACa                                                                                                                                  <  2:1529562/150‑1 (MQ=255)
                   tCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACa                                                                                                                                  <  1:999489/150‑1 (MQ=255)
                    cGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACAt                                                                                                                                 >  1:1271174/1‑150 (MQ=255)
                     ggCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATa                                                                                                                                >  2:2047654/1‑150 (MQ=255)
                     ggCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATa                                                                                                                                >  1:1795884/1‑150 (MQ=255)
                      gCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCgg                                                                                                                                          >  2:1147974/1‑139 (MQ=255)
                       cAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATaaa                                                                                                                              <  1:1702407/150‑1 (MQ=255)
                         gagTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAaga                                                                                                                            <  1:208188/150‑1 (MQ=255)
                           gTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagaga                                                                                                                          >  2:1028513/1‑150 (MQ=255)
                            ttttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagagag                                                                                                                         <  1:838080/150‑1 (MQ=255)
                              ttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGagga                                                                                                                       <  1:1577146/150‑1 (MQ=255)
                              ttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGagga                                                                                                                       >  2:236144/1‑150 (MQ=255)
                                        ggcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGATGGAGCCaga                                                                                                             >  2:211530/1‑150 (MQ=255)
                                         aaaactcTGCGCCAATCAATGAGTGATTTAACAGTAAATGATATAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCagag                                                                                                            <  1:1226503/144‑1 (MQ=255)
                                          cAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGc                                                                                                           >  1:2340291/1‑150 (MQ=255)
                                             aCCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTcg                                                                                                        >  2:686421/1‑150 (MQ=255)
                                               cTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTcgcg                                                                                                      >  1:1556689/1‑150 (MQ=255)
                                                  cgcCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCCCCGGTCATCCCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCgcgcgc                                                                                                   >  1:2390242/1‑150 (MQ=255)
                                                   gcCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGcc                                                                                                  <  1:1421370/150‑1 (MQ=255)
                                                    ccAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCg                                                                                                 <  2:787735/150‑1 (MQ=255)
                                                    ccAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCg                                                                                                 <  1:409157/150‑1 (MQ=255)
                                                      aatcaatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCTATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGAt                                                                                               >  2:2224957/1‑150 (MQ=255)
                                                       atcaatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATc                                                                                              <  1:1015019/150‑1 (MQ=255)
                                                        tcaatGAGGGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCa                                                                                             <  2:813512/150‑1 (MQ=255)
                                                           atGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAAtt                                                                                          <  1:1542637/150‑1 (MQ=255)
                                                            tGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtg                                                                                         <  2:127679/150‑1 (MQ=255)
                                                            tGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtg                                                                                         <  1:308995/150‑1 (MQ=255)
                                                             gAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgt                                                                                        >  1:1988945/1‑150 (MQ=255)
                                                              aGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgtg                                                                                       <  2:1445687/150‑1 (MQ=255)
                                                                 gATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATc                                                                                              <  1:836048/140‑1 (MQ=255)
                                                                 gATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGttt                                                                                    >  2:1431441/1‑150 (MQ=255)
                                                                        cAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgcc                                                                             <  1:1488188/150‑1 (MQ=255)
                                                                        cAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgcc                                                                             >  1:404493/1‑150 (MQ=255)
                                                                           tAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCcgcg                                                                          <  1:1107109/150‑1 (MQ=255)
                                                                           tAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCcgcg                                                                          <  2:1102627/150‑1 (MQ=255)
                                                                              atgatgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTCCGCCGCGCAt                                                                       <  1:1882335/150‑1 (MQ=255)
                                                                               tgatgTAATTCACCGGGCATCACGGTGCGAAAATCTTTGATTTGCTGGCTATGGGCCTGCACCCATTTGCAATGGGTTGCGTGCATGACATAAAGAGAGGATGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTGTTCGCCGCGCATc                                                                      <  2:984175/150‑1 (MQ=255)
                                                                                 atgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATcac                                                                    >  1:910501/1‑150 (MQ=255)
                                                                                    tAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACAtt                                                                 >  2:270605/1‑150 (MQ=255)
                                                                                     aaTTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATtg                                                                >  1:1263915/1‑150 (MQ=255)
                                                                                      aTTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATtgt                                                               <  2:1785305/150‑1 (MQ=255)
                                                                                       ttCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATtgtg                                                              >  1:7496/1‑150 (MQ=255)
                                                                                        tCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGTAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTgg                                                             <  2:1555887/150‑1 (MQ=255)
                                                                                        tCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTgg                                                             <  1:753814/150‑1 (MQ=255)
                                                                                          aCCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGtt                                                           <  2:878061/150‑1 (MQ=255)
                                                                                           ccGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTa                                                          <  2:428536/150‑1 (MQ=255)
                                                                                           ccGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTa                                                          <  2:1566951/150‑1 (MQ=255)
                                                                                             ggTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATc                                                        <  1:1130205/150‑1 (MQ=255)
                                                                                                catcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTc                                                     <  2:171585/150‑1 (MQ=255)
                                                                                                  tcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCAt                                                   >  2:1764246/1‑150 (MQ=255)
                                                                                                   caCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATg                                                  <  2:1248041/150‑1 (MQ=255)
                                                                                                   caCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATg                                                  <  2:1089403/150‑1 (MQ=255)
                                                                                                     cGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgag                                                >  1:24701/1‑150 (MQ=255)
                                                                                                       gTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgagag                                              <  1:552760/150‑1 (MQ=255)
                                                                                                       gTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgagag                                              >  1:542772/1‑150 (MQ=255)
                                                                                                          cGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGaca                                           >  2:574442/1‑150 (MQ=255)
                                                                                                          cGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGaca                                           <  1:236144/150‑1 (MQ=255)
                                                                                                             aaaTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACAt                                        <  1:1856125/150‑1 (MQ=255)
                                                                                                             aaaTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACAt                                        <  2:2340291/150‑1 (MQ=255)
                                                                                                             aaaTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACAt                                        >  2:1193116/1‑150 (MQ=255)
                                                                                                                  gTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATcc                                   >  1:919051/1‑150 (MQ=255)
                                                                                                                    tGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCgg                                 <  1:1026533/150‑1 (MQ=255)
                                                                                                                    tGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCgg                                 <  1:1034666/150‑1 (MQ=255)
                                                                                                                    tGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCgg                                 <  2:228344/150‑1 (MQ=255)
                                                                                                                     taTTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCggg                                <  2:1600673/149‑1 (MQ=255)
                                                                                                                     taTTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCggg                                <  2:642739/149‑1 (MQ=255)
                                                                                                                     gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGCTCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCggg                                >  1:2143669/1‑150 (MQ=255)
                                                                                                                     gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCggg                                >  2:738245/1‑150 (MQ=255)
                                                                                                                               cTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAg                      >  2:2175646/1‑150 (MQ=255)
                                                                                                                               cTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAg                      <  2:496044/150‑1 (MQ=255)
                                                                                                                                tATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGa                     >  1:1048594/1‑150 (MQ=255)
                                                                                                                                 aTCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGat                    <  1:1343773/150‑1 (MQ=255)
                                                                                                                                  tCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGCAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGata                   <  2:1282333/150‑1 (MQ=255)
                                                                                                                                    ggCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATAtt                 >  1:1989543/1‑150 (MQ=255)
                                                                                                                                     gCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTg                <  1:2041629/150‑1 (MQ=255)
                                                                                                                                     cccggaaCCCATTGCAAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTg                >  1:1160903/7‑150 (MQ=255)
                                                                                                                                      ccggaaCCCATTGCAAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGt               >  1:1771871/6‑150 (MQ=255)
                                                                                                                                      ccTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGt               >  2:2285923/1‑150 (MQ=255)
                                                                                                                                        tGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCg             >  2:1029642/1‑150 (MQ=255)
                                                                                                                                         gaaCCCATTGCAAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGc            >  1:494412/3‑150 (MQ=255)
                                                                                                                                         gCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGc            <  2:1988945/150‑1 (MQ=255)
                                                                                                                                         gCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGc            >  2:2005753/1‑150 (MQ=255)
                                                                                                                                            cccATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTAGCCCCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCaa         >  2:1397511/1‑150 (MQ=255)
                                                                                                                                              cATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACt       >  2:1578894/1‑150 (MQ=255)
                                                                                                                                              cATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACt       >  1:475560/1‑150 (MQ=255)
                                                                                                                                                   gCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACGt  >  2:453420/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACGT  >  NC_000913/3873902‑3874198

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