BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	858,924	T→G	5.5%	R46R (CGA→CGC)	ybiU ←	DUF1479 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	858,924	0	T	G	5.5%	605.7 / 20.9	174	R46R (CGA→CGC)	ybiU	DUF1479 family protein
Reads supporting (aligned to +/- strand): ref base T (88/76); new base G (8/2); total (96/78)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 174 total aligned reads displayed.

AAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACC  >  NC_000913/858781‑859068

aaaTGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGt                                                                                                                                             <  2:1478311/149‑1 (MQ=255)

 aaTGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTg                                                                                                                                            >  1:1348177/1‑149 (MQ=255)

  aTGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTgg                                                                                                                                           >  2:1198569/1‑149 (MQ=255)

  aTGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGATTTGAGTGCGTTGATTTCAGCCACTCGCGTgg                                                                                                                                           >  2:1117515/1‑149 (MQ=255)

      ccTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCGGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAAt                                                                                                                                       >  2:1368785/1‑149 (MQ=255)

      ccTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAAt                                                                                                                                       <  4:449599/149‑1 (MQ=255)

      ccTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAAt                                                                                                                                       <  2:1410564/149‑1 (MQ=255)

        ttttATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGt                                                                                                                                     >  4:1262769/1‑149 (MQ=255)

            aTCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATc                                                                                                                                 >  2:155420/1‑149 (MQ=255)

            aTCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATc                                                                                                                                 <  3:1308079/149‑1 (MQ=255)

            aTCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATc                                                                                                                                 >  2:607207/1‑149 (MQ=255)

            aTCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATc                                                                                                                                 <  2:27000/149‑1 (MQ=255)

             tCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCg                                                                                                                                >  3:692904/1‑149 (MQ=255)

               aCCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCt                                                                                                                              >  4:582657/1‑149 (MQ=255)

               aCCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCt                                                                                                                              >  1:1409290/1‑149 (MQ=255)

                 cGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTa                                                                                                                            <  4:31255/149‑1 (MQ=255)

                    acaACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCt                                                                                                                         >  2:262639/1‑149 (MQ=255)

                    acaACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCt                                                                                                                         >  3:969458/1‑149 (MQ=255)

                      aaCCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGa                                                                                                                       <  3:1251327/149‑1 (MQ=255)

                      aaCCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGa                                                                                                                       >  1:43305/1‑149 (MQ=255)

                      aaCCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGa                                                                                                                       <  3:281163/149‑1 (MQ=255)

                       aCCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGAt                                                                                                                      <  1:475430/149‑1 (MQ=255)

                       aCCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGAt                                                                                                                      <  3:1683686/149‑1 (MQ=255)

                          gcgACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTaa                                                                                                                   >  1:1659179/1‑149 (MQ=255)

                           cgACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTaaa                                                                                                                  <  4:773401/149‑1 (MQ=255)

                            gACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGt                                                                                                                                             >  2:1074770/1‑121 (MQ=255)

                            gACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGt                                                                                                                                             <  1:1074770/121‑1 (MQ=255)

                              cGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGAt                                                                                                               <  3:464581/149‑1 (MQ=255)

                               gTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATc                                                                                                              >  1:739887/1‑149 (MQ=255)

                                tttAATCTGTTCGCGCTGCTCTACAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATct                                                                                                             >  2:568640/1‑149 (MQ=255)

                                   aaTCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgct                                                                                                          >  3:1438716/1‑149 (MQ=255)

                                    aTCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgctg                                                                                                         >  1:844437/1‑149 (MQ=255)

                                    aTCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgctg                                                                                                         >  4:871800/1‑149 (MQ=255)

                                    aTCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgc                                                                                                           >  3:1503092/1‑147 (MQ=255)

                                    aTCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgc                                                                                                           <  4:1503092/147‑1 (MQ=255)

                                     tCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTgg                                                                                                                                           <  4:1487887/114‑1 (MQ=255)

                                     tCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTgg                                                                                                        >  2:988215/1‑149 (MQ=255)

                                     tCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTgg                                                                                                        >  1:728351/1‑149 (MQ=255)

                                     tCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTgg                                                                                                        >  1:211020/1‑149 (MQ=255)

                                     tCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTAATTTCAGCCACTCGCGTgg                                                                                                                                           >  3:1487887/1‑114 (MQ=255)

                                              gctgctCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgct                                                                                                          <  4:1160466/138‑1 (MQ=255)

                                              gctgctCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCa                                                                                               >  4:1534021/1‑149 (MQ=255)

                                                   tctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTg                                                                                          >  3:1601924/1‑149 (MQ=255)

                                                      gCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTgcgc                                                                                       >  3:974755/1‑149 (MQ=255)

                                                          tAACATGACCTGCCTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGc                                                                                   >  2:1599282/1‑149 (MQ=255)

                                                           aaCATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCa                                                                                  <  4:204027/149‑1 (MQ=255)

                                                            aCATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCAt                                                                                                                                  <  4:802830/100‑1 (MQ=255)

                                                            aCATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCAt                                                                                                                                  >  3:802830/1‑100 (MQ=255)

                                                                gACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGcgcg                                                                             >  1:1640149/1‑149 (MQ=255)

                                                                 aCCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTaaa                                                                                                                  >  3:763814/1‑111 (MQ=255)

                                                                 aCCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGt                                                                                                    >  3:831820/1‑125 (MQ=255)

                                                                  ccTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTg                                                                           >  2:619789/1‑149 (MQ=255)

                                                                  ccTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTg                                                                           <  1:821047/149‑1 (MQ=255)

                                                                  ccTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTg                                                                           >  3:1104204/1‑149 (MQ=255)

                                                                   cTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTgc                                                                                         >  1:31752/1‑134 (MQ=255)

                                                                   cTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGc                                                                          <  1:86152/149‑1 (MQ=255)

                                                                       tttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                      >  1:995568/1‑149 (MQ=255)

                                                                            tatCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGa                                                                 <  2:1649840/149‑1 (MQ=255)

                                                                             atCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGAc                                                                >  4:1621611/1‑149 (MQ=255)

                                                                              tCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACg                                                               <  4:1318226/149‑1 (MQ=255)

                                                                                ggCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGAc                                                                >  2:339134/1‑146 (MQ=255)

                                                                                ggCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGAc                                                                <  1:339134/146‑1 (MQ=255)

                                                                                     aaGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgct                                                        >  2:871308/1‑149 (MQ=255)

                                                                                     aaGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgct                                                        <  1:154989/149‑1 (MQ=255)

                                                                                       gACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgctgc                                                      >  4:1470505/1‑149 (MQ=255)

                                                                                       gACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgctgc                                                      >  3:1218078/1‑149 (MQ=255)

                                                                                        aCAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGCCGGATAgctgct                                                     >  4:1148196/1‑149 (MQ=255)

                                                                                          aGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCttt                                                   <  4:1360248/149‑1 (MQ=255)

                                                                                            cACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTtgt                                                 <  3:1621611/149‑1 (MQ=255)

                                                                                            cACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTtgt                                                 <  2:1560446/149‑1 (MQ=255)

                                                                                              ccGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCTATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGa                                               <  3:1629448/149‑1 (MQ=255)

                                                                                                ggCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATc                                             <  1:35959/149‑1 (MQ=255)

                                                                                                         gCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCaa                                    <  3:1470505/149‑1 (MQ=255)

                                                                                                          cATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAAt                                   <  3:1534021/149‑1 (MQ=255)

                                                                                                          cATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAAt                                   >  2:1207727/1‑149 (MQ=255)

                                                                                                          cATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAAt                                   <  1:1117515/149‑1 (MQ=255)

                                                                                                           atggCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCTTCGCTTAGCTGATTAAAGACCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCATGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAAtg                                  >  4:1074037/4‑149 (MQ=255)

                                                                                                           aTCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAAtg                                  >  1:1016950/1‑149 (MQ=255)

                                                                                                              gCCCTGTGCTTTGAGTGCGTTGGTTTCAGCCACTCGCGTGGCAATGTCAt                                                                                                                                  >  1:1452336/1‑50 (MQ=255)

                                                                                                                  tgtgCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCt                           <  2:405376/149‑1 (MQ=255)

                                                                                                                   gtgCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTg                          <  1:1025329/149‑1 (MQ=255)

                                                                                                                   gtgCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTg                          >  2:296174/1‑149 (MQ=255)

                                                                                                                   gtgCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTg                          <  1:568640/149‑1 (MQ=255)

                                                                                                                         tGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTTCTGGATGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCAGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaaaa                    <  4:755112/149‑1 (MQ=255)

                                                                                                                         tGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaaaa                    <  2:879349/149‑1 (MQ=255)

                                                                                                                         tGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaaaa                    <  2:43305/149‑1 (MQ=255)

                                                                                                                         tGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaaaa                    <  2:209265/149‑1 (MQ=255)

                                                                                                                            gTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTa                 >  1:1692778/1‑149 (MQ=255)

                                                                                                                                gTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTgg                                                                                                        <  1:521344/58‑1 (MQ=255)

                                                                                                                                 ttGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGc            >  3:531490/1‑149 (MQ=255)

                                                                                                                                     tttCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATa        >  2:671101/1‑149 (MQ=255)

                                                                                                                                      ttCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATaa       <  4:969458/149‑1 (MQ=255)

                                                                                                                                      ttCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATaa       <  1:1163935/149‑1 (MQ=255)

                                                                                                                                       attgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTgg                         <  4:162669/127‑1 (MQ=255)

                                                                                                                                       attgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGaa              <  1:264302/138‑1 (MQ=255)

                                                                                                                                       attgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac      >  2:1451923/4‑149 (MQ=255)

                                                                                                                                        ttgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAacc     >  1:1685983/3‑149 (MQ=255)

                                                                                                                                         aGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAacca    >  1:639908/1‑149 (MQ=255)

                                                                                                                                          gCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCt                                                                                           >  4:475246/1‑61 (MQ=255)

                                                                                                                                           ccACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAaccacc  >  2:240907/1‑149 (MQ=255)

AAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACC  >  NC_000913/858781‑859068

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