BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,396,058	Δ2 bp	coding (200‑201/3801 nt)	yhdP ←	DUF3971‑AsmA2 domains protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,396,057	0	G	.	98.7%	2136.3 / 20.5	459	coding (202/3801 nt)	yhdP	DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand): ref base G (3/3); new base . (233/220); total (236/223)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
*	NC_000913	3,396,058	0	C	.	98.5%	2125.6 / 24.0	458	coding (201/3801 nt)	yhdP	DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand): ref base C (5/2); new base . (231/220); total (236/222)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 491 total aligned reads displayed.

GGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAAGCGTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGCAGGCCGCTAACCAGCAGGGC  >  NC_000913/3395908‑3396201

||

ggTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAAgt                                                                                                                                                 >  7:1376797/1‑150 (MQ=255)

 gTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAAgtc                                                                                                                                                <  4:512717/151‑3 (MQ=255)

   gCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAAgtcgg                                                                                                                                              <  6:1113080/151‑5 (MQ=255)

    cGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAAggcggg                                                                                                                                             <  2:777108/151‑6 (MQ=255)

      aaaGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAAgtcgggcc                                                                                                                                           <  3:1656984/151‑8 (MQ=255)

          cgcAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCaaaa                                                                                                                                     <  6:1444962/151‑1 (MQ=255)

              gCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCt                                                                                                                                 <  2:586641/151‑1 (MQ=255)

              gCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCt                                                                                                                                 >  6:1233700/1‑151 (MQ=255)

                tGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGc                                                                                                                               >  8:769676/1‑151 (MQ=255)

                  ccAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAAGCGTCGGGCCAAAATTCTGc                                                                                                                               >  1:1411961/1‑151 (MQ=255)

                  ccAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGcc                                                                                                                              >  2:375479/1‑150 (MQ=255)

                  ccAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGcc                                                                                                                              >  1:1685534/1‑150 (MQ=255)

                   cAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGcc                                                                                                                              >  6:866581/1‑149 (MQ=255)

                     gAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAgct                                                                                                                          >  8:417095/1‑151 (MQ=255)

                      aaaGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAgctg                                                                                                                         >  2:1823787/1‑151 (MQ=255)

                         gTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAgctggct                                                                                                                      <  1:526611/151‑1 (MQ=255)

                           gAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCCTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTCCCAgctggctg                                                                                                                     >  7:843873/1‑150 (MQ=255)

                               gcgcgAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGc                                                                                                                <  8:217192/150‑1 (MQ=255)

                                cgcgAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCt                                                                                                               >  4:1379584/1‑151 (MQ=255)

                                cgcgAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCt                                                                                                               <  7:1627789/151‑1 (MQ=255)

                                 gcgAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTg                                                                                                              >  6:1845534/1‑151 (MQ=255)

                                 gcgAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTg                                                                                                              >  5:542004/1‑151 (MQ=255)

                                 gcgAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTg                                                                                                              >  4:244194/1‑151 (MQ=255)

                                  cgAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTg                                                                                                              >  8:599969/1‑150 (MQ=255)

                                   gAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGAc                                                                                                            >  1:1004268/1‑151 (MQ=255)

                                    aaaCTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACt                                                                                                           <  6:962778/151‑1 (MQ=255)

                                     aaCTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTg                                                                                                          <  2:484857/151‑1 (MQ=255)

                                     aaCTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTg                                                                                                          >  4:1092702/1‑151 (MQ=255)

                                       cTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAAATCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGCCTGAGAGCTGACTGGc                                                                                                        >  1:1025784/1‑151 (MQ=255)

                                        tGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCt                                                                                                       >  4:465787/1‑151 (MQ=255)

                                          ccAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTCAGAGCTGACTGGCTTc                                                                                                     >  6:765153/1‑151 (MQ=255)

                                           cAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCt                                                                                                    <  7:433467/151‑1 (MQ=255)

                                            aGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTa                                                                                                   <  5:1564697/151‑1 (MQ=255)

                                              cgcATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTAcc                                                                                                 >  1:1150571/1‑151 (MQ=255)

                                                cATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCgg                                                                                               <  3:106062/151‑1 (MQ=255)

                                                  tatGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGc                                                                                              >  6:313419/1‑150 (MQ=255)

                                                       aaCAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGcc                                                                                         >  5:995714/1‑150 (MQ=255)

                                                         cAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGt                                                                                      >  3:6326/1‑151 (MQ=255)

                                                          aGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTTCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTc                                                                                     <  5:1357107/151‑1 (MQ=255)

                                                          aGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTc                                                                                     <  7:1884762/151‑1 (MQ=255)

                                                           ggCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTcg                                                                                    <  4:130387/151‑1 (MQ=255)

                                                              tctGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCgg                                                                                 >  4:1734489/1‑151 (MQ=255)

                                                              gctGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCgg                                                                                 <  5:286910/150‑1 (MQ=255)

                                                                 gCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGAt                                                                               <  1:1858145/150‑1 (MQ=255)

                                                                 gCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGAt                                                                               <  1:898803/150‑1 (MQ=255)

                                                                   cAGACATCCACCCCCAGAGTAACGCGTATAACCGAAAATTCGCCGCCATCTTTTAGTACTGCACGGAACTCGTGTGCTACTA‑‑GTCGGGCCAAAATTCTGCAAGCTGGGAACGAGCTCAATCGCTTCTACCCGCATGCCAGTCGCGGAATCt                                                                            >  8:1160510/1‑151 (MQ=255)

                                                                    aGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCt                                                                            >  2:1322283/1‑150 (MQ=255)

                                                                    aGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTGCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTa                                                                           <  2:1017567/151‑1 (MQ=255)

                                                                     gACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTAt                                                                          <  1:154449/151‑1 (MQ=255)

                                                                     gACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTAt                                                                          >  5:1762267/1‑151 (MQ=255)

                                                                      aCATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTAtt                                                                         <  4:941690/151‑1 (MQ=255)

                                                                        aTCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTAtttt                                                                       <  7:1287097/151‑1 (MQ=255)

                                                                           cAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTttgtt                                                                    <  8:1799826/151‑1 (MQ=255)

                                                                            aGCGCCAGAGTAACGCGTTTAACCGAAAAATCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTttgttg                                                                   <  5:679383/151‑1 (MQ=255)

                                                                                 cAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCCTGCCAGTCGCGGGGTCTATTTTGTTGAGGAt                                                              >  8:1531382/1‑151 (MQ=255)

                                                                                  agagTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGAtt                                                             >  8:1351689/1‑151 (MQ=255)

                                                                                    agTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTc                                                           <  2:1184057/151‑1 (MQ=255)

                                                                                     gTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTcc                                                          <  3:1736481/151‑1 (MQ=255)

                                                                                         cgcgTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGAc                                                      >  5:1006938/1‑151 (MQ=255)

                                                                                         cgcgTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGAc                                                      <  7:766635/151‑1 (MQ=255)

                                                                                         cgcgTGTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGCGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCGGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGAc                                                      <  8:306390/151‑1 (MQ=255)

                                                                                            gTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTGTATTTTTTTGAGGATTTCCGGACGcc                                                   >  7:722191/1‑151 (MQ=255)

                                                                                            gTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGcc                                                   >  5:1475251/1‑151 (MQ=255)

                                                                                            gTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGcc                                                   <  1:364439/151‑1 (MQ=255)

                                                                                              ttAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCa                                                  <  2:1713688/150‑1 (MQ=255)

                                                                                              ttAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCa                                                  <  4:1872084/150‑1 (MQ=255)

                                                                                               tAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAgg                                                >  1:666478/1‑151 (MQ=255)

                                                                                                aaCCGAAACTTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAgg                                                <  3:582006/150‑1 (MQ=255)

                                                                                                aaCCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAgg                                                <  5:1266208/150‑1 (MQ=255)

                                                                                                  ccGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGt                                             <  3:1856316/151‑1 (MQ=255)

                                                                                                    gAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCa                                           <  4:876288/151‑1 (MQ=255)

                                                                                                       aaTTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCa                                           >  2:841498/1‑148 (MQ=255)

                                                                                                        aTTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCa                                           >  2:1497081/1‑147 (MQ=255)

                                                                                                          tCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCg                                                         <  5:1643704/131‑1 (MQ=255)

                                                                                                          tCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGATGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGAAATCTATTTTTTTGAGGATTTCCGGACGCCACGCGtgcacatgc                                     >  3:835048/1‑143 (MQ=255)

                                                                                                           cgccgcCATCTTTTAGTTCTGCACGGATGTCGTGGGCTTCAA‑‑GTCGGCCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCg                                    <  3:372159/151‑1 (MQ=255)

                                                                                                              cgccATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGt                                  >  7:516173/1‑150 (MQ=255)

                                                                                                              cgccATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTa                                 <  4:697282/151‑1 (MQ=255)

                                                                                                               gccATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTaa                                <  6:449963/151‑1 (MQ=255)

                                                                                                               gccATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTa                                 >  1:1192861/1‑150 (MQ=255)

                                                                                                                ccATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑C‑TCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAACTTGCGGTaa                                >  5:379509/1‑150 (MQ=255)

                                                                                                                 cATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTaa                                >  2:1841068/1‑149 (MQ=255)

                                                                                                                 cATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAg                              <  5:861061/151‑1 (MQ=255)

                                                                                                                    cTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCa                            >  7:1770881/1‑150 (MQ=255)

                                                                                                                    cTTTTAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCa                            >  6:1566154/1‑150 (MQ=255)

                                                                                                                     ttttAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAAt                          >  8:1725003/1‑151 (MQ=255)

                                                                                                                     ttttAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGAGGCGGTAAAGCAAt                          >  7:1159071/1‑151 (MQ=255)

                                                                                                                      tttAGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAAt                          >  1:752225/1‑150 (MQ=255)

                                                                                                                         aGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGc                      >  3:1259630/1‑151 (MQ=255)

                                                                                                                         aGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGc                      <  4:591658/151‑1 (MQ=255)

                                                                                                                         aGTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGc                      >  7:1554740/1‑151 (MQ=255)

                                                                                                                          gTTCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTACGGACGCCAGGCGTCAAGATGCCGTAAAGCAATACGc                      >  4:1584942/1‑150 (MQ=255)

                                                                                                                            tCTGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGCAg                    <  3:1444931/150‑1 (MQ=255)

                                                                                                                              tGCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGCAGGcc                 >  6:388875/1‑151 (MQ=255)

                                                                                                                               gCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGCAGGcc                 <  2:651617/150‑1 (MQ=255)

                                                                                                                               gCACGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGCAGGCCg                >  5:1144928/1‑151 (MQ=255)

                                                                                                                                 aCGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGCAGGCCGCt              <  6:1755209/151‑1 (MQ=255)

                                                                                                                                  cGGATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGGCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGCAGGCCGCTa             <  2:90145/151‑1 (MQ=255)

                                                                                                                                   ggATGTCGTGTGCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGAGTGGCTTCTACCCGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGCAGGCCGCTaa            <  6:64113/151‑1 (MQ=255)

                                                                                                                                           tgtgCTTCAA‑‑GTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGCAGGCCGCTAACCAGCAgg    >  6:1310108/1‑151 (MQ=255)

||

GGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCATCTTTTAGTTCTGCACGGATGTCGTGTGCTTCAAGCGTCGGGCCAAAATTCTGCCAGCTGGCTGAGAGCTGACTGGCTTCTACCGGCATGCCAGTCGCGGATTCTATTTTGTTGAGGATTTCCGGACGCCAGGCGTCAAGATGCGGTAAAGCAATACGCAGGCCGCTAACCAGCAGGGC  >  NC_000913/3395908‑3396201

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