Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,874,048 T→G 6.1% K143Q (AAA→CAA)  dgoK ← 2‑oxo‑3‑deoxygalactonate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,874,0480TG 6.1% 408.6 / 12.7 141K143Q (AAA→CAA) dgoK2‑oxo‑3‑deoxygalactonate kinase
Reads supporting (aligned to +/- strand):  ref base T (59/73);  new base G (3/6);  total (62/79)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

Only 100 of 145 total aligned reads displayed.

CTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACG  >  NC_000913/3873899‑3874197
                                                                                                                                                     |                                                                                                                                                     
  cAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGc                                                                                                                                                     <  2:997286/150‑1 (MQ=255)
  cAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGc                                                                                                                                                     <  2:1436858/150‑1 (MQ=255)
    aGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCaa                                                                                                                                                   <  1:1062965/150‑1 (MQ=255)
    aGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCaa                                                                                                                                                   >  1:400058/1‑150 (MQ=255)
      gCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATg                                                                                                                                                 >  1:761428/1‑150 (MQ=255)
       cccgCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCTATCAATGAGTGATTTCACAGTAAATGATGTTATTCACCGGTGATCACGGTGCGAAAATCGTTGATTTGGTGGCTATCGGCGTGCACCCATTTTCAATgg                                                                                                                                                <  2:1787176/147‑1 (MQ=255)
        cAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATggg                                                                                                                                               >  2:503811/1‑150 (MQ=255)
        cAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATggg                                                                                                                                               >  2:661448/1‑150 (MQ=255)
        cAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATggg                                                                                                                                               <  1:1378230/150‑1 (MQ=255)
        cAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATATAACAGTAAATGATGAAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCGGCACCCATTGCAAATggg                                                                                                                                               <  1:1414047/150‑1 (MQ=255)
         aGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGt                                                                                                                                              <  1:33578/150‑1 (MQ=255)
            tgtgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTcc                                                                                                                                           >  2:1113609/1‑150 (MQ=255)
             gtgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCg                                                                                                                                          >  1:1557710/1‑150 (MQ=255)
              tgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCgg                                                                                                                                         >  2:204911/1‑150 (MQ=255)
              tgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCgg                                                                                                                                         <  2:1645007/150‑1 (MQ=255)
               gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACGGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCggg                                                                                                                                        <  1:201956/150‑1 (MQ=255)
               gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCggg                                                                                                                                        >  2:1368358/1‑150 (MQ=255)
               gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCggg                                                                                                                                        >  1:1351404/1‑150 (MQ=255)
               gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGAATTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgcaggc                                                                                                                                        <  1:1448128/150‑7 (MQ=255)
                aaGGCTTCTGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGc                                                                                                                                       <  1:1639556/150‑1 (MQ=255)
                aaGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGc                                                                                                                                       >  1:117015/1‑150 (MQ=255)
                 aGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCa                                                                                                                                      <  2:836730/150‑1 (MQ=255)
                 aGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCa                                                                                                                                      >  1:1303297/1‑150 (MQ=255)
                    cATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGa                                                                                                                                   <  1:1190688/150‑1 (MQ=255)
                    cATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGa                                                                                                                                   >  1:1126840/1‑150 (MQ=255)
                     aTCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGAc                                                                                                                                  <  1:325651/150‑1 (MQ=255)
                           agagTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAag                                                                                                                            <  1:1320447/150‑1 (MQ=255)
                            gagTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAaga                                                                                                                           <  2:1477233/150‑1 (MQ=255)
                            gagTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAaga                                                                                                                           >  1:233137/1‑150 (MQ=255)
                               ttttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagagag                                                                                                                        <  2:1873487/150‑1 (MQ=255)
                               ttttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagagag                                                                                                                        <  2:932612/150‑1 (MQ=255)
                                tttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGagg                                                                                                                       <  2:926556/150‑1 (MQ=255)
                                 ttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACTGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGagga                                                                                                                      <  2:536703/150‑1 (MQ=255)
                                   ccTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaag                                                                                                                    >  1:625326/1‑150 (MQ=255)
                                   ccTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaag                                                                                                                    >  2:1013736/1‑150 (MQ=255)
                                       cggcggCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGc                                                                                                                >  1:1612937/1‑150 (MQ=255)
                                          cggcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCag                                                                                                             >  1:675885/1‑150 (MQ=255)
                                            gcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCagag                                                                                                           >  2:554820/1‑150 (MQ=255)
                                             cAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGc                                                                                                          >  2:613497/1‑150 (MQ=255)
                                               aaCCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTc                                                                                                        >  2:110569/1‑150 (MQ=255)
                                                   tGCGCCAATCAATGAGTGATTTAACAGTAATTGATTTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTAGCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTcgcgc                                                                                                    <  1:158581/150‑1 (MQ=255)
                                                      gcCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGcc                                                                                                 >  1:1392945/1‑150 (MQ=255)
                                                       ccAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCg                                                                                                <  1:1114010/150‑1 (MQ=255)
                                                       ccAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCg                                                                                                >  1:1832414/1‑150 (MQ=255)
                                                          atcaatGAGTGATTTAACAGTAAATGATGTAATTCCCCGGTCATCACGGTGCGAAACGCGTTGAGTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGCGAGGAAGGAGCCAGAGCTCGCGCGCCGATc                                                                                             <  1:1704214/150‑1 (MQ=255)
                                                          atcaatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATc                                                                                             >  1:1867948/1‑150 (MQ=255)
                                                            caatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCaa                                                                                           >  1:842043/1‑150 (MQ=255)
                                                             aatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAAt                                                                                          >  2:980660/1‑150 (MQ=255)
                                                                 aGTGTTTTAACAGTAAATGAGGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATTGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgtg                                                                                      <  1:770686/150‑1 (MQ=255)
                                                                 aGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgtg                                                                                      <  2:1685358/150‑1 (MQ=255)
                                                                  gTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTgtgt                                                                                     >  1:1287576/1‑150 (MQ=255)
                                                                   tGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTCTCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGtt                                                                                    <  2:1869982/150‑1 (MQ=255)
                                                                    gATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGTCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGttt                                                                                   >  1:1436896/1‑150 (MQ=255)
                                                                    gATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGttt                                                                                   >  1:1855676/1‑150 (MQ=255)
                                                                        tAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCCCCCATTTGCAATGGGTACCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGAACAAATGTGTttcgtc                                                                               >  1:287201/1‑150 (MQ=255)
                                                                            tgTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGAGATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGAGCAATTGTGTTTCTTcgccg                                                                           >  2:185671/2‑150 (MQ=255)
                                                                                aaTGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCgcgcc                                                                       >  1:1876078/1‑149 (MQ=255)
                                                                                aaTGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCa                                                                       <  1:974626/150‑1 (MQ=255)
                                                                                 atgatgTAATTCACCGGTCATCACGGTGCGAAAATCGTGGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCAt                                                                      <  2:1587296/150‑1 (MQ=255)
                                                                                  ttatgTAATTCACCGGTCATCACGTTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGAATTACATAAAGAGACGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATc                                                                     <  2:1623705/148‑1 (MQ=255)
                                                                                            cACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGt                                                           <  1:1553949/150‑1 (MQ=255)
                                                                                                 gTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCg                                                      <  2:1512728/150‑1 (MQ=255)
                                                                                                 gTCAGCACGGTGCGAAAATCGTTGATTTGCTGGCTTTCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCg                                                      <  1:334289/150‑1 (MQ=255)
                                                                                                  tcatcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGt                                                     <  1:1630605/150‑1 (MQ=255)
                                                                                                  tcatcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGt                                                     <  1:1555293/150‑1 (MQ=255)
                                                                                                    atcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGTAGCCAGAGCTCGCGCGCAGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCa                                                   <  1:126657/150‑1 (MQ=255)
                                                                                                    atcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCa                                                   <  2:1044208/150‑1 (MQ=255)
                                                                                                     tcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCAt                                                  >  1:538752/1‑150 (MQ=255)
                                                                                                      caCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATg                                                 >  1:616069/1‑150 (MQ=255)
                                                                                                      caCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATg                                                 <  2:568978/150‑1 (MQ=255)
                                                                                                      caCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATg                                                 <  2:1832414/150‑1 (MQ=255)
                                                                                                      caCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATg                                                 <  2:117015/150‑1 (MQ=255)
                                                                                                       aCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATga                                                <  1:164148/150‑1 (MQ=255)
                                                                                                           tGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACACATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgagaga                                            >  2:545416/1‑150 (MQ=255)
                                                                                                             cGAAAATCGGTGATTTGCTCGGTTTCGGCCTGCACGCTTTTGCCTTGGGTTCCGGGCATGACATAAAGACCGGAAGGAGCCAGAGCTCGCGCGCCGATCATTTGTGTTGCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGaca                                          <  2:364495/150‑1 (MQ=255)
                                                                                                               aaaaTCTTTGATTTGCTGGCTATCGGACTGCACCAATTTGCAATGGGTTCCGGGCATGACATATAGTGCGGAAGGAGCAAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCTTGAGAGAcaca                                        <  1:1147206/150‑1 (MQ=255)
                                                                                                                aaaTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACAt                                       >  2:1381064/1‑150 (MQ=255)
                                                                                                                       tGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCgg                                >  1:749047/1‑150 (MQ=255)
                                                                                                                       tGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCgg                                <  2:233137/150‑1 (MQ=255)
                                                                                                                        gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCggg                               <  2:400058/150‑1 (MQ=255)
                                                                                                                              cTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGTGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGaataat                         <  1:4160/150‑1 (MQ=255)
                                                                                                                               tGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATc                        <  1:653012/150‑1 (MQ=255)
                                                                                                                               tGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATc                        <  2:1351404/150‑1 (MQ=255)
                                                                                                                                 gCTATCGGCATGCACCCATTTGCAATGGGTTCGGGGCATGACATAAAGAGAGGATGGAGACAGAGGTCGCGAGCCGATCAATTGTGTTTCTTCGCCGCGCATCTCATTGTGGTTTTCTTCATGAGAGACTCATAATCCGGGAATAATCCa                      <  2:233772/150‑1 (MQ=255)
                                                                                                                                  cTATCGGCCTGCACGCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAg                     <  2:1576002/150‑1 (MQ=255)
                                                                                                                                  cTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCCCTTTGTGGTTATCGTCATGAGAGACACATCATCCGGGAATAATCCAg                     >  1:140748/1‑150 (MQ=255)
                                                                                                                                  cTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAg                     <  2:538752/150‑1 (MQ=255)
                                                                                                                                  cTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAg                     <  1:661448/150‑1 (MQ=255)
                                                                                                                                   tATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGa                    <  1:110569/150‑1 (MQ=255)
                                                                                                                                       ggCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGGCGCTCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATAtt                >  2:298646/1‑150 (MQ=255)
                                                                                                                                       ggCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATAtt                >  2:1198368/1‑150 (MQ=255)
                                                                                                                                           tGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCg            <  2:1855676/150‑1 (MQ=255)
                                                                                                                                           tGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCg            <  1:223302/150‑1 (MQ=255)
                                                                                                                                            gCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGc           <  1:1351846/150‑1 (MQ=255)
                                                                                                                                            gCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGc           >  2:834270/1‑150 (MQ=255)
                                                                                                                                             cACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGcc          >  1:602846/1‑150 (MQ=255)
                                                                                                                                               cccATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCaa        >  1:1317581/1‑150 (MQ=255)
                                                                                                                                                 cATTTGCAATGGGTTCCGGGAATGACATAAAGAGAGGAAGGAGCCAGAGCTCTCTCGACGATCAATTGAGTTTCTTCGCCTCGCATCACATTGTGGTTATCGTCATGAGACACACATAATCCCGGAATAATACAGATATTGTCGCCAACt      >  2:472721/1‑150 (MQ=255)
                                                                                                                                                    ttGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCATCTGAc   <  2:1557710/150‑1 (MQ=255)
                                                                                                                                                     tGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACg  <  2:213708/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACG  >  NC_000913/3873899‑3874197

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