BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	858,927	C→A	5.3%	T45T (ACG→ACT)	ybiU ←	DUF1479 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	858,927	0	C	A	5.3%	460.4 / 10.5	129	T45T (ACG→ACT)	ybiU	DUF1479 family protein
Reads supporting (aligned to +/- strand): ref base C (60/62); new base A (5/2); total (65/64)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

Only 100 of 132 total aligned reads displayed.

GCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCAT  >  NC_000913/858786‑859074

gCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCa                                                                                                                                               <  4:211366/148‑1 (MQ=255)

 ccTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCaa                                                                                                                                              <  4:759112/148‑1 (MQ=255)

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

        tCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATc                                                                                                                                       >  4:834990/1‑148 (MQ=255)

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

         cACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGc                                                                                                                                     >  1:18224/1‑149 (MQ=255)

           ccGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATg                                                                                                                                            >  3:207663/1‑140 (MQ=255)

           ccGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATg                                                                                                                                            <  4:207663/140‑1 (MQ=255)

           ccGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCt                                                                                                                                    >  4:716589/1‑148 (MQ=255)

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

             gCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTa                                                                                                                                  >  3:388616/1‑148 (MQ=255)

                      cgACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTaa                                                                                                                         <  4:334261/148‑1 (MQ=255)

                      cgACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTaa                                                                                                                         <  3:361337/148‑1 (MQ=255)

                          gTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGAt                                                                                                                     <  3:717448/148‑1 (MQ=255)

                          gTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATc                                                                                                                    <  2:846401/149‑1 (MQ=255)

                           tttAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATct                                                                                                                   <  1:617227/149‑1 (MQ=255)

                            ttAATCTGTTTGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATct                                                                                                                   >  3:617928/1‑148 (MQ=255)

                             tAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctg                                                                                                                  <  3:130732/148‑1 (MQ=255)

                               aTCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATctgc                                                                                                                 >  2:799984/1‑147 (MQ=255)

                                  tGTTCGCGCTGGTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGAc                                                                                                            >  1:798277/1‑149 (MQ=255)

                                   gTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGAc                                                                                                            >  3:932557/1‑148 (MQ=255)

                                    ttCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGt                                                                                                          >  2:135613/1‑149 (MQ=255)

                                             ctctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCt                                                                                                 >  2:1046015/1‑149 (MQ=255)

                                              tctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTg                                                                                                >  2:1010277/1‑149 (MQ=255)

                                               ctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTgc                                                                                               <  1:847937/149‑1 (MQ=255)

                                               ctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTgc                                                                                               <  2:1012330/149‑1 (MQ=255)

                                               ctGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTg                                                                                                <  4:1028203/148‑1 (MQ=255)

                                                     tAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGc                                                                                         >  1:135499/1‑149 (MQ=255)

                                                       aCATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCCGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAg                                                                                       >  1:836053/1‑149 (MQ=255)

                                                            atcTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGCCGTCGCCAAGCTGCGCCCGCAGCGCGt                                                                                  <  2:458760/147‑1 (MQ=255)

                                                            aCCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGcgcg                                                                                   >  3:1037701/1‑148 (MQ=255)

                                                            aCCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGt                                                                                  <  1:784323/149‑1 (MQ=255)

                                                              cTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGc                                                                                >  1:474042/1‑149 (MQ=255)

                                                              cTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGc                                                                                >  2:862580/1‑149 (MQ=255)

                                                                  tttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTc                                                                             >  4:366801/1‑148 (MQ=255)

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

                                                                   ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                            >  4:293432/1‑148 (MQ=255)

                                                                   ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                            >  4:44386/1‑148 (MQ=255)

                                                                   ttGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCa                                                                            >  4:575351/1‑148 (MQ=255)

                                                                    tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCAt                                                                           <  4:617928/148‑1 (MQ=255)

                                                                    tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATc                                                                          <  2:133772/149‑1 (MQ=255)

                                                                    tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATc                                                                          <  2:881584/149‑1 (MQ=255)

                                                                    tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATc                                                                          <  2:1046117/149‑1 (MQ=255)

                                                                    tGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATc                                                                          <  2:798277/149‑1 (MQ=255)

                                                                        atCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGAc                                                                      <  2:434525/149‑1 (MQ=255)

                                                                         tCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGAc                                                                      <  3:786369/148‑1 (MQ=255)

                                                                         tCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACg                                                                     >  2:546041/1‑149 (MQ=255)

                                                                           ggCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGa                                                                   <  2:18224/149‑1 (MQ=255)

                                                                            gCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGAt                                                                  >  1:134063/1‑149 (MQ=255)

                                                                             cATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGAt                                                                  <  3:575351/148‑1 (MQ=255)

                                                                             cATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGAt                                                                  <  3:716589/148‑1 (MQ=255)

                                                                             cATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATa                                                                 >  2:532822/1‑149 (MQ=255)

                                                                             aaTAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATa                                                                 >  2:912091/2‑149 (MQ=255)

                                                                                aaGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCCGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgct                                                              <  2:836053/149‑1 (MQ=255)

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

                                                                                  gACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCAGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAgctgc                                                            >  2:497758/1‑149 (MQ=255)

                                                                                     aGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGa                                                                   <  2:418877/139‑1 (MQ=255)

                                                                                     aGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGa                                                                   >  1:418877/1‑139 (MQ=255)

                                                                                     aGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCtt                                                          <  4:360001/148‑1 (MQ=255)

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

                                                                                         ccGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGa                                                     <  1:927099/149‑1 (MQ=255)

                                                                                            gCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCg                                                  <  2:741153/149‑1 (MQ=255)

                                                                                               aGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggca                                               <  1:122771/149‑1 (MQ=255)

                                                                                               aGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggca                                               <  1:270991/149‑1 (MQ=255)

                                                                                               aGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggc                                                <  4:853096/148‑1 (MQ=255)

                                                                                                gACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggca                                               >  4:1025429/1‑148 (MQ=255)

                                                                                                  cGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCggcagg                                             <  4:388616/148‑1 (MQ=255)

                                                                                                     cATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCaa                                          <  3:834990/148‑1 (MQ=255)

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

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

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

                                                                                                           ccTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCg                                   >  1:892467/1‑149 (MQ=255)

                                                                                                             tgtgCTTTGAGTGCGTTGATTTCAGCCACTCGCGTAGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGc                                  >  4:1000250/1‑148 (MQ=255)

                                                                                                              gtgCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCt                                 >  4:487291/1‑148 (MQ=255)

                                                                                                                  tttGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaa                            >  1:61577/1‑149 (MQ=255)

                                                                                                                       gTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTaaaa                          >  1:1028145/1‑146 (MQ=255)

                                                                                                                           gTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAg                   >  1:848321/1‑149 (MQ=255)

                                                                                                                           gTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGGGTCGCTGGTAAAAGTAGaa                    >  3:839960/1‑148 (MQ=255)

                                                                                                                            ttGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGa                                                                       <  4:181834/96‑1 (MQ=255)

                                                                                                                            ttGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGa                                                                       >  3:181834/1‑96 (MQ=255)

                                                                                                                            ttGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAg                   >  3:347509/1‑148 (MQ=255)

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

                                                                                                                             tGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGAt                                                    <  3:58677/114‑1 (MQ=255)

                                                                                                                             tGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGAt                                                    >  4:58677/1‑114 (MQ=255)

                                                                                                                             tGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGc                  >  4:909526/1‑148 (MQ=255)

                                                                                                                              gATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCa                <  2:61577/149‑1 (MQ=255)

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

                                                                                                                                  tCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac            >  2:63303/1‑149 (MQ=255)

                                                                                                                                  attgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac            >  1:284506/4‑149 (MQ=255)

                                                                                                                                   ttgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTc                                                                             >  4:253415/3‑83 (MQ=255)

                                                                                                                                   ttgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac            >  4:255919/3‑148 (MQ=255)

                                                                                                                                   ttgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac            >  3:405209/3‑148 (MQ=255)

                                                                                                                                   cAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAac            >  3:227546/1‑148 (MQ=255)

                                                                                                                                    tgCCACGCGAGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGAt                                                    <  3:963927/106‑1 (MQ=255)

                                                                                                                                    tgCCACGCGAGTGGCAAAGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGAt                                                    >  4:963927/2‑107 (MQ=255)

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

                                                                                                                                    aGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAacc           >  3:148724/1‑148 (MQ=255)

                                                                                                                                           cgcgTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCtctc    <  3:293432/148‑1 (MQ=255)

                                                                                                                                           cgcgTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTCGTAAAAGTACAAGCCATAACCATCTCTCa   <  2:892467/149‑1 (MQ=255)

                                                                                                                                             cgTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCAt  >  3:341077/1‑148 (MQ=255)

GCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTGGCAATGTCATCGCTTAGCTGATTAAAGATCTGCTGGACGTCGCCAAGCTGCGCCCGCAGCGCGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCAT  >  NC_000913/858786‑859074

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