Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,874,048 T→G 5.4% 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.4% 306.9 / 7.3 93K143Q (AAA→CAA) dgoK2‑dehydro‑3‑deoxygalactonokinase
Reads supporting (aligned to +/- strand):  ref base T (46/42);  new base G (2/3);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACG  >  NC_000913/3873901‑3874197
                                                                                                                                                   |                                                                                                                                                     
cAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGc                                                                                                                                                     >  2:1536358/1‑150 (MQ=255)
     ccAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATgg                                                                                                                                                >  1:447699/1‑150 (MQ=255)
         cTGTGAAGGCATGGGCAGAGTTTTCCTGCGGCGTCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgc                                                                                                                                            <  2:796275/150‑3 (MQ=255)
          tgtgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTcc                                                                                                                                           >  1:1335772/1‑150 (MQ=255)
           gtgAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCg                                                                                                                                          >  2:1338568/1‑150 (MQ=255)
             gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCggg                                                                                                                                        <  1:922647/150‑1 (MQ=255)
             gAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgcaggc                                                                                                                                        <  2:999246/150‑7 (MQ=255)
              aaGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGc                                                                                                                                       <  2:957286/150‑1 (MQ=255)
              aaGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTGCAAATGGGtgcaggcc                                                                                                                                       <  1:925216/150‑8 (MQ=255)
                   aTCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGAc                                                                                                                                  >  2:468235/1‑150 (MQ=255)
                      ggCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATa                                                                                                                               <  2:16826/150‑1 (MQ=255)
                      ggCAGAGTTTTCCTGCGGCGACAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATa                                                                                                                               >  1:1077821/1‑150 (MQ=255)
                         agagTTTTCCTGCGGCGGCTAACGTGCGGCAATCAATGAGTGATTTAACAGTACATGATGTAATTCAGCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGTCTATCGGCCTGCACCCATTTGCTATGGGTTCCGGGCATGACATACag                                                                                                                            <  2:800040/150‑1 (MQ=255)
                             ttttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAagagag                                                                                                                        <  2:1110053/150‑1 (MQ=255)
                               ttCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCTCCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGagga                                                                                                                      >  1:1121488/1‑150 (MQ=255)
                                tCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaa                                                                                                                     >  2:404076/1‑150 (MQ=255)
                                 ccTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaag                                                                                                                    <  1:176027/150‑1 (MQ=255)
                                 ccTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGaggaag                                                                                                                    <  2:1248986/150‑1 (MQ=255)
                                     cggcggCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGc                                                                                                                <  1:227037/150‑1 (MQ=255)
                                        cggcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCag                                                                                                             <  2:41358/150‑1 (MQ=255)
                                         ggcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGGGAGGAAGGAGCCaga                                                                                                            >  1:150611/1‑150 (MQ=255)
                                          gcAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCagag                                                                                                           <  1:782173/150‑1 (MQ=255)
                                                cTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTcgcg                                                                                                     >  2:466594/1‑150 (MQ=255)
                                                  gcgcCAAGCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAACTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAACAGAGGAAGGAGCCAGAGCTcgcgcg                                                                                                   <  2:1428790/150‑1 (MQ=255)
                                                   cgcCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCgcgcgc                                                                                                  >  1:1336713/1‑150 (MQ=255)
                                                    gcCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGcc                                                                                                 >  2:53645/1‑150 (MQ=255)
                                                    gcCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGcc                                                                                                 <  1:978533/150‑1 (MQ=255)
                                                      caatcaatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGa                                                                                               >  1:1045795/1‑150 (MQ=255)
                                                       aatcaatGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGAt                                                                                              >  2:235906/1‑150 (MQ=255)
                                                           aatTAGAGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCTTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAAt                                                                                          >  2:1396493/1‑150 (MQ=255)
                                                              gAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgt                                                                                       >  1:217553/1‑150 (MQ=255)
                                                              gAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgt                                                                                       >  1:187272/1‑150 (MQ=255)
                                                              gAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgt                                                                                       >  1:181673/1‑150 (MQ=255)
                                                               aGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgtg                                                                                      >  1:315225/1‑150 (MQ=255)
                                                               aGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATtgtg                                                                                      <  2:539746/150‑1 (MQ=255)
                                                                   aTTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTttc                                                                                  >  1:892989/1‑150 (MQ=255)
                                                                        aCAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgc                                                                             <  2:725406/150‑1 (MQ=255)
                                                                        aCAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgc                                                                             >  2:1091597/1‑150 (MQ=255)
                                                                         cAGTAAATGATGTAATTCACCGGTCATCGCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgcc                                                                            >  2:936042/1‑150 (MQ=255)
                                                                         cAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgcc                                                                            >  2:1417786/1‑150 (MQ=255)
                                                                          aGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTcgccg                                                                           <  2:1162779/150‑1 (MQ=255)
                                                                             aaaTGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCgcgc                                                                        >  2:851317/1‑150 (MQ=255)
                                                                              aaTGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCa                                                                       >  2:150754/1‑150 (MQ=255)
                                                                              aaTGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCa                                                                       <  2:247737/150‑1 (MQ=255)
                                                                               atgatgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCAt                                                                      >  2:421969/1‑150 (MQ=255)
                                                                                tgatgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATc                                                                     <  2:1145087/150‑1 (MQ=255)
                                                                                tgatgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATc                                                                     <  2:228557/150‑1 (MQ=255)
                                                                                 gatgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATca                                                                    <  1:73037/150‑1 (MQ=255)
                                                                                  atgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATcac                                                                   <  1:147704/150‑1 (MQ=255)
                                                                                  atgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATcac                                                                   >  1:498539/1‑150 (MQ=255)
                                                                                   tgTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATcaca                                                                  <  1:1536358/150‑1 (MQ=255)
                                                                                      aaTTCACCGGTCATCACGTTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCAATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTATCCGCGCATAACATtg                                                               >  2:929218/1‑150 (MQ=255)
                                                                                      aaTTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATtg                                                               >  2:497567/1‑150 (MQ=255)
                                                                                          cACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCTTTGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGt                                                           <  1:795588/150‑1 (MQ=255)
                                                                                               gTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATTGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCg                                                      >  2:170513/1‑150 (MQ=255)
                                                                                                 catcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTc                                                    <  1:150754/150‑1 (MQ=255)
                                                                                                  atcaCGTTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCa                                                   <  1:678980/150‑1 (MQ=255)
                                                                                                   tcaCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCAt                                                  >  1:185768/1‑150 (MQ=255)
                                                                                                    caCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATg                                                 >  2:302138/1‑150 (MQ=255)
                                                                                                    caCGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATg                                                 >  2:654894/1‑150 (MQ=255)
                                                                                                       ggTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGGGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgaga                                              <  2:150611/150‑1 (MQ=255)
                                                                                                         tGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAACAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATgagaga                                            >  2:1318189/1‑150 (MQ=255)
                                                                                                           cGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGaca                                          >  1:381755/1‑150 (MQ=255)
                                                                                                           cGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGaca                                          <  2:1335772/150‑1 (MQ=255)
                                                                                                            gAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACCTTGTGGTTATCGTCATGAGAGacac                                         <  2:1256865/150‑1 (MQ=255)
                                                                                                             aaaaTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGAcaca                                        <  1:1107073/150‑1 (MQ=255)
                                                                                                             aaaaTCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGAcaca                                        >  2:1342240/1‑150 (MQ=255)
                                                                                                                   gTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATcc                                  >  2:1123742/1‑150 (MQ=255)
                                                                                                                     tGATTTGCTGGCTATCGGCCTGCTCCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCgg                                <  2:1121488/150‑1 (MQ=255)
                                                                                                                     tGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATAGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCgg                                >  2:1440885/1‑150 (MQ=255)
                                                                                                                      gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCggg                               <  2:185768/150‑1 (MQ=255)
                                                                                                                      gATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCggg                               <  1:1338568/150‑1 (MQ=255)
                                                                                                                       aTTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGa                              >  1:1092566/1‑150 (MQ=255)
                                                                                                                            cTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGaataat                         <  1:322421/150‑1 (MQ=255)
                                                                                                                             tGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTGTTCGCCTCGCATCACATTGTGGTTATCGGCAAGAGACACACATAGTACGGGaggaatg                        >  1:661027/1‑144 (MQ=255)
                                                                                                                              ggCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATcc                       <  2:181673/150‑1 (MQ=255)
                                                                                                                              ggCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATcc                       <  2:187272/150‑1 (MQ=255)
                                                                                                                               gCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCACCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCa                      <  2:1336713/150‑1 (MQ=255)
                                                                                                                                cTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAg                     <  1:1091142/150‑1 (MQ=255)
                                                                                                                                cTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCCCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAg                     <  2:447699/150‑1 (MQ=255)
                                                                                                                                 tATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGa                    >  1:347658/1‑150 (MQ=255)
                                                                                                                                  aTCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGat                   <  1:1091597/150‑1 (MQ=255)
                                                                                                                                   tCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGata                  <  1:1417786/150‑1 (MQ=255)
                                                                                                                                      cccggaaCCCATTGCAAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTg               >  2:831291/7‑150 (MQ=255)
                                                                                                                                       ccggaaCCCATTGCAAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGt              >  1:468558/6‑150 (MQ=255)
                                                                                                                                        cTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTc             >  2:1471059/1‑150 (MQ=255)
                                                                                                                                         tGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCg            <  1:468235/150‑1 (MQ=255)
                                                                                                                                           cACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGcc          >  2:573907/1‑150 (MQ=255)
                                                                                                                                           cACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGcc          <  2:892989/150‑1 (MQ=255)
                                                                                                                                             cccATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCaa        >  1:96991/1‑150 (MQ=255)
                                                                                                                                             cccATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCaa        >  1:1272747/1‑150 (MQ=255)
                                                                                                                                             cccATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCACCaa        >  1:387138/1‑150 (MQ=255)
                                                                                                                                                aTTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTg     <  1:497567/150‑1 (MQ=255)
                                                                                                                                                 tttGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCACCAACTGa    <  2:387138/150‑1 (MQ=255)
                                                                                                                                                  ttGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGAc   <  1:404076/150‑1 (MQ=255)
                                                                                                                                                   tGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGATTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACg  >  1:992369/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCACCGGTCATCACGGTGCGAAAATCGTTGATTTGCTGGCTATCGGCCTGCACCCATTTGCAATGGGTTCCGGGCATGACATAAAGAGAGGAAGGAGCCAGAGCTCGCGCGCCGATCAATTGTGTTTCTTCGCCGCGCATCACATTGTGGTTATCGTCATGAGAGACACATAATCCGGGAATAATCCAGATATTGTCGCCAACTGACG  >  NC_000913/3873901‑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: