BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	443,022	G→A	6.6%	intergenic (+25/+29)	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,022	0	G	A	6.6%	616.9 / 22.3	192	intergenic (+25/+29)	thiI/yajL	tRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone
Reads supporting (aligned to +/- strand): ref base G (77/102); new base A (7/6); total (84/108)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

Only 100 of 192 total aligned reads displayed.

TCAATGGCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCAGCGTCGGGAAGCCGGTCATATTACCAATCGGGAAG  >  NC_000913/442750‑443288

tCAATGGGTTCGGCCCGAACGACGGGATCACCGCTCGTCGTTCTATCGATGAACAGGACGATAAGCCACCGCAAGTCGAAGGGATGGATGTGGGTTCTCGTCCGTTCTATAAACTGAGCACCAAATGTGTCGATCTCGCCCAGAACAAACCCTGGCTGCTGGGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGt                                                                                                                                                                                                                                                                           <  4:280793/274‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                                              ttGGCGATCTCGCCCAGAACAAAACCTGGCTGCTGTGGTGTTAGCGCGGGGTGATCAGCCGCCTGAACGCGCTCTATCTGCACGAGCTGGGCTTTTACAATGAGAAAGTATCTGGCCCGAAATTTGAGGTTGTTACGTCTGATTTGGTc                                                                                                                                                                                                                                                                          <  6:1282935/149‑1 (MQ=255)

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

                                                                                                                                  cGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGAt                                                                                                                                                                                                                                                                      <  7:636800/149‑1 (MQ=255)

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

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

                                                                                                                                    aTCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGc                                                                                                                                                                                                                                                                    >  5:1134484/1‑149 (MQ=255)

                                                                                                                                      ctcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGa                                                                                                                                                                                                                                                                  <  6:1352611/149‑1 (MQ=255)

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

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

                                                                                                                                       tcGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGAc                                                                                                                                                                                                                                                                 >  5:1087946/1‑149 (MQ=255)

                                                                                                                                          aCCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGtt                                                                                                                                                                                                                                                              >  5:469865/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)

                                                                                                                                               aaCAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCg                                                                                                                                                                                                                                                         >  5:1030246/1‑149 (MQ=255)

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

                                                                                                                                                   aaaCCTGGCTGCTGTGGTGTGAGCGGGTGGTTGTGGGCCGGCTGCAGGAGCTCTATCTGCGCGAGGAGGGTTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATTCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCgg                                                                                                                        <  4:88314/274‑1 (MQ=255)

                                                                                                                                                   aaaCCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATc                                                                                                                                                                                                                                                     >  5:126909/1‑149 (MQ=255)

                                                                                                                                                      ccTGGCTCCTGTGGTGTCAGTACGGGAATACGAGTCGTCTGAAGGCGCTCTATCTGCGGGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCTTTACTTTGTGGTTTACCCGTCGCATCTGGTCAGATGCGACGTTTGCGGGATccta                                                                                                                                                                                                                                                  <  6:123647/149‑3 (MQ=255)

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

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

                                                                                                                                                             gctgTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGaca                                                                                                                                                                                                                                                <  4:36162/144‑1 (MQ=255)

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

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

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

                                                                                                                                                                     gtgAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACCCTACTCGtaataa                                                                                                                                                                                                                                   >  7:877716/1‑149 (MQ=255)

                                                                                                                                                                     gtgAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaataa                                                                                                                                                                                                                                   <  6:1134484/149‑1 (MQ=255)

                                                                                                                                                                            cgGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATaaa                                                                                                                                                                                                                            <  8:414862/149‑1 (MQ=255)

                                                                                                                                                                            cgGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATaaa                                                                                                                                                                                                                            <  8:251020/149‑1 (MQ=255)

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

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

                                                                                                                                                                                 tgatgaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCt                                                                                                                                                                                                                       <  8:195092/149‑1 (MQ=255)

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

                                                                                                                                                                                        ccGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCa                                                                                                                                                                                                                >  7:607741/1‑149 (MQ=255)

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

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

                                                                                                                                                                                               cAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAg                                                                                                                                                                                                         >  5:159634/1‑149 (MQ=255)

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

                                                                                                                                                                                                     ctctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATtat                                                                                                                                                                                                                               <  6:746504/121‑1 (MQ=255)

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

                                                                                                                                                                                                     ctctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGa                                                                                                                                                                                                   <  6:159634/149‑1 (MQ=255)

                                                                                                                                                                                                     ctctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGa                                                                                                                                                                                                   <  6:126909/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                                                                               gcgAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACAACTGGTGATggc                                                                                                                                                                                               <  5:399456/143‑3 (MQ=255)

                                                                                                                                                                                                                 gAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATg                                                                                                                                                                                       <  7:784890/149‑1 (MQ=255)

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

                                                                                                                                                                                                                    cAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCg                                                                                                                                                                                                                                         >  4:310258/1‑96 (MQ=255)

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

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

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

                                                                                                                                                                                                                      gggCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCtt                                                                                                                                                                                  >  7:958194/1‑149 (MQ=255)

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

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

                                                                                                                                                                                                                          tttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTc                                                                                                                                                                              >  7:178336/1‑149 (MQ=255)

                                                                                                                                                                                                                            tAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCAc                                                                                                                                                                            <  7:334198/149‑1 (MQ=255)

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

                                                                                                                                                                                                                              aCAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGc                                                                                                                                                                          >  7:675918/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:4887/149‑1 (MQ=255)

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

                                                                                                                                                                                                                                         gTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGGCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAg                                                                                                                                                               >  7:440780/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                           atatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGt                                                                                                                                                             <  8:334144/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)

                                                                                                                                                                                                                                              tCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCccaac                                                                                                                                                                     >  7:764696/1‑138 (MQ=255)

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

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

                                                                                                                                                                                                                                                gCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATa                                                                                                                                                        >  6:1324135/1‑149 (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)

                                                                                                                                                                                                                                                        tttGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAg                                                                                                                                                <  5:5351/149‑1 (MQ=255)

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

                                                                                                                                                                                                                                                              gTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCaaa                                                                                                                                          <  7:843265/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                                                                                                                                       tCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGc                                                                                                                                 >  7:368633/1‑149 (MQ=255)

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

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

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

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

TCAATGGCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCAGCGTCGGGAAGCCGGTCATATTACCAATCGGGAAG  >  NC_000913/442750‑443288

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