Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 231,589 C→T 5.5% H156H (CAC→CAT)  yafD → endo/exonuclease/phosphatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,5890CT 5.5% 233.6 / 2.4 73H156H (CAC→CAT) yafDendo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (42/27);  new base T (2/2);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

CTTCTGGCGTAATGACCCTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTT  >  NC_000913/231441‑231736
                                                                                                                                                    |                                                                                                                                                   
cTTCTGGCGTAATGACCCTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATacac                                                                                                                                                     <  4:274004/149‑1 (MQ=255)
 ttCTGGCGTAATGACCCTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACg                                                                                                                                                    >  5:47071/1‑149 (MQ=255)
     ggCGTAATGACCCTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGt                                                                                                                                                >  5:357621/1‑149 (MQ=255)
     ggCGTAATGACCCTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGt                                                                                                                                                >  8:20691/1‑149 (MQ=255)
        gTAATGACCCTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCaa                                                                                                                                             >  4:20327/1‑149 (MQ=255)
         tAATGACCCTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAAc                                                                                                                                            >  4:264033/1‑149 (MQ=255)
                ccTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTTACGGTCTATACATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTc                                                                                                                                     >  1:335118/1‑149 (MQ=255)
                 cTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCt                                                                                                                                    >  7:214614/1‑149 (MQ=255)
                  tttCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTg                                                                                                                                   >  3:95538/1‑149 (MQ=255)
                     cggcggCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTggg                                                                                                                                 <  2:6213/148‑1 (MQ=255)
                      ggcggcACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGc                                                                                                                                <  1:82080/148‑1 (MQ=255)
                      ggcggcACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGc                                                                                                                                <  7:217874/148‑1 (MQ=255)
                         ggcACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTgg                                                                                                                            <  7:20691/149‑1 (MQ=255)
                          gcACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGa                                                                                                                           <  7:25410/149‑1 (MQ=255)
                          gcACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGa                                                                                                                           <  3:307451/149‑1 (MQ=255)
                             caTCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGt                                                                                                                        >  3:267324/1‑149 (MQ=255)
                                ccAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCta                                                                                                                     >  2:2904/1‑149 (MQ=255)
                                 cAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCta                                                                                                                     <  4:44763/148‑1 (MQ=255)
                                  aGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCtata                                                                                                                   <  3:264033/149‑1 (MQ=255)
                                     gtATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGt                                                                                                                 >  5:19126/1‑148 (MQ=255)
                                      tATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACATGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGt                                                                                                                 >  4:277095/1‑147 (MQ=255)
                                      tATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACATGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTaa                                                                                                               >  4:281755/1‑149 (MQ=255)
                                          gctgcCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCa                                                                                                            <  1:2904/148‑1 (MQ=255)
                                           ctgcCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAg                                                                                                           <  1:69238/148‑1 (MQ=255)
                                                cGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGGTACtt                                                                                                     <  5:357135/149‑1 (MQ=255)
                                                      gTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCG‑CTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTg                                                                                              >  1:335812/1‑149 (MQ=255)
                                                         aaCGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGc                                                                                            <  8:214614/149‑1 (MQ=255)
                                                            gagaACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTTCCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCATCTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGAt                                                                                         <  8:288464/149‑1 (MQ=255)
                                                                 cccATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACTCCGTCAACTTCAGTCTGGTCGTGGATGTCTATAGTAAGCAGTTAATTCCTATTGGCGATCAGa                                                                                     >  6:63796/1‑148 (MQ=255)
                                                                         gCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTcac                                                                             <  6:357621/148‑1 (MQ=255)
                                                                         gCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGAAGTTACTTCCTATTGGCGATCAGATAGCTcac                                                                             >  7:274207/1‑148 (MQ=255)
                                                                          cGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTcacca                                                                           <  3:20327/149‑1 (MQ=255)
                                                                           gTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTcacca                                                                           >  4:224875/1‑148 (MQ=255)
                                                                                gCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAg                                                                        <  6:19126/146‑1 (MQ=255)
                                                                                gCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCgg                                                                     >  8:195825/1‑149 (MQ=255)
                                                                                gCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGGTAATATACACGCCGTCAACTTCAGTCTGGGCGGGGATGTCTATAGTAAGAAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCgg                                                                     >  7:106332/1‑149 (MQ=255)
                                                                                  gAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACATGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGcc                                                                   <  3:281755/149‑1 (MQ=255)
                                                                                  gAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACATGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGcc                                                                   <  3:277095/149‑1 (MQ=255)
                                                                                   aaGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGccc                                                                  >  5:14073/1‑149 (MQ=255)
                                                                                     gTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCgg                                                                <  4:267324/149‑1 (MQ=255)
                                                                                      tCGGCATTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGt                                                               >  4:161629/1‑149 (MQ=255)
                                                                                        ggCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGa                                                                                                                           <  4:83582/87‑1 (MQ=255)
                                                                                        ggCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGa                                                                                                                           >  3:83582/1‑87 (MQ=255)
                                                                                         gCACTGGTGACGGTATATCCATTACCTGACACCCGCATGTTGATGGGGGTTAATATACACGCCGTCACCGTCAGTCTGGGCGTGGATGTAGATAGTAAGCGGTTACTTACTATTGTAGATCAGATAGCTCACCACATCGGCCCGGTCAt                                                            >  5:197755/1‑149 (MQ=255)
                                                                                            ttGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTAt                                                         <  3:161629/148‑1 (MQ=255)
                                                                                               gTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGAt                                                                                    <  3:167315/119‑1 (MQ=255)
                                                                                               gTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGAt                                                                                    >  4:167315/1‑119 (MQ=255)
                                                                                                tGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCg                                                     >  1:272479/1‑149 (MQ=255)
                                                                                                          tCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTc                                            <  6:47071/148‑1 (MQ=255)
                                                                                                           ccATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCac                                          >  8:102623/1‑148 (MQ=255)
                                                                                                            cATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGc                                                                                            <  7:341271/98‑1 (MQ=255)
                                                                                                            cATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGc                                                                                            >  8:341271/1‑98 (MQ=255)
                                                                                                             aTTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCa                                                                                       <  1:91053/102‑1 (MQ=255)
                                                                                                             aTTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCa                                                                                       >  2:91053/1‑102 (MQ=255)
                                                                                                                  cTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACcaca                                                                         <  4:67931/111‑1 (MQ=255)
                                                                                                                  cTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACcaca                                                                         >  3:67931/1‑111 (MQ=255)
                                                                                                                  cTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCAt                                                            <  4:222823/124‑1 (MQ=255)
                                                                                                                  cTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCAt                                                            >  3:222823/1‑124 (MQ=255)
                                                                                                                   tGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTg                                    >  8:13022/1‑147 (MQ=255)
                                                                                                                    gACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGgaga                                                <  1:244832/134‑1 (MQ=255)
                                                                                                                    gACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGgaga                                                >  2:244832/1‑134 (MQ=255)
                                                                                                                    gACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAg                                 >  7:36349/1‑149 (MQ=255)
                                                                                                                    gACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAg                                 >  7:341351/1‑149 (MQ=255)
                                                                                                                       acCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCg                              >  6:170500/1‑149 (MQ=255)
                                                                                                                       acCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCg                              >  3:267965/1‑149 (MQ=255)
                                                                                                                           gCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAg                           >  2:78662/1‑148 (MQ=255)
                                                                                                                            ccTGTTGATGGTGGTTAATATACACGCAGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTACGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGaa                         <  2:335118/149‑1 (MQ=255)
                                                                                                                              tGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAg                        >  6:343142/1‑148 (MQ=255)
                                                                                                                                 tGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGAt                     >  1:206621/1‑148 (MQ=255)
                                                                                                                                     ggtggtTAATATACACGCCGTCAACTTCAGTCTGGGCGGGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTTTGGCGGGAGATTTCATTGCCTGGAGCCGTAGAAGGATGAAcg                >  4:279216/1‑149 (MQ=255)
                                                                                                                                       tggtTAATATCCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgcg              >  6:102992/1‑149 (MQ=255)
                                                                                                                                         gtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATAAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGt             >  3:5432/1‑148 (MQ=255)
                                                                                                                                          ttAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTta           >  4:281279/1‑149 (MQ=255)
                                                                                                                                              tataCACGCCGTCAACTTCAGTCTGTGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGCCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGCCGGGAGATTTCAATGACTGGAGCCGAAGAAGGATGAACGCTATATATc       >  4:199462/1‑149 (MQ=255)
                                                                                                                                                   acGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCttt  >  6:170064/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTTCTGGCGTAATGACCCTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTT  >  NC_000913/231441‑231736

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