Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 343,523 G→A A214T (GCA→ACA)  yahK → broad specificity NADPH‑dependent aldehyde reductase, Zn‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913343,5230GA98.4% 1312.2 / 5.0 383A214T (GCA→ACA) yahKbroad specificity NADPH‑dependent aldehyde reductase, Zn‑containing
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base A (190/186);  total (192/191)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 383 total aligned reads displayed.

TCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCA  >  NC_000913/343376‑343670
                                                                                                                                                   |                                                                                                                                                   
tCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCGCGCGATGGGGGCACATGTGGTGGCATTTAGCACTTCTGAGGCAAAACGCGAAGCGCCac                                                                                                                                                  >  6:499903/1‑150 (MQ=255)
  gccgcTACGTCACTGGCAGGCCGGGCCGTGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACaaaa                                                                                                                                                <  5:566614/151‑1 (MQ=255)
  gccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACaaaa                                                                                                                                                <  2:1582157/151‑1 (MQ=255)
     gcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGcc                                                                                                                                             >  5:597017/1‑151 (MQ=255)
        aCGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATTTGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCt                                                                                                                                           <  5:305420/150‑1 (MQ=255)
        aCGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCt                                                                                                                                           <  3:814686/150‑1 (MQ=255)
        aCGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGccc                                                                                                                                            <  6:1700385/149‑1 (MQ=255)
          gTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATGTGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTggg                                                                                                                                        <  3:1275934/151‑1 (MQ=255)
          gTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCTAACCGCGAAGCGACAAAAGCCCTggg                                                                                                                                        >  8:169801/1‑151 (MQ=255)
            cagtGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTgggg                                                                                                                                       <  8:418735/147‑1 (MQ=255)
            cACTTGCAGGCCGGGCCGTGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGTGATTAAGCTGGCCCACGCGATGTGTGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTgggg                                                                                                                                       <  1:1338014/150‑1 (MQ=255)
            cACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTgggg                                                                                                                                       <  3:1520409/150‑1 (MQ=255)
              cTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGTACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGcc                                                                                                                                    >  1:241086/1‑151 (MQ=255)
                ggCAGGCCGGGCCGGGTAAAATAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATTTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGcc                                                                                                                                    <  4:514700/149‑1 (MQ=255)
                 gCAGTCCGGGCCGTGTAAAAAAGTGGGCGTGGTGGGGATCGGCGGTCTGTGACATATGGGGATTAAGCTGGCCCACGCTATGGTGGCACATGTTGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGAt                                                                                                                                 <  2:651992/151‑1 (MQ=255)
                   aGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATg                                                                                                                                >  2:1126481/1‑150 (MQ=255)
                    ggccgggccgGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGCAGCGACAAAAGCCCTGGGGGCCGATGa                                                                                                                               >  5:670830/1‑150 (MQ=255)
                      ccgggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAgt                                                                                                                            >  5:383792/1‑151 (MQ=255)
                      ccgggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAgt                                                                                                                            <  8:260128/151‑1 (MQ=255)
                      ccgggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAgt                                                                                                                            >  2:762170/1‑151 (MQ=255)
                          gccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCGGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAgttgtt                                                                                                                        <  7:262183/151‑1 (MQ=255)
                            cgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTa                                                                                                                       >  8:1452977/1‑150 (MQ=255)
                                 aaaaaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCa                                                                                                                  <  6:1423867/150‑1 (MQ=255)
                                 aaaaaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCAc                                                                                                                 <  7:1339105/151‑1 (MQ=255)
                                 aaaaaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCAc                                                                                                                 <  1:730377/151‑1 (MQ=255)
                                  aaaaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACg                                                                                                                >  1:1007094/1‑151 (MQ=255)
                                        tGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAACAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGc                                                                                                          >  3:231227/1‑151 (MQ=255)
                                         gggCGGGGTCGGCATCGGCGGTGTGGGAGATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGc                                                                                                          <  3:1359430/150‑1 (MQ=255)
                                            cgggggcGGCAGCGGCGGTCTGGGACATATGGGGCTTAAGCTGGCCCACGCGAGGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGAt                                                                                                      <  6:866058/145‑1 (MQ=255)
                                            cGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGAt                                                                                                      <  3:109925/151‑1 (MQ=255)
                                            cGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGAt                                                                                                      >  2:1178164/1‑151 (MQ=255)
                                             gTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATg                                                                                                     >  5:440875/1‑151 (MQ=255)
                                              tGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATg                                                                                                     >  6:370666/1‑150 (MQ=255)
                                               ggTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATg                                                                                                     >  3:512331/1‑149 (MQ=255)
                                                 tCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATg                                                                                                     >  3:1398245/1‑147 (MQ=255)
                                                  cGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATg                                                                                                <  3:1419020/151‑1 (MQ=255)
                                                  cGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATg                                                                                                >  7:705241/1‑151 (MQ=255)
                                                  cGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATg                                                                                                <  4:1525611/151‑1 (MQ=255)
                                                   ggCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATgg                                                                                               >  7:562645/1‑151 (MQ=255)
                                                   ggCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATgg                                                                                               >  8:749156/1‑151 (MQ=255)
                                                     cATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATggcg                                                                                             <  3:592842/151‑1 (MQ=255)
                                                       tCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATggcggc                                                                                           >  7:20589/1‑151 (MQ=255)
                                                        cggcggTCTGGTAGATATGGGGATTAAGGTGTCCCAAGGGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCTAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATggcggc                                                                                           <  6:752819/150‑1 (MQ=255)
                                                         ggcggTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTc                                                                                         >  5:160652/1‑151 (MQ=255)
                                                         ggcggTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGGGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTc                                                                                         <  1:1651388/151‑1 (MQ=255)
                                                           cggTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCAt                                                                                       >  1:1446148/1‑151 (MQ=255)
                                                               cTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTg                                                                                    >  5:1114789/1‑150 (MQ=255)
                                                                 gggACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGaa                                                                                  <  5:955320/150‑1 (MQ=255)
                                                                 gggACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGa                                                                                   >  2:248899/1‑149 (MQ=255)
                                                                  ggACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAaga                                                                                <  4:1653173/151‑1 (MQ=255)
                                                                  ggACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAaga                                                                                <  4:248342/151‑1 (MQ=255)
                                                                   gACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAag                                                                                 >  5:904357/1‑149 (MQ=255)
                                                                    aCATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGt                                                                              >  7:1223296/1‑151 (MQ=255)
                                                                     cATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGTCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGtt                                                                             <  5:187342/151‑1 (MQ=255)
                                                                     cATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAACGCCCGGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATCAGATGGCGGCTCCTCTGAAGACtt                                                                             >  4:73619/1‑151 (MQ=255)
                                                                     cATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGtt                                                                             <  2:1369959/151‑1 (MQ=255)
                                                                        atGGGTATTATGCTGGCCCACGCTATGGGGGCAGATGTGGTGTCATTTTCCACTTCTGAGTCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATTAAGTTGTTAACTCACTCAATGCCGATGAGATGGCGGCTCATCTTAAGAGTTTc                                                                           <  1:564575/150‑1 (MQ=255)
                                                                        atGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTc                                                                           <  5:1250862/150‑1 (MQ=255)
                                                                        atGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCg                                                                          >  7:300137/1‑151 (MQ=255)
                                                                              aTTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCCCGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTc                                                                     >  2:908637/1‑150 (MQ=255)
                                                                              aTTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTc                                                                     >  1:629942/1‑150 (MQ=255)
                                                                               ttAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCAt                                                                   >  7:984925/1‑151 (MQ=255)
                                                                                tAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCGCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGGTGATCTTAAGAGTTTCGATTTCAtt                                                                  >  6:159176/1‑151 (MQ=255)
                                                                                tAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCAtt                                                                  >  8:1576743/1‑151 (MQ=255)
                                                                                  aGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCAtttt                                                                >  6:984552/1‑151 (MQ=255)
                                                                                  aGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCAtttt                                                                <  6:589633/151‑1 (MQ=255)
                                                                                    cTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTc                                                                                         >  7:1104341/1‑124 (MQ=255)
                                                                                      ggCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGa                                                              <  5:1551177/149‑1 (MQ=255)
                                                                                       gCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAAt                                                            >  4:1414711/1‑150 (MQ=255)
                                                                                       gCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCAAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAAt                                                            >  7:149739/1‑150 (MQ=255)
                                                                                        cccACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATAc                                                          <  8:645791/151‑1 (MQ=255)
                                                                                           aCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGt                                                       <  5:1349917/151‑1 (MQ=255)
                                                                                           aCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGt                                                       >  1:1214016/1‑151 (MQ=255)
                                                                                            cgcgATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTa                                                      <  4:1070695/151‑1 (MQ=255)
                                                                                             gcgATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAg                                                     <  4:1126554/151‑1 (MQ=255)
                                                                                               gATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCt                                                   >  6:1069269/1‑151 (MQ=255)
                                                                                                 tGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCt                                                                                     >  3:662808/1‑115 (MQ=255)
                                                                                                  gggggCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTgc                                                 <  5:1006974/150‑1 (MQ=255)
                                                                                                     ggCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCca                                             <  1:946440/151‑1 (MQ=255)
                                                                                                     ggCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCca                                             <  1:384722/151‑1 (MQ=255)
                                                                                                        acaTGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACAt                                          <  8:1218874/151‑1 (MQ=255)
                                                                                                        acaTGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACAt                                          <  3:117958/151‑1 (MQ=255)
                                                                                                        acaTGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACAt                                          >  2:564929/1‑151 (MQ=255)
                                                                                                         caTGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACAt                                          >  6:1225282/1‑150 (MQ=255)
                                                                                                          aTGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATaa                                        <  5:763049/151‑1 (MQ=255)
                                                                                                              ggtggCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtctc                                    <  1:1125580/151‑1 (MQ=255)
                                                                                                              ggtggCATTTACCACTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtctc                                    >  7:892287/1‑151 (MQ=255)
                                                                                                                    aTTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAt                              >  8:480904/1‑151 (MQ=255)
                                                                                                                           aCTTCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacc                        >  2:1012763/1‑150 (MQ=255)
                                                                                                                              tCTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCt                    >  6:1081857/1‑151 (MQ=255)
                                                                                                                               cTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCtt                   >  7:79431/1‑151 (MQ=255)
                                                                                                                               cTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCtt                   >  7:457265/1‑151 (MQ=255)
                                                                                                                               cTGAGGCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCtt                   >  3:1236818/1‑151 (MQ=255)
                                                                                                                                    gCAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGGTGCGCCACATAATCTCGCCGAGTTTACCACCTTGCTg               >  4:849624/1‑150 (MQ=255)
                                                                                                                                     cAAAACGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGa              >  6:416520/1‑150 (MQ=255)
                                                                                                                                      aaaaCGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGa              >  2:1532497/1‑149 (MQ=255)
                                                                                                                                        aaCGCGAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCg          >  2:344027/1‑151 (MQ=255)
                                                                                                                                          cgcgAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTg        >  4:1314949/1‑151 (MQ=255)
                                                                                                                                           gcgAAGCGACAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACAACCTTGCTGAAGCGTg        >  8:1083931/1‑150 (MQ=255)
                                                                                                                                             gAAGCGACAAAAGCGCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATg     <  4:881636/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCA  >  NC_000913/343376‑343670

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