Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 231,803 C→T P228S (CCG→TCG)  yafD → endo/exonuclease/phosphatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,8030CT96.9% 490.1 / NA 127P228S (CCG→TCG) yafDendo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base T (52/71);  total (56/71)
Fisher's exact test for biased strand distribution p-value = 3.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

CACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAG  >  NC_000913/231659‑231948
                                                                                                                                                |                                                                                                                                                 
caccacAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGctcgctcg                                                                                                                                             <  4:588079/151‑1 (MQ=255)
  ccacAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGa                                                                                                                                            <  1:1118245/150‑1 (MQ=255)
  ccacAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGAt                                                                                                                                           <  4:158722/151‑1 (MQ=255)
  ccacAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGAGCCGTCGGGCGTTTGGTCGCTCGCTCGAt                                                                                                                                           <  1:1153634/151‑1 (MQ=255)
     caGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGAtttt                                                                                                                                        <  2:294893/151‑1 (MQ=255)
        cGGCCCGGTCATTATGGCGGGGGGTTTTAATGCCCGGGGGCGGAGAAGGGGGAACGCGTTTTATTGCTTTTTGCGGGGAATGTCGCTGCGCCCGGGGCGGTTTTCCGGTGATCAGCGCCGGCGGGGGGTTTGTCGCCCGCCCGGTtttttt                                                                                                                                     >  1:1822354/1‑151 (MQ=255)
           cccGGTCATTATGGCGGGAGATTTTAAAGCCTTGGGGCGGAGAAGGAAGAACGCGGTTTATCGCCTTTCGCGGGGAAAGTCGCTGCGCCCCGGGCGTTTTTTCGGTGATTAGCGGCGTCCGGGGGTTGGGCGCTCGCTCGgttttgtttt                                                                                                                                   >  1:806596/1‑150 (MQ=255)
           cccGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTc                                                                                                                                  <  4:1469693/151‑1 (MQ=255)
           cccGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTc                                                                                                                                  <  1:1408770/151‑1 (MQ=255)
             cGGTCATTATGGCGGGAGATTTTAATGCCTGGGGCCGTAGAAGGAGGAACGCGTTATATCGCTTTTCGCGGGGAATGTCGCTGCGCCCGGGGGGGTTTACCGATGATCAACGCCGGCGGGCGGTTGGTCGCTCGCTCGGTtttttttttt                                                                                                                                 >  1:1661821/1‑148 (MQ=255)
             cGGTCATTATGGCGGGAGATTTTAAAGCCTGGGGGCCGAGAAGGATGAAAACGTTATAATGCTTTGCGCGGGAAAGGTCGCCGCGCCAGGGGGGGTTTTTTGAAGATTAGCGCCGTCGGGCGGTTGGTCGCCCGCTCGGTtttttttttt                                                                                                                                 >  3:191372/1‑148 (MQ=255)
             cGGTCATTATGGCGGGAGATTTTAAAGCCTGGGGACGGGGAAGGGTGAACGCGTGATAACGCTTTGCGCGGGGCATGTCGCTGCCCCCGGGGCGGTTTTCCGCTGGTCAACGCCGCTGGGGGTGTGGTCGCCCGCTccggtttttttt                                                                                                                                   >  3:280929/1‑137 (MQ=255)
             cGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTa                                                                                                                                <  2:342086/151‑1 (MQ=255)
                 cATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCg                                                                                                                             <  4:167064/150‑1 (MQ=255)
                  aTTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTg                                                                                                                           <  1:63068/151‑1 (MQ=255)
                       ggCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTg                                                                                                                      >  1:398066/1‑151 (MQ=255)
                         cGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGa                                                                                                                     >  4:1188611/1‑150 (MQ=255)
                          gggAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAAc                                                                                                                   <  3:489122/151‑1 (MQ=255)
                            gagaTTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGt                                                                                                                 <  1:479100/151‑1 (MQ=255)
                              gaTTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCa                                                                                                               <  1:853046/151‑1 (MQ=255)
                                 ttCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGt                                                                                                             <  4:547791/150‑1 (MQ=255)
                                 ttCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGt                                                                                                             <  4:450160/150‑1 (MQ=255)
                                   cAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGaa                                                                                                          <  4:1792720/151‑1 (MQ=255)
                                   cAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGa                                                                                                           <  1:1916984/150‑1 (MQ=255)
                                   cAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGa                                                                                                           >  1:91015/1‑150 (MQ=255)
                                     aTGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGc                                                                                                        >  3:210385/1‑151 (MQ=255)
                                                   aGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGtt                                                                                           >  2:296773/1‑150 (MQ=255)
                                                    gAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTAc                                                                                         <  4:1138821/151‑1 (MQ=255)
                                                    gAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTAc                                                                                         >  4:302019/1‑151 (MQ=255)
                                                     aaGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcg                                                                                        >  2:48298/1‑151 (MQ=255)
                                                        gATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcgcgc                                                                                     <  3:308961/151‑1 (MQ=255)
                                                         aTGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcgcgcg                                                                                    >  4:1496657/1‑151 (MQ=255)
                                                           gAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACgcgcgc                                                                                   <  4:221880/150‑1 (MQ=255)
                                                           gAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCt                                                                                  >  4:172041/1‑151 (MQ=255)
                                                             aCGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTc                                                                                <  4:1759323/151‑1 (MQ=255)
                                                               gcgTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCg                                                                              <  2:407876/151‑1 (MQ=255)
                                                                cgTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGa                                                                             <  4:111487/151‑1 (MQ=255)
                                                                    atatCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATcac                                                                         >  4:587760/1‑151 (MQ=255)
                                                                    atatCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATcac                                                                         >  4:1353662/1‑151 (MQ=255)
                                                                       tCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAAt                                                                      >  4:1638519/1‑151 (MQ=255)
                                                                         gCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATcc                                                                    <  1:858929/151‑1 (MQ=255)
                                                                         gCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATcc                                                                    <  4:1901640/151‑1 (MQ=255)
                                                                          cTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCg                                                                   >  4:98422/1‑151 (MQ=255)
                                                                           tttGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGc                                                                  >  1:531971/1‑151 (MQ=255)
                                                                            ttGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCt                                                                 >  3:494881/1‑151 (MQ=255)
                                                                             tGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCt                                                                 >  1:1230641/1‑150 (MQ=255)
                                                                             tGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCt                                                                 <  2:90353/150‑1 (MQ=255)
                                                                              gcgcGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGATa                                                                <  1:202296/150‑1 (MQ=255)
                                                                                gcgGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTc                                                             >  3:372405/1‑151 (MQ=255)
                                                                                       aTGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAAt                                                      <  2:372416/151‑1 (MQ=255)
                                                                                       aTGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAAt                                                      <  1:1433309/151‑1 (MQ=255)
                                                                                         gTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGATTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAAt                                                      >  1:1379524/1‑149 (MQ=255)
                                                                                          tCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTc                                                    >  4:262835/1‑150 (MQ=255)
                                                                                           cGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAg                                                  <  1:1356743/151‑1 (MQ=255)
                                                                                           agCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTc                                                    >  4:70394/2‑149 (MQ=255)
                                                                                            gCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGt                                                 >  4:329655/1‑151 (MQ=255)
                                                                                            gCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGt                                                 <  4:191808/151‑1 (MQ=255)
                                                                                               gcgcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTccc                                              >  1:375819/1‑151 (MQ=255)
                                                                                               gcgcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTccc                                              <  1:1861179/151‑1 (MQ=255)
                                                                                               gcgcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTc                                                <  4:1238948/149‑1 (MQ=255)
                                                                                                 gcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCgg                                            >  1:1068426/1‑151 (MQ=255)
                                                                                                 gcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCgg                                            <  1:354214/151‑1 (MQ=255)
                                                                                                   cAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCa                                          <  4:205506/151‑1 (MQ=255)
                                                                                                    aGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTCGAATTCAGTCCCGGCaa                                         <  1:122572/151‑1 (MQ=255)
                                                                                                     ggTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCa                                          >  4:492669/1‑149 (MQ=255)
                                                                                                     ggTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCa                                          >  4:1245065/1‑149 (MQ=255)
                                                                                                        gCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCt                                     <  1:1918006/151‑1 (MQ=255)
                                                                                                        gCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCt                                     >  1:997774/1‑151 (MQ=255)
                                                                                                          gTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGa                                   <  4:1795800/151‑1 (MQ=255)
                                                                                                          gTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGa                                   <  1:1000762/151‑1 (MQ=255)
                                                                                                             ttACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGata                                 <  2:400284/150‑1 (MQ=255)
                                                                                                              tACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataa                                >  3:34678/1‑150 (MQ=255)
                                                                                                               aCCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataaa                               <  4:936058/150‑1 (MQ=255)
                                                                                                                 cGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataaataa                            <  4:529132/151‑1 (MQ=255)
                                                                                                                  gatgatCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataaata                             >  4:823891/1‑149 (MQ=255)
                                                                                                                     gatCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTa                        <  3:216139/151‑1 (MQ=255)
                                                                                                                      atCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTAt                       <  1:43018/151‑1 (MQ=255)
                                                                                                                       tCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGGAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATg                      >  3:251045/1‑151 (MQ=255)
                                                                                                                        cAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATg                      <  1:1012247/150‑1 (MQ=255)
                                                                                                                        cAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATg                      <  1:367793/150‑1 (MQ=255)
                                                                                                                        cAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGt                     >  1:1554873/1‑151 (MQ=255)
                                                                                                                          gcgcCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCa                   <  1:1353110/151‑1 (MQ=255)
                                                                                                                            gcCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAgg                 >  1:134092/1‑151 (MQ=255)
                                                                                                                            gcCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAgg                 <  1:411247/151‑1 (MQ=255)
                                                                                                                            gcCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAgg                 <  1:181656/151‑1 (MQ=255)
                                                                                                                             ccGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGt                <  4:1823069/151‑1 (MQ=255)
                                                                                                                             ccGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGt                <  1:1342502/151‑1 (MQ=255)
                                                                                                                             ccGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGt                <  1:171991/151‑1 (MQ=255)
                                                                                                                                tCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTg             >  1:1184002/1‑151 (MQ=255)
                                                                                                                                 ggggCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGc            >  1:1045514/2‑151 (MQ=255)
                                                                                                                                 cGGGCGTTTGGTCGCTCGCACGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGc            <  2:46558/151‑1 (MQ=255)
                                                                                                                                   ggCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGcc           >  1:1268826/1‑150 (MQ=255)
                                                                                                                                   ggCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGcc           >  1:587246/1‑150 (MQ=255)
                                                                                                                                   ggCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCca          <  4:940502/151‑1 (MQ=255)
                                                                                                                                    gCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCcac         >  4:80578/1‑151 (MQ=255)
                                                                                                                                      gTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCcac         >  1:480669/1‑149 (MQ=255)
                                                                                                                                      gTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCcac         >  1:308355/1‑149 (MQ=255)
                                                                                                                                        ttGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAgg      <  4:605889/150‑1 (MQ=255)
                                                                                                                                          ggTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCa   <  1:1678752/151‑1 (MQ=255)
                                                                                                                                           gTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAg  <  1:47805/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAG  >  NC_000913/231659‑231948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑