Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,138,006 T→G 8.3% T344P (ACG→CCG)  gldA ← glycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,138,0060TG 8.3% 327.5 / 16.9 126T344P (ACG→CCG) gldAglycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (57/58);  new base G (3/8);  total (60/66)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

Only 100 of 128 total aligned reads displayed.

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACGGT  >  NC_000913/4137857‑4138153
                                                                                                                                                     |                                                                                                                                                   
cGTAGGCCTGTTAAGCGTAGCGCTTCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGt                                                                                                                                                     <  1:335025/150‑1 (MQ=255)
 gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTTATCTGGCGTc                                                                                                                                                    >  2:1357736/1‑150 (MQ=255)
 gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTc                                                                                                                                                    <  1:1210669/150‑1 (MQ=255)
     gCCTGATAAGCTTAGCGCATCAGTCAATTTTGATTTAAAACTAAAGGACAAGCAGTGAGTTTTGAGTAGGTTATTAACACTCTTGCAGGAAACGCTGACCGTACTGGTAGTCTACCAGCAGAGCGGCGTAAAACTGATCTGGCGTAgcgc                                                                                                                                                >  2:1245498/1‑150 (MQ=255)
      ccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCAACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTTGCGTCGcccc                                                                                                                                               >  2:1273350/1‑150 (MQ=255)
       cTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGcgccg                                                                                                                                              <  2:85711/150‑1 (MQ=255)
        tGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGcgccgc                                                                                                                                             >  1:1104874/1‑150 (MQ=255)
         gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                                                                                                                            <  2:94936/150‑1 (MQ=255)
         gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                                                                                                                            <  2:75093/150‑1 (MQ=255)
         gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                                                                                                                            <  1:1237878/150‑1 (MQ=255)
         gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                                                                                                                            <  1:513108/150‑1 (MQ=255)
              gCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatc                                                                                                                                       <  2:1072330/150‑4 (MQ=255)
               cGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGCTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatca                                                                                                                                      <  1:65447/150‑5 (MQ=255)
                gaaGCGCATCAGGCAATTTTGCGGTCAAACTCCCGGACCCGCCGGGAGTTTGGAGTAGGGTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatcag                                                                                                                                     <  1:151818/148‑6 (MQ=255)
                gTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtg                                                                                                                                     <  2:641316/150‑1 (MQ=255)
                 tAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGATAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgt                                                                                                                                    <  2:756268/150‑1 (MQ=255)
                 tAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgt                                                                                                                                    <  1:1707123/150‑1 (MQ=255)
                 tAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatcagg                                                                                                                                    <  2:1720708/150‑7 (MQ=255)
                  aGCGCATCAGGCAATTTTGCGTTCAAACTCCCAGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACAGGTCGGCTACCAGCAGAGCTGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgtt                                                                                                                                   >  2:563100/1‑150 (MQ=255)
                   gcgcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgttg                                                                                                                                  <  1:943189/150‑1 (MQ=255)
                     gcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAACCGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGTCGTCGCGCCGCCAGGCATGTtgtg                                                                                                                                >  2:281090/1‑150 (MQ=255)
                           ggCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGt                                                                                                                          >  1:1412627/1‑150 (MQ=255)
                                   tGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCAGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCtt                                                                                                                  >  2:1401000/1‑150 (MQ=255)
                                     cGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCt                                                                                                                >  2:1232603/1‑150 (MQ=255)
                                       ttCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGc                                                                                                              <  1:674757/150‑1 (MQ=255)
                                       tccAAACTCCCGGCTTGTCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGGGAATGGTTTCACCtgctgc                                                                                                              >  1:90236/3‑150 (MQ=255)
                                         cAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGcac                                                                                                            <  1:775663/150‑1 (MQ=255)
                                          aaaCTGCCGGTCCAGCCGGGAGTTTGGAGTAGGTTATTCCCTCTCTTGCAGGAAACGCTGACCGTTCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTTATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGcaca                                                                                                           <  1:1586636/150‑1 (MQ=255)
                                          aaaCTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGcaca                                                                                                           >  1:454590/1‑150 (MQ=255)
                                             cTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCGCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGc                                                                                                        >  1:889972/1‑150 (MQ=255)
                                             cTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAACCCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGc                                                                                                        >  2:843057/1‑150 (MQ=255)
                                                    aCAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCt                                                                                                 <  2:369575/150‑1 (MQ=255)
                                                     cAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTTCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGACAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTg                                                                                                <  2:1433329/150‑1 (MQ=255)
                                                     cAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAGGTTGTGAATGGTTTCACCTACTGCACATGCCGCTTCTg                                                                                                >  1:1453116/1‑150 (MQ=255)
                                                       aGCCGGGAGGTTGGAGTAGGTTACTCCCACTCTCGCAGGAAACGCGGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTGGGGAATGGTTTCACCTTCTGCACATGCCGCTTCTGcc                                                                                              <  1:1305631/150‑1 (MQ=255)
                                                          cGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCcaca                                                                                           >  2:39227/1‑150 (MQ=255)
                                                              aGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTc                                                                                       >  2:596637/1‑150 (MQ=255)
                                                               gTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCg                                                                                      >  2:689230/1‑150 (MQ=255)
                                                                 ttGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCa                                                                                    >  2:469762/1‑150 (MQ=255)
                                                                 ttGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGGCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCa                                                                                    <  1:306310/150‑1 (MQ=255)
                                                                  tGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCAt                                                                                   <  2:532027/150‑1 (MQ=255)
                                                                  tGGAGTAGGTAATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCAt                                                                                   <  1:1515628/150‑1 (MQ=255)
                                                                    gAGTAGGTTATTCCCACTCTTGCAGGAAAAGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCAttt                                                                                 >  1:488626/1‑150 (MQ=255)
                                                                         ggTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGcc                                                                            <  1:1018786/150‑1 (MQ=255)
                                                                              ttCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGAc                                                                       <  2:1429722/150‑1 (MQ=255)
                                                                                  cACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAtct                                                                   <  1:733087/150‑1 (MQ=255)
                                                                                    ctctTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAtcttc                                                                 <  1:927514/150‑1 (MQ=255)
                                                                                       ttGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCttt                                                              <  2:1126531/150‑1 (MQ=255)
                                                                                        tGCAGGAAACGCTGACCGTACTTGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCCCGCCGCCAGGCATGTTGTGAATGGTATAACCTTCTGCACATGCCGCTTCTGACACAATTCGCATTTTCGACGGGACATCTTCTTTa                                                             >  2:1383103/1‑150 (MQ=255)
                                                                                        tGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTa                                                             <  1:1403583/150‑1 (MQ=255)
                                                                                         gCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTat                                                            <  1:653627/150‑2 (MQ=255)
                                                                                          cAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAat                                                           >  2:659955/1‑150 (MQ=255)
                                                                                           aGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAata                                                          <  2:1591832/150‑1 (MQ=255)
                                                                                           aGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAata                                                          <  1:1568159/150‑1 (MQ=255)
                                                                                           aGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCACGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTGTAata                                                          >  1:174471/1‑150 (MQ=255)
                                                                                            ggAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCGCCATTGGCATTTTCGCCGGGCCATCTTCTTTAatat                                                         >  1:1383594/1‑150 (MQ=255)
                                                                                              aaaCGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATcc                                                       <  2:1000085/150‑1 (MQ=255)
                                                                                                 cGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCCCCTTCAGCAAAAGCCGCTTCTGCCACAATTCGAATTTTCGCCGCGACATCTTCTTTAATATCCAGt                                                    >  1:1449751/1‑150 (MQ=255)
                                                                                                       ccGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCg                                              >  2:348040/1‑150 (MQ=255)
                                                                                                       ccGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCg                                              <  2:1724431/150‑1 (MQ=255)
                                                                                                          tACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAGTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgagc                                           >  2:864560/1‑149 (MQ=255)
                                                                                                           aCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTTCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGagag                                          <  2:1103536/150‑1 (MQ=255)
                                                                                                           aCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGagag                                          <  1:281090/150‑1 (MQ=255)
                                                                                                             tGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGtt                                        >  1:515116/1‑150 (MQ=255)
                                                                                                              ggTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGttt                                       >  2:670875/1‑149 (MQ=255)
                                                                                                               gTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTAt                                      >  1:1578857/1‑150 (MQ=255)
                                                                                                               gTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTAt                                      <  2:382772/150‑1 (MQ=255)
                                                                                                                tCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGAGATCTTCTTTAATATCCAGTTGAGCGAGAGTTAtt                                     >  2:621992/1‑150 (MQ=255)
                                                                                                                gcGGCTAGCAGCAGAGCGGCGTACACCTGATCGGGCGTCGCGCCGCCTGGCATGTTGTGATTGGTTTCACCTGTTGGACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTAtt                                     <  1:31605/149‑1 (MQ=255)
                                                                                                                  ggCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTgg                                   <  2:488626/150‑1 (MQ=255)
                                                                                                                  ggCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGAGCCGCCAGGCATGTTTTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATATTCTTTAATATCCAGTTGAGCGAGAGTTATTgg                                   <  1:1273350/150‑1 (MQ=255)
                                                                                                                  ggCGACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTgg                                   <  1:1220999/150‑1 (MQ=255)
                                                                                                                   gCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGc                                  >  1:222274/1‑150 (MQ=255)
                                                                                                                   gCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCATCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGc                                  <  1:1232603/150‑1 (MQ=255)
                                                                                                                    cTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCa                                 <  2:1412627/150‑1 (MQ=255)
                                                                                                                    cTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATACGCATTTTCGCCGCGACATCTACTTTAATATCCAGTAGCGCGAGACTTATTGGCa                                 >  2:81317/1‑150 (MQ=255)
                                                                                                                      aCCAGCAGCGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCaaa                               >  1:1164807/1‑150 (MQ=255)
                                                                                                                        cagcagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAAcc                             <  2:1578857/150‑1 (MQ=255)
                                                                                                                         agcagAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCt                            >  2:978934/1‑150 (MQ=255)
                                                                                                                            agagCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTAcc                         >  2:1582459/1‑150 (MQ=255)
                                                                                                                            agagCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTAcc                         <  2:1731118/150‑1 (MQ=255)
                                                                                                                               gcggcgTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAGTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCa                      >  1:1514188/1‑150 (MQ=255)
                                                                                                                                cggcgTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCAt                     <  1:563100/150‑1 (MQ=255)
                                                                                                                                cggcgTAAACCTGATCTGGCGGCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCAt                     <  1:1141901/150‑1 (MQ=255)
                                                                                                                                  gcgTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGACACAATTCGCATTTTGGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGAATgg                   >  2:477501/1‑150 (MQ=255)
                                                                                                                                   cgTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGc                  <  2:147602/150‑1 (MQ=255)
                                                                                                                                        aacatgccTGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGCCATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAgg             >  1:358915/8‑150 (MQ=255)
                                                                                                                                         ccTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAggg            >  1:199162/1‑150 (MQ=255)
                                                                                                                                         acatgccTGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCGTTTTCGCCGGGGCATCTTCGTTAATATCCAGTTGAGCGAGAGTTATTGGCACACCTACCGCGTGGCTAAggg            >  1:467734/7‑150 (MQ=255)
                                                                                                                                          cTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGc           <  1:469762/150‑1 (MQ=255)
                                                                                                                                            taTCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAg         <  2:313690/149‑1 (MQ=255)
                                                                                                                                            gATCTGGCGTCGCGCCGCCAGGCATGTTGGGAATGGTTTCACCTTCTGCACCTGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAg         <  2:1311467/150‑1 (MQ=255)
                                                                                                                                              tCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCt       >  1:996679/1‑150 (MQ=255)
                                                                                                                                              ccTGGCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAATCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCt       >  2:294805/2‑150 (MQ=255)
                                                                                                                                                 ggcggcGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACg    >  2:953820/1‑150 (MQ=255)
                                                                                                                                                 ggCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACg    <  1:546818/150‑1 (MQ=255)
                                                                                                                                                 ggCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACg    >  2:1098370/1‑150 (MQ=255)
                                                                                                                                                 ggCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACg    >  2:1726130/1‑150 (MQ=255)
                                                                                                                                                   cGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACGGt  >  2:863400/1‑150 (MQ=255)
                                                                                                                                                   cGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACGGt  >  2:1002210/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACGGT  >  NC_000913/4137857‑4138153

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: