BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	443,023	T→C	6.9%	intergenic (+26/+28)	thiI → / ← yajL	tRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	443,023	0	T	C	6.9%	597.2 / 22.3	187	intergenic (+26/+28)	thiI/yajL	tRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone
Reads supporting (aligned to +/- strand): ref base T (73/99); new base C (7/6); total (82/105)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 198 total aligned reads displayed.

TCAATGGCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCAGCGTCGGGAAGCCGGTCATATTACCAATCGGGAAG  >  NC_000913/442750‑443288

               ctaccgacgTGATCCTCGATAGCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGc                                                                                                                                                                                                                                                           <  2:283989/271‑1 (MQ=255)

                     actgtagcgtagATATCCGTTCTAACGTTGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTGTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATc                                                                                                                                                                                                                                                     <  4:297281/264‑1 (MQ=255)

                        tGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGa                                                                                                                                                                                                                                                  <  2:240327/275‑1 (MQ=255)

                                     ccGTTCGATGGATGAAAAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaat                                                                                                                                                                                                                                     <  4:98220/275‑1 (MQ=255)

                                      ctttttATTGATGAACAGGAAGGAAAGCGACGTAGAGGTGAAGGTGTTGATTGTGTTTCAAGGCCGTTCTATAAACTGAGCGCCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaat                                                                                                                                                                                                                                     <  4:254963/269‑1 (MQ=255)

                                                                                                                    aGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACCACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTAGTGATGCAACTTCACAGGCTTTTTCACGCCCAACCAACAGGCCGaaaa                                                                                                                                                       >  2:298829/1‑274 (MQ=255)

                                                                                                                             tttGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGt                                                                                                                                                                                                                                                                           <  8:240859/149‑1 (MQ=255)

                                                                                                                             tttGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGt                                                                                                                                                                                                                                                                           >  6:841187/1‑149 (MQ=255)

                                                                                                                             tttGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGt                                                                                                                                                                                                                                                                           >  7:1454792/1‑149 (MQ=255)

                                                                                                                              ttGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTc                                                                                                                                                                                                                                                                          <  7:1145262/149‑1 (MQ=255)

                                                                                                                                  cGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGCGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGAt                                                                                                                                                                                                                                                                      <  7:575609/149‑1 (MQ=255)

                                                                                                                                  cGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGAt                                                                                                                                                                                                                                                                      <  5:1099160/149‑1 (MQ=255)

                                                                                                                                   gATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATg                                                                                                                                                                                                                                                                     >  5:952643/1‑149 (MQ=255)

                                                                                                                                   gATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATg                                                                                                                                                                                                                                                                     >  6:1278481/1‑149 (MQ=255)

                                                                                                                                   gATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTc                                                                                                                                       >  2:111691/1‑275 (MQ=255)

                                                                                                                                      ctcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGa                                                                                                                                                                                                                                                                  <  7:620475/149‑1 (MQ=255)

                                                                                                                                       tcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTACTTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGAc                                                                                                                                                                                                                                                                 <  8:1378483/149‑1 (MQ=255)

                                                                                                                                       tcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGAc                                                                                                                                                                                                                                                                 >  8:4887/1‑149 (MQ=255)

                                                                                                                                           ccAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGttt                                                                                                                                                                                                                                                             <  6:944604/149‑1 (MQ=255)

                                                                                                                                             aGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGc                                                                                                                                                                                                                                                           >  6:580288/1‑149 (MQ=255)

                                                                                                                                              gAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGcc                                                                                                                                                                                                                                                          <  6:338800/149‑1 (MQ=255)

                                                                                                                                                aCAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGc                                                                                                                                                                                                                                                        >  5:908655/1‑149 (MQ=255)

                                                                                                                                                   aaaCCTGGCTGGNGTGGNGTGAGCGCGGGGTGATGAGCCGTNTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaat                                                                                                                                                                                                                                     <  4:346741/165‑1 (MQ=255)

                                                                                                                                                         ggCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATcc                                                                                                                                                                                                                                                    >  2:36162/1‑144 (MQ=255)

                                                                                                                                                         ggCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGacac                                                                                                                                                                                                                                               >  7:1062507/1‑149 (MQ=255)

                                                                                                                                                                gtggtgTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGt                                                                                                                                                                                                                                        >  7:953337/1‑149 (MQ=255)

                                                                                                                                                                  ggtgTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaa                                                                                                                                                                                                                                      <  5:515078/149‑1 (MQ=255)

                                                                                                                                                                  gctgtgAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaa                                                                                                                                                                                                                                      <  7:1189778/147‑1 (MQ=255)

                                                                                                                                                                   gtgtGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaat                                                                                                                                                                                                                                     <  8:700307/149‑1 (MQ=255)

                                                                                                                                                                   gtgtGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaat                                                                                                                                                                                                                                     >  5:812375/1‑149 (MQ=255)

                                                                                                                                                                   gtgtGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaat                                                                                                                                                                                                                                     >  8:334198/1‑149 (MQ=255)

                                                                                                                                                                     gtgAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaataa                                                                                                                                                                                                                                   >  5:421658/1‑149 (MQ=255)

                                                                                                                                                                         gcgcGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATtat                                                                                                                                                                                                                               >  8:270835/1‑149 (MQ=255)

                                                                                                                                                                               ggTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATcc                                                                                                                                                                                                                         <  6:1108905/149‑1 (MQ=255)

                                                                                                                                                                               ggTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATcc                                                                                                                                                                                                                         <  8:1285708/149‑1 (MQ=255)

                                                                                                                                                                                        ccGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATcc                                                                                                                                                                                                                         <  5:1056424/140‑1 (MQ=255)

                                                                                                                                                                                        ccGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATcc                                                                                                                                                                                                                         >  6:1056424/1‑140 (MQ=255)

                                                                                                                                                                                           tCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGCTAATTTTCAGACCCCAGTCGTACGCTTCACCCGGCCCCTCGCCGGTCAGCAATTTTACCCGTGCATCCCAcacg                                                                                >  4:74965/1‑274 (MQ=255)

                                                                                                                                                                                            cTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCAc                                                                                                                                                                                                            >  5:998637/1‑149 (MQ=255)

                                                                                                                                                                                              gCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCa                                                                                                                                                                                                          <  8:931026/149‑1 (MQ=255)

                                                                                                                                                                                                 ggCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGtt                                                                                                                                                                                                       <  7:865138/149‑1 (MQ=255)

                                                                                                                                                                                                 ggCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGtt                                                                                                                                                                                                       <  8:1454792/149‑1 (MQ=255)

                                                                                                                                                                                                     ctctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGa                                                                                                                                                                                                   <  7:270835/149‑1 (MQ=255)

                                                                                                                                                                                                     ctctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGa                                                                                                                                                                                                   <  8:911548/149‑1 (MQ=255)

                                                                                                                                                                                                      tctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGAt                                                                                                                                                                                                  <  6:812375/149‑1 (MQ=255)

                                                                                                                                                                                                      tctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGAt                                                                                                                                                                                                  <  8:987660/149‑1 (MQ=255)

                                                                                                                                                                                                      tctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATCATAAATCCCTGCCGCCATCACCAGTtgt                                                                                                                                                                                                     >  8:1254702/1‑146 (MQ=255)

                                                                                                                                                                                                      tctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATCATAAATCCCTGCCGCCATCACCAGTtgt                                                                                                                                                                                                     <  7:1254702/146‑1 (MQ=255)

                                                                                                                                                                                                        tATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGc                                                                                                                                                                                                <  5:1479796/149‑1 (MQ=255)

                                                                                                                                                                                                          tCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGa                                                                                                                                                                                                   <  7:645950/144‑1 (MQ=255)

                                                                                                                                                                                                          tCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCa                                                                                                                                                                                              >  7:77545/1‑149 (MQ=255)

                                                                                                                                                                                                                 gAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATg                                                                                                                                                                                       >  7:828105/1‑149 (MQ=255)

                                                                                                                                                                                                                  aGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATgg                                                                                                                                                                                      <  6:76912/149‑1 (MQ=255)

                                                                                                                                                                                                                     aGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtgt                                                                                                                                                                                                     <  6:385484/131‑1 (MQ=255)

                                                                                                                                                                                                                     aGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtgt                                                                                                                                                                                                     >  5:385484/1‑131 (MQ=255)

                                                                                                                                                                                                                     aGGGCTTTAAAAATTTGAAGGTTTATCGCCCTTAATTTTTGGTTTTTACGCTGCACCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTcggcg                                                     <  4:124510/275‑1 (MQ=255)

                                                                                                                                                                                                                        gCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaat                                                                                                                                                                                                                                     <  2:310258/96‑1 (MQ=255)

                                                                                                                                                                                                                        gCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCtttt                                                                                                                                                                                <  5:580288/149‑1 (MQ=255)

                                                                                                                                                                                                                         cTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCttttt                                                                                                                                                                               <  8:526040/149‑1 (MQ=255)

                                                                                                                                                                                                                             aaCAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACg                                                                                                                                                                           >  7:969274/1‑149 (MQ=255)

                                                                                                                                                                                                                                 aTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgccg                                                                                                                                                                                                                   >  2:16578/1‑105 (MQ=255)

                                                                                                                                                                                                                                  tgtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCccaa                                                                                                                                                                      >  8:943118/1‑149 (MQ=255)

                                                                                                                                                                                                                                   gtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCccaac                                                                                                                                                                     <  6:1034432/149‑1 (MQ=255)

                                                                                                                                                                                                                                     gAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCAt                                                                                                                                                                                                               <  4:16578/105‑1 (MQ=255)

                                                                                                                                                                                                                                         gTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAg                                                                                                                                                               <  7:712831/149‑1 (MQ=255)

                                                                                                                                                                                                                                         gTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAg                                                                                                                                                               <  7:1000358/149‑1 (MQ=255)

                                                                                                                                                                                                                                           atatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCTGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGt                                                                                                                                                             <  6:690348/149‑1 (MQ=255)

                                                                                                                                                                                                                                           atatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGt                                                                                                                                                             <  8:1088845/149‑1 (MQ=255)

                                                                                                                                                                                                                                            tatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTc                                                                                                                                                            >  6:462226/1‑149 (MQ=255)

                                                                                                                                                                                                                                             atCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCg                                                                                                                                                           >  8:609736/1‑149 (MQ=255)

                                                                                                                                                                                                                                              tCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCccaac                                                                                                                                                                     <  8:764696/138‑1 (MQ=255)

                                                                                                                                                                                                                                               cGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGAt                                                                                                                                                         <  8:659254/149‑1 (MQ=255)

                                                                                                                                                                                                                                                gCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCcaaccaac                                                                                                                                                                 <  5:136728/140‑1 (MQ=255)

                                                                                                                                                                                                                                                gCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCcaaccaac                                                                                                                                                                 >  6:136728/1‑140 (MQ=255)

                                                                                                                                                                                                                                                gCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATa                                                                                                                                                        >  7:436746/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  ccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt                                                                                                                                                                                   <  5:1243767/120‑1 (MQ=255)

                                                                                                                                                                                                                                                  ccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt                                                                                                                                                                                   >  6:1243767/1‑120 (MQ=255)

                                                                                                                                                                                                                                                   cGTACTTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCGGACACTACTCGTAATAATTATAAATCCCTGCAGCCATCACCAGTTGTGATGCCACTTCAGGGGCGTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCAGCGTCGGGAAGcc                        <  2:74965/274‑1 (MQ=255)

                                                                                                                                                                                                                                                   cGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACANGTCGATAATTTTCAGACCAAAGTCGAAAGCTGTACCCGGCCCCTGGCTGGTCATCAATTTTACCCGTTCATCCCNCACCACGCCCTTGTCCCGCCCTGTTTCGGCNNGGATTTTTTTTTTccgttgttctttat                         >  3:48929/1‑260 (MQ=255)

                                                                                                                                                                                                                                                     tAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAAtttt                                                                                                                                                   >  5:1423629/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       cttttTGGTTTTTACGTAGAATCTGGTCAGATGCGACGTTTGCCGCAGCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAAtt                                               <  1:48929/246‑1 (MQ=255)

                                                                                                                                                                                                                                                       aTTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCa                                                                                                                                                 <  8:105379/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       aTTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCa                                                                                                                                                 <  6:1164337/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         ttGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGa                                                                                                                                               <  8:953337/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              gTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaaa                                                                                                                                          <  6:132988/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      gTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAt                                                                                                                                                                                        >  5:295777/1‑95 (MQ=255)

                                                                                                                                                                                                                                                                      gTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAt                                                                                                                                                                                        <  6:295777/95‑1 (MQ=255)

                                                                                                                                                                                                                                                                      gTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAg                                                                                                                                  >  7:354864/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      gTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAg                                                                                                                                  >  7:113247/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       tCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGc                                                                                                                                 >  5:71858/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       tCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACAACTGGTGAt                                                                                                                                                                                                  <  6:1314248/84‑1 (MQ=255)

                                                                                                                                                                                                                                                                       tCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGa                                                                                                                                                          >  5:758974/1‑124 (MQ=255)

                                                                                                                                                                                                                                                                       tCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGa                                                                                                                                                          <  6:758974/124‑1 (MQ=255)

                                                                                                                                                                                                                                                                       tCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTCCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGc                                                                                                                                 >  8:1423336/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                            tCTGGTCAGATGCGACGTTTGTCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACGGGTCGATAATTTTCAGACCAAAGt                                                                                                                                        >  7:472275/1‑137 (MQ=255)

                                                                                                                                                                                                                                                                            tCTGGTCAGATGCGACGTTTGTCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGt                                                                                                                                        <  8:472275/137‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 tCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaa                                                                                                                                           <  6:167591/129‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 tCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaa                                                                                                                                           >  5:167591/1‑129 (MQ=255)

                                                                                                                                                                                                                                                                                 tCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGc                                                                                                                       >  8:174985/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                                 tCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGc                                                                                                                       >  5:1294245/1‑149 (MQ=255)

TCAATGGCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCAGCGTCGGGAAGCCGGTCATATTACCAATCGGGAAG  >  NC_000913/442750‑443288

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