Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,593,431 A→G 5.2% V431A (GTC→GCC) ‡ lgoT ← putative L‑galactonate:H+ symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,593,4310AG 5.2% 342.6 / 7.0 115V431A (GTC→GCC) ‡lgoTputative L‑galactonate:H+ symporter
Reads supporting (aligned to +/- strand):  ref base A (44/64);  new base G (2/4);  total (47/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

Only 100 of 116 total aligned reads displayed.

GAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACG  >  NC_000913/4593294‑4593580
                                                                                                                                         |                                                                                                                                                     
taaCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAga                                                                                                                                           <  2:1692433/149‑1 (MQ=255)
 aaCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgcg                                                                                                                                          <  1:1368096/150‑3 (MQ=255)
 aaCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgat                                                                                                                                          <  2:803914/150‑1 (MQ=255)
 aaCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGAGCCCCCGCTGCGGTTGCGTCACCGCAgcg                                                                                                                                          <  1:1183881/150‑3 (MQ=255)
    tGGGAAACATTTTAATCGCCATGAAAATGCCCGATCAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatgat                                                                                                                                       <  2:1249555/150‑1 (MQ=255)
    tGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgcgggg                                                                                                                                       <  2:851726/150‑6 (MQ=255)
    cgggAAACTTTTTAATCGCCCTGAATCTGCCCGATAACCAAAATGTTATCGGGCATAAGGTGACTAATCTTTACGTGGGTCGTTGATGGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatgat                                                                                                                                       <  2:1243238/149‑1 (MQ=255)
     gggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTTTCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACTCCTGTGCCCCCGCTGCGGTGACGCAACCGCTGATGAtt                                                                                                                                      >  1:360500/1‑150 (MQ=255)
     gggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGAtt                                                                                                                                      <  2:976877/150‑1 (MQ=255)
     gggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGAtt                                                                                                                                      <  2:389897/150‑1 (MQ=255)
     gggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGATGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGAtt                                                                                                                                      >  2:1887165/1‑150 (MQ=255)
          aCATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGc                                                                                                                                 <  2:1429676/150‑1 (MQ=255)
           cattgaaaTCGCCATGAAAATGCCCGATAAGCAAAATGTAAGCGCGCATAAGGAGATTAATCTTTACGTGGGTCGTTGCTCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGcc                                                                                                                                <  1:372030/144‑1 (MQ=255)
            aTTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCa                                                                                                                               >  2:1160286/1‑150 (MQ=255)
            aTTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCa                                                                                                                               >  2:1172891/1‑150 (MQ=255)
             ttttAATCGCCATGAAAATGCCCGATATGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAg                                                                                                                              >  1:1789783/1‑150 (MQ=255)
                 aaTCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACgg                                                                                                                          <  2:399900/150‑1 (MQ=255)
                     gCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATg                                                                                                                      <  1:475367/150‑1 (MQ=255)
                      ccATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGCCGCAACCGCAGATGATTAGTGCCAGACGGAATGa                                                                                                                     >  1:682951/1‑150 (MQ=255)
                         tGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCTGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCACCCGCAGATGATTAGTGCCAGACCGAATTAgtg                                                                                                                  >  2:371914/1‑150 (MQ=255)
                         tGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtg                                                                                                                  >  1:1329848/1‑150 (MQ=255)
                            aactGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCATTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtgggt                                                                                                               <  1:987640/147‑1 (MQ=255)
                                 cccGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACTAGGAAGATGTACGCCAGTGCCCACGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTAt                                                                                                          >  1:1591429/1‑150 (MQ=255)
                                  ccGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATc                                                                                                         <  1:29186/150‑1 (MQ=255)
                                   cGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATca                                                                                                        <  2:1888577/150‑1 (MQ=255)
                                    caTAAGCAAAATGTTATCGGGAATAAGGAGCTTAATCTTTACTTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGGCCCCGCTGCGGGGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATcaa                                                                                                       <  2:1786893/149‑1 (MQ=255)
                                     aTAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGGTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATcaac                                                                                                      <  2:559656/150‑1 (MQ=255)
                                      tAAGCAAACTGTTATCGGGCATAAGGAGATTAATCTTAACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATcaaca                                                                                                     <  1:116576/150‑1 (MQ=255)
                                       acgCAAAATGTTATCGGGCATAAGGAGTTTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCTCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATcaacaa                                                                                                    <  1:1072306/148‑1 (MQ=255)
                                        aGCAAAATTTAATCGGGAATAAGGAGAATAATCTTAACGTGGGTCGTTGAACGACTGACGAACCAGGAAGATGTACGCCAGAGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAAt                                                                                                   <  1:1760222/150‑1 (MQ=255)
                                        aGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAAt                                                                                                   <  2:23564/150‑1 (MQ=255)
                                         gCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                  >  2:1841253/1‑150 (MQ=255)
                                         gCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                  >  1:840993/1‑150 (MQ=255)
                                         gCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                  >  1:673362/1‑150 (MQ=255)
                                         gCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                  >  2:1166480/1‑150 (MQ=255)
                                              aTGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACTAAACCGCAGATGATTAGTGCCAGACGGAATGAGTCGGTTTTGTCAACAAAAAAAcc                                                                                             >  1:946638/1‑150 (MQ=255)
                                              aTGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAAcc                                                                                             <  2:938159/150‑1 (MQ=255)
                                                 ttATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGt                                                                                          >  2:1869123/1‑150 (MQ=255)
                                                  tATCGGGCATAAGGAGATTAATCTTTACGTTGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCGGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTa                                                                                         <  1:987092/150‑1 (MQ=255)
                                                   aTCGGGCATAAGGAGATAAATCTTTACGTGGGTCGTTGATCGGCTGACGAATCAGGAAGATGTACGCCAGTGCCCCCTCTGCGGTGACGCAAACGCAGATGATTAGTGACAGACGCAATGAGTGGGTGGTATCAACAATAAAACCAGTaa                                                                                        <  1:355226/150‑1 (MQ=255)
                                                     cGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATg                                                                                      <  2:1571565/150‑1 (MQ=255)
                                                           tAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGgc                                                                                <  2:129388/150‑1 (MQ=255)
                                                            aaGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGgcg                                                                               >  2:820008/1‑150 (MQ=255)
                                                            aaGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGgcg                                                                               <  1:1121091/150‑1 (MQ=255)
                                                            aaGGAGATTAATCTTTACGTGGGTCGGTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGgcg                                                                               <  1:371304/150‑1 (MQ=255)
                                                             aGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGgcgc                                                                              >  2:1498903/1‑150 (MQ=255)
                                                             aGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCACCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTCGTATCAACAATAAAACCAGTAATGATCGtcgc                                                                              >  1:1632625/1‑150 (MQ=255)
                                                              ggAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCa                                                                             >  1:1819484/1‑150 (MQ=255)
                                                               tagaTTAAGCTTTACGTGGGTCGTTGATCGGCATACGAACCAGGAAGATTTACGCCAGTGCCCCCGCTGTGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGAAGCAACAATAAAACCAGAACTGATCGGCGCaa                                                                            <  1:311706/149‑1 (MQ=255)
                                                                 gaTTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCCACCGCAGATGATTAGTGCCAGCCGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGAaaaa                                                                          >  2:72565/1‑150 (MQ=255)
                                                                 gaTTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAACCCAGTAATGATCGGCGCaaaa                                                                          >  2:759822/1‑150 (MQ=255)
                                                                 gaTTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCaaaa                                                                          <  2:513208/150‑1 (MQ=255)
                                                                  aTTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAag                                                                         <  2:1299984/150‑1 (MQ=255)
                                                                     aaTCTTTACGTGGGTCGTTGATCTGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGAGCGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAgg                                                                      >  1:1471473/1‑150 (MQ=255)
                                                                     aaTCTTAACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAgg                                                                      <  2:994397/150‑1 (MQ=255)
                                                                        cTTTACGTGGGACGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAAACGCAGAAGATTAGAGCCAGCCGGAATGACTGGGTGGTATCCACAATACACCCCGTAATGATCGGCGCAAAAGCGgcgc                                                                   >  2:1083742/1‑150 (MQ=255)
                                                                          ttACGTGGGTCGTGGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAg                                                                 <  2:336787/150‑1 (MQ=255)
                                                                               tGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGaa                                                            <  2:32858/150‑1 (MQ=255)
                                                                                 gttCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGc                                                          <  1:246414/148‑1 (MQ=255)
                                                                                  gTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCACCCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCt                                                         >  1:1301628/1‑150 (MQ=255)
                                                                                   gcGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTg                                                        <  2:1709274/149‑1 (MQ=255)
                                                                                    cGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTgg                                                       <  2:1426657/150‑1 (MQ=255)
                                                                                     gTTGATCGGCTGACGTACCAGGATGATGGACCCCAGGGCCCCGGCTGCGGGTACGCAACCGAATATGAGTTGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAATAGAGGCGCAGATGAAGCTGGc                                                      <  1:490583/150‑1 (MQ=255)
                                                                                       tgctCGGCTGACGAACCAGGAAGATGTACGCCAGAGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCaa                                                    <  1:1799410/147‑1 (MQ=255)
                                                                                         aTCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAg                                                  >  2:209905/1‑150 (MQ=255)
                                                                                              cTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTg                                             <  2:1110471/150‑1 (MQ=255)
                                                                                              cTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGGGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTg                                             <  2:1167768/150‑1 (MQ=255)
                                                                                               tGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                            <  2:197281/150‑1 (MQ=255)
                                                                                                gACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGa                                           <  1:1417270/150‑1 (MQ=255)
                                                                                                 aCGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGAt                                          <  1:1163501/150‑1 (MQ=255)
                                                                                                  cGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATa                                         >  2:1825562/1‑150 (MQ=255)
                                                                                                      ccAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGc                                     <  2:1718633/150‑1 (MQ=255)
                                                                                                       cAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGcc                                    >  2:1476299/1‑150 (MQ=255)
                                                                                                         ggAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCa                                  <  1:1060386/150‑1 (MQ=255)
                                                                                                              aTGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCACCCCGAc                             >  2:1840806/1‑150 (MQ=255)
                                                                                                               tGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACg                            >  2:1667153/1‑150 (MQ=255)
                                                                                                                gTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGc                           <  1:1732027/150‑1 (MQ=255)
                                                                                                                gTACGCCAGTCCCCCCGCGGCGGGGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGTTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGc                           <  2:360500/150‑1 (MQ=255)
                                                                                                                  aCGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAg                         <  2:1840281/150‑1 (MQ=255)
                                                                                                                   cGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGCCGGAGt                        >  2:814223/1‑150 (MQ=255)
                                                                                                                   cGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGCCGCAGt                        >  1:729986/1‑150 (MQ=255)
                                                                                                                   cGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGt                        >  1:686876/1‑150 (MQ=255)
                                                                                                                   cGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGt                        >  1:1846373/1‑150 (MQ=255)
                                                                                                                        gTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGc                   <  2:218923/150‑1 (MQ=255)
                                                                                                                         tGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCg                  <  1:833194/150‑1 (MQ=255)
                                                                                                                           cccccGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCgag                <  1:424182/150‑1 (MQ=255)
                                                                                                                            ccccGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCgaga               >  1:1587386/1‑150 (MQ=255)
                                                                                                                            ccccGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCgaga               >  2:111626/1‑150 (MQ=255)
                                                                                                                            ccccGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGATCTGGATACTGCCCACCGACGCAGTCATGCgaga               >  1:196296/1‑150 (MQ=255)
                                                                                                                            ccccGCGGCGGTGACGCAACCGCAGAAGATTAGTGCCAGTCGGAATGAGTGTGTGGTATCAACAAAAATACCAGTAAAGATCGGCGCAAACGAGGCGCATATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCgaga               <  1:1727573/150‑1 (MQ=255)
                                                                                                                             cccGCTGAGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGaa              >  2:909502/1‑150 (MQ=255)
                                                                                                                               atcTGCGGTTGCGTCACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGc            >  1:1080628/3‑150 (MQ=255)
                                                                                                                                tcTGCGGTTGCGTCACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGCCGCAGTCATGCGAGAAGcc           >  1:1938920/2‑149 (MQ=255)
                                                                                                                                gCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCa           >  2:208717/1‑150 (MQ=255)
                                                                                                                                 cTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCaa          >  1:1047029/1‑150 (MQ=255)
                                                                                                                                  tGCGGTGACGCAACCGCAGATGATTAGTGCCAGACTGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAAc         <  1:174684/150‑1 (MQ=255)
                                                                                                                                    gggTGACGCAACCGCAGATGATTAGTGCCATACGGAATGTGTGGGTGGTATCAACAATAAAACGACTATTGATCGGCGCACAAGAGTCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCTTGCGTGAAGCAACTg       <  2:11449/149‑1 (MQ=255)
                                                                                                                                     ggTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGc      >  1:1744860/1‑150 (MQ=255)
                                                                                                                                      gTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCg     >  2:1550019/1‑150 (MQ=255)
                                                                                                                                         aCGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACg  >  1:1480684/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                                     
GAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACG  >  NC_000913/4593294‑4593580

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: