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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,874,0480TG 5.9% 271.9 / 7.2 85K143Q (AAA→CAA) dgoK2‑dehydro‑3‑deoxygalactonokinase
Reads supporting (aligned to +/- strand):  ref base T (38/42);  new base G (2/3);  total (40/45)
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.99e-01

CTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACG  >  NC_000913/3873899‑3874197
                                                                                                                                                     |                                                                                                                                                     
cTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCAttt                                                                                                                                                       <  1:50879/150‑1 (MQ=255)
cTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTAACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCAttt                                                                                                                                                       >  2:51817/1‑150 (MQ=255)
 tCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTg                                                                                                                                                      <  1:53254/150‑1 (MQ=255)
 tCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTg                                                                                                                                                      <  1:1104677/150‑1 (MQ=255)
  cAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGc                                                                                                                                                     <  1:285258/150‑1 (MQ=255)
       ccAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGAGTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATgg                                                                                                                                                <  2:1299850/150‑1 (MQ=255)
        cAGCTGTGAAGGCATCCGCAGAGTTTTCCTGCGGCGTCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATggg                                                                                                                                               <  2:819247/150‑1 (MQ=255)
           cTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgc                                                                                                                                            <  1:1599136/150‑3 (MQ=255)
            tgtgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTcc                                                                                                                                           >  1:164764/1‑150 (MQ=255)
            tgtgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgca                                                                                                                                           <  1:1339219/150‑4 (MQ=255)
             gtgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCg                                                                                                                                          >  1:1246350/1‑150 (MQ=255)
               gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCggg                                                                                                                                        >  2:584817/1‑150 (MQ=255)
               gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCggg                                                                                                                                        <  1:1347492/150‑1 (MQ=255)
               gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCggg                                                                                                                                        >  2:1101524/1‑150 (MQ=255)
               gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCggg                                                                                                                                        <  1:1002607/150‑1 (MQ=255)
               gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgcaggc                                                                                                                                        <  1:1021934/150‑7 (MQ=255)
                   gCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATg                                                                                                                                    >  2:510747/1‑150 (MQ=255)
                    cATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGa                                                                                                                                   >  2:29773/1‑150 (MQ=255)
                    cATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGa                                                                                                                                   <  1:968772/150‑1 (MQ=255)
                        ggCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATATGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATa                                                                                                                               <  2:1457905/150‑1 (MQ=255)
                              gTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagaga                                                                                                                         <  1:827794/150‑1 (MQ=255)
                                 ttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGagga                                                                                                                      <  1:471493/150‑1 (MQ=255)
                                  tCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaa                                                                                                                     >  1:4238/1‑150 (MQ=255)
                                   ccTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaag                                                                                                                    >  1:968039/1‑150 (MQ=255)
                                    cTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaagg                                                                                                                   >  1:117938/1‑150 (MQ=255)
                                    cTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCAACCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaagg                                                                                                                   <  2:1434570/150‑1 (MQ=255)
                                     tGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaagga                                                                                                                  >  2:300921/1‑150 (MQ=255)
                                        ggcggcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGcc                                                                                                               <  1:547671/150‑1 (MQ=255)
                                           ggcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCaga                                                                                                            <  1:412441/150‑1 (MQ=255)
                                                  cTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTcgcg                                                                                                     <  1:867715/150‑1 (MQ=255)
                                                   tGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGGTcgcgc                                                                                                    <  2:1371713/150‑1 (MQ=255)
                                                   tGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTcgcgc                                                                                                    >  1:1582321/1‑150 (MQ=255)
                                                    gcgcCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTcgcgcg                                                                                                   >  1:1547539/1‑150 (MQ=255)
                                                         aatcaatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGAt                                                                                              >  1:431394/1‑150 (MQ=255)
                                                          atcaatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATc                                                                                             >  1:1288364/1‑150 (MQ=255)
                                                           tcaatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCa                                                                                            >  2:650628/1‑150 (MQ=255)
                                                             aatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAAt                                                                                          >  1:944645/1‑150 (MQ=255)
                                                                gAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGa                                                                                                                                   >  2:111753/1‑106 (MQ=255)
                                                                gAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGa                                                                                                                                   <  1:111753/106‑1 (MQ=255)
                                                                gAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgt                                                                                       >  1:1365071/1‑150 (MQ=255)
                                                                 aGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgtg                                                                                      <  2:52028/150‑1 (MQ=255)
                                                                  gTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTgtgt                                                                                     <  2:431394/150‑1 (MQ=255)
                                                                       ttAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTttctt                                                                                <  1:1101524/150‑1 (MQ=255)
                                                                        tAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTttcttc                                                                               <  1:764975/150‑1 (MQ=255)
                                                                         aaCAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcg                                                                              <  1:218364/150‑1 (MQ=255)
                                                                          aCAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgc                                                                             >  2:608764/1‑150 (MQ=255)
                                                                           cAGTAACTGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCAGGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgcc                                                                            <  1:1457964/150‑1 (MQ=255)
                                                                           cAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgcc                                                                            >  1:198559/1‑150 (MQ=255)
                                                                           cAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgcc                                                                            >  2:1355729/1‑150 (MQ=255)
                                                                           cAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgcc                                                                            >  2:1090021/1‑150 (MQ=255)
                                                                            aGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgccg                                                                           <  1:1404171/150‑1 (MQ=255)
                                                                            aGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgccg                                                                           >  2:1044432/1‑150 (MQ=255)
                                                                            aGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGATCGCGAGCCGATCAATTGTGTTTCTTcgcag                                                                           >  2:562101/1‑148 (MQ=255)
                                                                             gTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgccgc                                                                          <  2:968039/150‑1 (MQ=255)
                                                                             gTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgccgc                                                                          >  2:1327105/1‑150 (MQ=255)
                                                                               aaaTGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCgcgc                                                                        <  2:495614/150‑1 (MQ=255)
                                                                                  tgatgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATc                                                                     <  2:567888/150‑1 (MQ=255)
                                                                                    atgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATcac                                                                   >  1:736303/1‑150 (MQ=255)
                                                                                          ttCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATtgtg                                                             <  1:510747/150‑1 (MQ=255)
                                                                                           tCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTgg                                                            <  1:902473/150‑1 (MQ=255)
                                                                                             aCCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGtt                                                          >  1:756404/1‑150 (MQ=255)
                                                                                                ggTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATc                                                       <  2:1582321/150‑1 (MQ=255)
                                                                                                 gTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCg                                                      <  1:1044432/150‑1 (MQ=255)
                                                                                                  tcatcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGt                                                     <  2:1365071/150‑1 (MQ=255)
                                                                                                     tcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCAt                                                  >  2:261064/1‑150 (MQ=255)
                                                                                                      caCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATg                                                 <  1:608764/150‑1 (MQ=255)
                                                                                                          gTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgagag                                             <  2:1617945/150‑1 (MQ=255)
                                                                                                              gAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGacac                                         >  2:1288936/1‑150 (MQ=255)
                                                                                                                        gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgag                                               >  1:1146015/1‑134 (MQ=255)
                                                                                                                        gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgag                                               <  2:1146015/134‑1 (MQ=255)
                                                                                                                        gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCggg                               <  2:944645/150‑1 (MQ=255)
                                                                                                                             gCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGaataa                          >  2:1529053/1‑150 (MQ=255)
                                                                                                                             gCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGaataa                          >  1:546864/1‑150 (MQ=255)
                                                                                                                                ggCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATcc                       <  2:736303/150‑1 (MQ=255)
                                                                                                                                    aTCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGat                   <  2:1288364/150‑1 (MQ=255)
                                                                                                                                     tCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGCGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACAAATAATCCGGGAATAATCCAGata                  >  1:1060136/1‑150 (MQ=255)
                                                                                                                                     tCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGata                  <  1:1029775/150‑1 (MQ=255)
                                                                                                                                      cGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGatat                 >  1:175675/1‑150 (MQ=255)
                                                                                                                                        cccggaaCCCATTGCAAATGGGTTCCGGGCATGACATAAAGATAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTg               >  2:45488/7‑150 (MQ=255)
                                                                                                                                           tGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGaataa                          >  1:923694/1‑136 (MQ=255)
                                                                                                                                           tGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGaataa                          <  2:923694/136‑1 (MQ=255)
                                                                                                                                           tGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCg            >  2:189287/1‑150 (MQ=255)
                                                                                                                                           tGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCg            <  2:1246350/150‑1 (MQ=255)
                                                                                                                                            gaaCCCATTGCAAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGc           >  1:1099336/3‑150 (MQ=255)
                                                                                                                                            gCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGc           <  1:536676/150‑1 (MQ=255)
                                                                                                                                               cccATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCaa        >  1:1303023/1‑150 (MQ=255)
                                                                                                                                                 cATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACt      >  1:1522298/1‑150 (MQ=255)
                                                                                                                                                   tttGCAATGGGTTCCGGGCATGACATAAAGCGAGGCAGGAGCCAGAGCTCGCGCGCCGAGCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGa    <  2:1051293/150‑1 (MQ=255)
                                                                                                                                                     tGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACg  <  1:1466873/150‑1 (MQ=255)
                                                                                                                                                     tGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACg  <  2:1352360/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 < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: