Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,609,406 T→C 5.9% F307S (TTC→TCC)  hyfG → hydrogenase 4, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,609,4060TC 5.9% 212.6 / 2.9 68F307S (TTC→TCC) hyfGhydrogenase 4, subunit
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TTGAAAATGCACTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCG  >  NC_000913/2609259‑2609552
                                                                                                                                                   |                                                                                                                                                  
ttGAAAATGCACTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTc                                                                                                                                                   >  3:111681/1‑149 (MQ=255)
     aaTGCACTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTCCcgcac                                                                                                                                              >  1:359143/1‑147 (MQ=255)
      aTGCACTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTg                                                                                                                                             >  7:178580/1‑149 (MQ=255)
          aCTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCgg                                                                                                                                          >  8:199680/1‑148 (MQ=255)
              gggATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGaa                                                                                                                                       <  4:3322/147‑1 (MQ=255)
                gATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTCCCGCACGCGGAaaaatt                                                                                                                                   <  2:359143/149‑3 (MQ=255)
                      ggTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGa                                                                                                                             <  7:181567/149‑1 (MQ=255)
                       gTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGGCTAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGAc                                                                                                                            <  6:302024/149‑1 (MQ=255)
                       gTGCCGAAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGAc                                                                                                                            <  8:170117/149‑1 (MQ=255)
                             cAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGc                                                                                                                      >  5:239816/1‑149 (MQ=255)
                                cGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCgg                                                                                                                    >  8:354332/1‑148 (MQ=255)
                                   gCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCcgc                                                                                                                                                >  4:123411/1‑117 (MQ=255)
                                   gCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCcgc                                                                                                                                                <  3:123411/117‑1 (MQ=255)
                                   gCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAAtt                                                                                                                >  6:345766/1‑149 (MQ=255)
                                    cacaTACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTg                                                                                                               <  1:308525/149‑1 (MQ=255)
                                        tACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAAt                                                                                                                 <  7:10769/143‑1 (MQ=255)
                                        tACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAAt                                                                                                                 >  8:10769/1‑143 (MQ=255)
                                             ttCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCATGCCATTTCGTTTGTTTTGATACAGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCgg                                                                                                       <  5:319853/148‑1 (MQ=255)
                                                gCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTcg                                                                                                   >  1:351898/1‑149 (MQ=255)
                                                gCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTcg                                                                                                   >  7:125326/1‑149 (MQ=255)
                                                gCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTc                                                                                                    <  4:56933/148‑1 (MQ=255)
                                                 cTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCTGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTcgc                                                                                                  <  7:199680/149‑1 (MQ=255)
                                                    gATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTcgc                                                                                                  >  7:173899/1‑146 (MQ=255)
                                                    gATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTcgc                                                                                                  <  8:173899/146‑1 (MQ=255)
                                                      ttCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTCCcgcac                                                                                                                                              <  7:296426/100‑3 (MQ=255)
                                                      ttCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTCCcgcac                                                                                                                                              >  8:296426/1‑98 (MQ=255)
                                                         tgctgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAAcc                                                                                          >  1:7393/1‑149 (MQ=255)
                                                         tgctgGAAGTCGAACGGCTACACAGTCATTTGCTTAACATTTGCCTCTCCTGCCATTTCGTTGGTTTTGATAACGGCTTTATGCAATTTTTCCGCGTGAGGGAAAAGTCGATGACGATGGCGGAATTGCTGATAGGGTAGCGTAAAAc                                                                                           >  4:203165/1‑148 (MQ=255)
                                                           ctgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTa                                                                                        <  4:111681/149‑1 (MQ=255)
                                                            tgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTAc                                                                                       >  2:88006/1‑149 (MQ=255)
                                                             tgAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCGTCTCCTGCCTTTTCGTTGGGTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACg                                                                                      <  6:49520/148‑1 (MQ=255)
                                                               aaGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGt                                                                                    <  6:268071/149‑1 (MQ=255)
                                                                 gTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCt                                                                                  <  7:350239/149‑1 (MQ=255)
                                                                    gAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGa                                                                                <  8:178580/148‑1 (MQ=255)
                                                                     aaCGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGaa                                                                               >  8:2035/1‑148 (MQ=255)
                                                                            acacAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGAGGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATtg                                                                       >  2:155951/1‑149 (MQ=255)
                                                                            acacAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATtg                                                                       >  2:40524/1‑149 (MQ=255)
                                                                               caGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATtggtg                                                                    <  7:354332/149‑1 (MQ=255)
                                                                                aGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCgg                                                                                                                                          >  1:283098/1‑78 (MQ=255)
                                                                                aGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTAATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATttgtgg                                                                   >  3:192630/1‑149 (MQ=255)
                                                                                  tCATTTGATGAACCTTGGCCTCTCCTGCCATTTCGTAGGATTAGATACCCTCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGCTGGCCGAATAGCTGATCGGGTCCCGTAAAACCTACGGTCTGGATCGGATTGgtggtg                                                                 >  4:57079/1‑149 (MQ=255)
                                                                                     tttGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTc                                                              <  1:226924/149‑1 (MQ=255)
                                                                                         cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCTGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGTTGGTGTTcgcc                                                           <  6:239816/148‑1 (MQ=255)
                                                                                         cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCt                                                                            >  6:321071/1‑131 (MQ=255)
                                                                                         cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCt                                                                            <  5:321071/131‑1 (MQ=255)
                                                                                         cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgccg                                                          >  8:325776/1‑149 (MQ=255)
                                                                                         cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgccg                                                          >  7:142531/1‑149 (MQ=255)
                                                                                          ttAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAGAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgccgc                                                         >  2:91096/1‑149 (MQ=255)
                                                                                            aaCCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgccgc                                                         <  2:351898/147‑1 (MQ=255)
                                                                                                   gCCTCTCCTGCCATTTCGTTGGTTTTGCTACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAATTCGATGACGATGGCGTATTTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATtctc                                                <  2:71380/149‑1 (MQ=255)
                                                                                                   gCCTCTCCTGCCATTTCGTTGGTGTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATtctc                                                <  8:390227/149‑1 (MQ=255)
                                                                                                     ctctCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACTGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCaa                                              <  8:125326/149‑1 (MQ=255)
                                                                                                         ccTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAagag                                          >  3:349552/1‑149 (MQ=255)
                                                                                                           tGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGc                                         >  6:93791/1‑148 (MQ=255)
                                                                                                             ccATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAAc                                      <  5:345766/149‑1 (MQ=255)
                                                                                                              cATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGAtt                                                                        >  6:287653/1‑114 (MQ=255)
                                                                                                              cATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGAtt                                                                        <  5:287653/114‑1 (MQ=255)
                                                                                                              cATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACg                                     >  3:168486/1‑149 (MQ=255)
                                                                                                                  tCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCTGGTAGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCt                                  <  3:203165/148‑1 (MQ=255)
                                                                                                                  tCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCt                                  <  1:68238/148‑1 (MQ=255)
                                                                                                                   cGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTg                                 <  1:40524/148‑1 (MQ=255)
                                                                                                                    gTTGGTTTTGATACCGGCTTTATGCAATTTATCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGc                                <  1:88006/148‑1 (MQ=255)
                                                                                                                             gATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGa                       <  8:142531/148‑1 (MQ=255)
                                                                                                                                   ggCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGt                <  7:325776/149‑1 (MQ=255)
                                                                                                                                       ttATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGGCTGCAAACGCTGAAACTGGTgcgc            >  2:219376/1‑149 (MQ=255)
                                                                                                                                           gCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTgcgc            >  3:57137/1‑145 (MQ=255)
                                                                                                                                            cAATTTTTCCGCGTGCGGGAAGAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGAt       <  1:91096/149‑1 (MQ=255)
                                                                                                                                            cAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGAt       >  3:323751/1‑149 (MQ=255)
                                                                                                                                                  ttCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgcgcg  <  1:219376/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGAAAATGCACTGGGGATTGAGGTGCCGCAACGAGCACATACTATTCGCTCGATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCG  >  NC_000913/2609259‑2609552

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: