Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,689,724 G→C 6.9% P180A (CCC→GCC)  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,7240GC 6.9% 226.4 / 9.9 87P180A (CCC→GCC) bcsZendo‑1,4‑D‑glucanase
Reads supporting (aligned to +/- strand):  ref base G (33/47);  new base C (4/2);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

CTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAA  >  NC_000913/3689577‑3689872
                                                                                                                                                   |                                                                                                                                                    
cTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGCTTGGTTTCGCGCTGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTgggg                                                                                                                                                     <  1:911032/149‑1 (MQ=255)
     gcTTGCTCATAGCGCACCCCGGCTGGCGAAAAACCGTTCGGCGCGCGTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGGCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTaaa                                                                                                                                                <  1:103926/148‑1 (MQ=255)
       tttCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAAc                                                                                                                                              <  1:537961/149‑1 (MQ=255)
       tttCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAAc                                                                                                                                              >  2:364923/1‑149 (MQ=255)
         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:43978/149‑1 (MQ=255)
                                                                           tttCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTc                                                                          <  1:846052/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                          <  1:206116/149‑1 (MQ=255)
                                                                                                                           ccAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCg                          <  2:658360/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:815595/149‑1 (MQ=255)
                                                                                                                                            tAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCttt         <  1:348460/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)
                                                                                                                                                   |                                                                                                                                                    
CTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACGGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCCCGCGCGATACGTTTTAGCAACGCGCTGCCGATGTCGGTATAACGCTGCTCTTTCCACAAA  >  NC_000913/3689577‑3689872

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