Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,689,733 G→C 6.4% R177G (CGT→GGT)  bcsZ ← endo‑1,4‑D‑glucanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,689,7330GC 6.4% 261.0 / 7.9 94R177G (CGT→GGT) bcsZendo‑1,4‑D‑glucanase
Reads supporting (aligned to +/- strand):  ref base G (36/52);  new base C (4/2);  total (40/54)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

TCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGC  >  NC_000913/3689586‑3689880
                                                                                                                                                   |                                                                                                                                                   
tctcATAGCGCACCCAGTCTGGCGAAAAACATTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGCTCCACGGCGCGCCAAAGCGGTTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGc                                                                                                                                                    >  1:41583/1‑149 (MQ=255)
 ctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGGTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGcc                                                                                                                                                   <  1:502590/149‑1 (MQ=255)
  tcATAGCGCACCCAGTCTGGCGAAAACCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCa                                                                                                                                                  <  1:723221/149‑1 (MQ=255)
   cATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGTGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGCGTTTTAACCCCAg                                                                                                                                                 <  2:959616/149‑1 (MQ=255)
    aTAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGc                                                                                                                                                >  1:707140/1‑149 (MQ=255)
       gcgcACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTTGTTTCGAGCATCGTGGTCCACGGAGCGCCAAAGCGGGTGAAATACTGCCCCATCGTCGCAGGCAGGTAGCTGGGATTAAAACGcaacctaa                                                                                                                                             >  1:1028155/1‑142 (MQ=255)
       gcgcACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGCGTTTTAACCCCAGCtac                                                                                                                                             <  1:873264/149‑3 (MQ=255)
          cACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGGGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTAt                                                                                                                                          >  1:209598/1‑149 (MQ=255)
          cACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTAt                                                                                                                                          >  1:658360/1‑149 (MQ=255)
             ccAGTCTGGCGAAAAACCTTTCGGGCCGGTTTCCAGCAATAATCGTTGATGGGTTGCGCGCAGCGTGGTTCTCGGCGCGCCAAAGCGGGTGAAATACTTCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCt                                                                                                                                       <  2:947621/149‑1 (MQ=255)
              cAGTCTGTCGAACAACCTTTCGGCGAGGTTTCCAGAAATAAACGTTGATTGGGTTCGCGCAGCGTGGTCAACGGCGCGAAAAAGCGGTTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTc                                                                                                                                      <  1:947446/149‑1 (MQ=255)
               aggcTGGCGAAAAACCTTTCGGGGCGGATGCCAGCAATAACCGTGGATTGGTCTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCa                                                                                                                                     <  1:774229/146‑1 (MQ=255)
               aGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCa                                                                                                                                     <  1:781468/149‑1 (MQ=255)
                 tCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGc                                                                                                                                   <  1:142933/149‑1 (MQ=255)
                 gcTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGc                                                                                                                                   <  1:1001175/148‑1 (MQ=255)
                     gCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCaaac                                                                                                                               >  1:516255/1‑148 (MQ=255)
                     gCGAAAAACCTTTCGGGCCGGTTTCCAGCAATAAACGTTGATTGGTATCGCGCAGCGTGCTCCACGGCGCGCCAAAGCGGGAGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCACCTGTTATACTCAGCaaaa                                                                                                                               >  2:321583/1‑149 (MQ=255)
                      cGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAAc                                                                                                                              <  1:1019035/149‑1 (MQ=255)
                       gAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAAcc                                                                                                                             >  2:1010667/1‑149 (MQ=255)
                             ccTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACttt                                                                                                                       <  2:812664/149‑1 (MQ=255)
                                ttCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAACCGCCAGCTGTTATCCTCCGCAACACCCCCTTTGcc                                                                                                                    >  2:945088/1‑149 (MQ=255)
                                  cGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCg                                                                                                                  <  1:403401/149‑1 (MQ=255)
                                    gtgCGGCTTCCACCAATAAACGTTGATTGGTTGCGCGCAGCGTGGGCCACGGCTCTCCAAAGCGGGTGAAATACTGCGCCAGCGTCTGCGGCAGGTAGCTGGGGTTATAACGCCAGCTGTTATCCTCAGAAAAACCCACTTTGCCCGGt                                                                                                                <  2:903946/147‑1 (MQ=255)
                                    gggCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGt                                                                                                                <  1:748247/149‑1 (MQ=255)
                                     ggCGGTTTCCAGCAATAAACGTGGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTTGGGTTAACACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTa                                                                                                               >  1:736749/1‑149 (MQ=255)
                                      gCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaa                                                                                                              >  1:404258/1‑149 (MQ=255)
                                       cGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGTAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaac                                                                                                             <  2:645323/149‑1 (MQ=255)
                                       cGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaac                                                                                                             <  1:474711/149‑1 (MQ=255)
                                       cGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGACAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaac                                                                                                             <  1:474698/149‑1 (MQ=255)
                                          tttCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaacaac                                                                                                          <  2:441178/149‑1 (MQ=255)
                                          tttCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGACGGCGGCAGGTAGCTGGGGTTAATCCGCCAGCTGTTATCCTCAGCCAACCCCACTTGGCCCGGTaacacc                                                                                                          >  1:740934/1‑147 (MQ=255)
                                            tCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCTCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAACCGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTTACAACAt                                                                                                        >  2:511152/1‑149 (MQ=255)
                                             ccAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCAGTTATCCTCAGCAAAACCCACTTTGCCCGGTAacaccaag                                                                                                       >  1:420365/1‑147 (MQ=255)
                                                  aaTAAACGTTGATTGGTTTCGCGCAGCGTTGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcc                                                                                                  <  2:404258/149‑1 (MQ=255)
                                                  aaTAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGGCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcc                                                                                                  <  1:220054/149‑1 (MQ=255)
                                                    tAAACGTGGATTGGTGTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGGGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAccca                                                                                                <  2:149703/149‑1 (MQ=255)
                                                      accGTTGAATGGATTCGCGCAGCGTGGTCCCCGGCGCGACAACGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGGTAAAAAGCCAGCTGTGAGCCTCCGCAACCCCCACTTTGGCCGGTCACTCCAAGGCAccccgc                                                                                              >  2:174783/3‑149 (MQ=255)
                                                       aCGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAccccgcc                                                                                             >  1:778395/1‑149 (MQ=255)
                                                          ttGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAccccgcccag                                                                                          >  1:543554/1‑149 (MQ=255)
                                                          ttGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcccagcccag                                                                                          <  2:225998/149‑1 (MQ=255)
                                                           tGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAATCCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAgg                                                                                         <  1:169999/149‑1 (MQ=255)
                                                             aTTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCTGCTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTATAAACGACAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGAAACAACATGGAACCCAGCCGAGGCa                                                                                       >  2:570318/1‑149 (MQ=255)
                                                                  tttCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAACCCCACTTTGCCCGGTAACCACATGGAACCCAGCCCCGACCCCGTc                                                                                  >  2:857728/1‑149 (MQ=255)
                                                                  tttCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTc                                                                                  <  1:846052/149‑1 (MQ=255)
                                                                  tttCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTc                                                                                  <  1:43978/149‑1 (MQ=255)
                                                                    tCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCATAACCCACTTTTCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTcac                                                                                <  2:272331/149‑1 (MQ=255)
                                                                    tCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTcac                                                                                <  1:986970/149‑1 (MQ=255)
                                                                     agcgcAGCGTGGTCCACGGCGCGCCAAAGCTGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGGTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTcacc                                                                               <  2:207206/148‑1 (MQ=255)
                                                                      gcgcAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAACCGCCAGCTGTTATCCTCAGCAAACCCCACTTTGCCCGGTAACAACATGGACCCCAGCCCAGGCACCGTcaccc                                                                              >  2:33999/1‑148 (MQ=255)
                                                                           gCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTtcc                                                                         <  2:707140/149‑1 (MQ=255)
                                                                             gTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTtcctc                                                                       <  1:460507/149‑1 (MQ=255)
                                                                                gTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATACTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCcg                                                                    >  1:556101/1‑149 (MQ=255)
                                                                                          cgcCAAATAGTGTGAAATTCTGCGCCAGAGTCTGCGGCAGGTATCTGGTGTTAAAACGCCATCTTTTATCCTCAGCAAAACCCACTTTGCCCGGTAAAAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCAGCGCTATACGt                                                          <  1:23523/149‑1 (MQ=255)
                                                                                            ccAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCCGCAAACCCCCCTGTGCCCGGTAACAACAGGGACCCCACCCCAGGCACAGTCACCACTTCCTCCCGCGCGATActgtt                                                        >  2:160510/1‑145 (MQ=255)
                                                                                            ccAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGttt                                                        >  2:212577/1‑149 (MQ=255)
                                                                                             cAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCCCTTTGCCCGGTAACAACATGGAACCCAGCCCCGGCACCGTCCCCACTTCCTCCCGCGCGATACGtttt                                                       >  1:633808/1‑149 (MQ=255)
                                                                                                  cGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCaa                                                  <  1:229734/149‑1 (MQ=255)
                                                                                                   gggTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGATGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAAc                                                 <  2:761364/149‑1 (MQ=255)
                                                                                                    ggTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAACCCCACTTTGCCCGGTAACAACATGGAACCCAGCCCCGGCACCGTCCCCCCTTCCTCCCGCGCGCTACGTTTTAGCCAcg                                                >  1:439462/1‑149 (MQ=255)
                                                                                                    ggTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCTACTTGGCCCGGTAACAACATGGAACCCAGCCCCGGCGCCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAAcg                                                >  2:106833/1‑149 (MQ=255)
                                                                                                    ggTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAAcg                                                <  1:877381/149‑1 (MQ=255)
                                                                                                       gAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGCACCCAGCCCAGGCACCGTCTCCACTTCCTCCCGCGCGAGACGTTTTAGCAACgcgc                                             >  1:414527/1‑149 (MQ=255)
                                                                                                          aTACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGc                                          <  2:664284/149‑1 (MQ=255)
                                                                                                           tACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAACTCCCACTTTGCCCGGTACCAACATGGATGCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGGTTTAGCAACGCGCTGcc                                         >  1:52037/1‑149 (MQ=255)
                                                                                                              tGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACGTCCTCCCGCGCGATACGTTTTAGTAACGCGCTGCCGAt                                      >  1:984908/1‑149 (MQ=255)
                                                                                                                 gcCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACACCGCCCAGGCACCGTCCCCACTGCCTCCCGCGCGGTACCTGTTAGCGACGCGCTGCCGaaggc                                   >  2:332073/1‑145 (MQ=255)
                                                                                                                  ccATCGTCGTCGGCAGGTAGCTGGGGTAAAAACGCCAGCTGATATCCTCAGCAAAACCCACTTTGCCCTGTAACATCATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCg                                  <  1:672900/149‑1 (MQ=255)
                                                                                                                  ccAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCg                                  <  2:658360/149‑1 (MQ=255)
                                                                                                                  ccAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCg                                  <  1:206116/149‑1 (MQ=255)
                                                                                                                  ccAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTAACCTCAGCAAAACCCACTTTCCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCg                                  <  1:590749/149‑1 (MQ=255)
                                                                                                                  ccAGAGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCTGCAAAACCCACTTTGCCCGGTTACAACATGGAACCCAGCCGAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACTCGCTGCCGATGTCg                                  <  1:321583/149‑1 (MQ=255)
                                                                                                                    aGCGTCGGCGGCAGGTAGCGGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGACACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGt                                <  1:857728/149‑1 (MQ=255)
                                                                                                                          ggcggcAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCCCCGGCCCCCCTTCCTCCCGCGCGATACGGTGTAGCAACGCGCTGCCAATGTCGGGAGAACg                          >  1:824034/1‑149 (MQ=255)
                                                                                                                            cggcAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACgct                        <  2:367010/149‑1 (MQ=255)
                                                                                                                               ataacaGCTGGCGTTTTAACCCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACgctgct                     >  2:622844/6‑149 (MQ=255)
                                                                                                                                 ggTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAACCCCAATTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGctct                   <  2:452828/149‑1 (MQ=255)
                                                                                                                                 aacaGCTGGCGTTTTAACCCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCCCTTCCTCCCGCGCGATACGTTTTAGCCACGCGCTGCCGATGTCGGTATAACGCTGctct                   >  2:950072/4‑149 (MQ=255)
                                                                                                                                   tAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCttt                 <  1:348460/149‑1 (MQ=255)
                                                                                                                                   tAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCttt                 <  1:815595/149‑1 (MQ=255)
                                                                                                                                   caGCTGGCGTTTTAACCCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACAGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGTTGTCGGTATAACGCTGCTCttt                 >  2:146476/2‑149 (MQ=255)
                                                                                                                                    aGCTTGGGTTAAAACGCCAGCTGTTATCCTCAGCAACAGCCACTTTGCCCGGAAACATCATGGAACCCAGCCCAGGCACCGTCACCACTTCGTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTc                <  2:210982/149‑1 (MQ=255)
                                                                                                                                    aGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTc                <  1:831775/149‑1 (MQ=255)
                                                                                                                                    aGCTGGCGTTTTAACCCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGACCCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTc                >  1:883011/1‑149 (MQ=255)
                                                                                                                                     gCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAACCCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACAGTCACCCCTTCCTCCCGCGCGATACTTTTTAGCAACGCGCTGCCGATGGCGGTATAGCGCTGCTCTTTcc               >  1:129544/1‑149 (MQ=255)
                                                                                                                                       gggggCTAAAACGCCAGCTGTTACCCTCATCAAAACCCACTGTGCCCGGTAACAACATGGAACCCAGCCCAGGGACCGTCACCACTTCCTCCCGCGCGCTACGTGTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCcac             <  2:825047/148‑1 (MQ=255)
                                                                                                                                        ggggTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCcaca            >  2:297908/1‑149 (MQ=255)
                                                                                                                                          ggTTAAAACGCCAGCTGTTATCCTCAGCAAACCCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACaaa          <  2:152090/149‑1 (MQ=255)
                                                                                                                                          ggTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCGGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACaaa          <  2:660700/149‑1 (MQ=255)
                                                                                                                                           gTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACGGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAAc         >  2:88938/1‑149 (MQ=255)
                                                                                                                                           gTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCCCCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAAc         >  1:721858/1‑149 (MQ=255)
                                                                                                                                                aaCGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTACCAGCATGGCACCCAGCCCAGGCTCCGGCACCACGTCCTCCCGCGCGATCCGTTCTCGCACCGCGCTGCCGAGGGCGGTAGATCGCGGCGCTGTCCACAAAGGcccc    >  2:755219/1‑149 (MQ=255)
                                                                                                                                                aaCGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGcccc    <  2:209598/149‑1 (MQ=255)
                                                                                                                                                aaCGCCAGCTGTTATCATCAGCAATTGCCACTTTGCCCGGTTACATCATGGAACCCAGCACAGGCACCGTCACCACTTCCTCCCGCGCTATACGTTTTAGCAACGCGCTGCCGATGTCTGTATAAAGCTGCTCTTTCCACAAACGcccc    <  2:41583/149‑1 (MQ=255)
                                                                                                                                                  cGCCAGCTGTTATCCTCAGGAAAACCCACTTTGCCCGGTAACAACATGGAAGCCAGCCCAGGCACCGTCTCGACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGc  <  2:740934/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAACGCCCCGC  >  NC_000913/3689586‑3689880

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: ‑