BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,057,150	C→T	5.5%	intergenic (+162/+28)	sibC → / ← serA	sRNA antisense regulator of toxic IbsC protein/D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,057,150	0	C	T	5.5%	761.8 / 21.4	231	intergenic (+162/+28)	sibC/serA	sRNA antisense regulator of toxic IbsC protein/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand): ref base C (97/120); new base T (6/7); total (103/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

Only 100 of 232 total aligned reads displayed.

ACTGACTGTTAATAAGCGCTGAAACTTATGAGTAACAGTACAATCAGTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTGGATGTGCATCAGACGACGCCCACCGTGCAGTGGCAGCGAGA  >  NC_000913/3056882‑3057422

                                                                         ttCTCCTTCAAGCCCCCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGCCTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGACTGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGCAGGGGAATTAGTACAGAAGACGGGCGCGAATGGTACCAGGAATAGCTTTGATTGCCAGc                                                                                                                                                                                                   >  2:135811/1‑275 (MQ=255)

                                                                                         cGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGCGAAGGGGAATTAGTACAGCAGCCGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCaa                                                                                                                                                                                   >  2:114550/1‑275 (MQ=255)

                                                                                                          tACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtca                                                                                                                                                                  >  2:50703/1‑275 (MQ=255)

                                                                                                                 aCAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCagag                                                                                                                                                                                                                                                               >  2:352724/1‑175 (MQ=255)

                                                                                                                     ccccATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                                                                                                                                                                                                           <  4:352724/175‑1 (MQ=255)

                                                                                                                         aTGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCt                                                                                                                                                                                                                                                                                 <  6:1229687/149‑1 (MQ=255)

                                                                                                                          tgctgcCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAAAGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATAGTGCCCGGAATCGCTGTCATTTCCTGAAGAGCTTTTTCGGACACGCCTTCGTCGGCGTCAATATCAATAACcaact                                                                                                                                                  >  4:25666/1‑272 (MQ=255)

                                                                                                                             tgccctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGccc                                                                                                                                                                                                                                                                             <  5:258368/149‑1 (MQ=255)

                                                                                                                             tgccctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                                                                                                                                                                                                                                                                             <  5:1412062/149‑1 (MQ=255)

                                                                                                                                cctgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGtt                                                                                                                                                                                                                                                                          <  7:747370/149‑1 (MQ=255)

                                                                                                                                  tgccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                                                                                                                                                                                                        >  7:939718/1‑149 (MQ=255)

                                                                                                                                    ccATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                                                                                                                                                                                                                      >  8:640187/1‑149 (MQ=255)

                                                                                                                                      aTCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtttt                                                                                                                                                                                                                                                                    >  5:604917/1‑149 (MQ=255)

                                                                                                                                       tCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                                                                                                                                                                                                   >  8:737739/1‑149 (MQ=255)

                                                                                                                                       tCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                                                                                                                                                                                                   <  8:1005693/149‑1 (MQ=255)

                                                                                                                                        cGTAAACCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCa                                                                                                                                                                                                                                                                  <  6:1477576/149‑1 (MQ=255)

                                                                                                                                         gTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCACTAACCACCTAACCCCTCTGGGCgg                                                                                                                                   >  2:238759/1‑275 (MQ=255)

                                                                                                                                          tAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCaga                                                                                                                                                                                                                                                                <  5:834494/149‑1 (MQ=255)

                                                                                                                                             aTCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGaa                                                                                                                                                                                                                                                             <  8:1092605/149‑1 (MQ=255)

                                                                                                                                             aTCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGaa                                                                                                                                                                                                                                                             <  5:1204431/149‑1 (MQ=255)

                                                                                                                                              tCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGAT‑AAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                                                                                                                                                                                                           <  5:1009579/149‑1 (MQ=255)

                                                                                                                                               ccccATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                                                                                                                                                                                                           <  6:1275187/149‑1 (MQ=255)

                                                                                                                                               ccccATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgg                                                                                                                                                                                                                                                           <  8:473408/149‑1 (MQ=255)

                                                                                                                                                 ccATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                                                                                                                                                                                                         >  6:813592/1‑149 (MQ=255)

                                                                                                                                                   aTTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGTGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                                                                                                                                                       >  5:80876/1‑149 (MQ=255)

                                                                                                                                                   aTTAAATAAACACAACGCATTGATCGGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                                                                                                                                                                                                       >  5:414035/1‑149 (MQ=255)

                                                                                                                                                    ttAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                                                                                                                                                                                                        >  2:39727/1‑131 (MQ=255)

                                                                                                                                                    ttAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAt                                                                                                                                                                                                                                                      <  7:440516/149‑1 (MQ=255)

                                                                                                                                                    ttAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAt                                                                                                                                                                                                                                                      <  8:110902/149‑1 (MQ=255)

                                                                                                                                                     taaataaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                                                                                                                                                                                                                     <  6:450333/149‑1 (MQ=255)

                                                                                                                                                        ataaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtttt                                                                                                                                                                                                                                                                    <  4:39727/131‑1 (MQ=255)

                                                                                                                                                         taaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTa                                                                                                                                                                                                                                                 <  5:631257/149‑1 (MQ=255)

                                                                                                                                                         taaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTa                                                                                                                                                                                                                                                 >  8:88521/1‑149 (MQ=255)

                                                                                                                                                         taaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTa                                                                                                                                                                                                                                                 <  6:604917/149‑1 (MQ=255)

                                                                                                                                                          aaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAc                                                                                                                                                                                                                                                >  8:274283/1‑149 (MQ=255)

                                                                                                                                                             cacaACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagc                                                                                                                                                                                                                                             <  5:57017/149‑1 (MQ=255)

                                                                                                                                                               caACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTCGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagcag                                                                                                                                                                                                                                           >  8:956195/1‑149 (MQ=255)

                                                                                                                                                                aaCGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTCGccc                                                                                                                                                                                                                                                                             <  8:1200567/114‑1 (MQ=255)

                                                                                                                                                                aaCGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTCGccc                                                                                                                                                                                                                                                                             >  7:1200567/1‑114 (MQ=255)

                                                                                                                                                                   gCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACgg                                                                                                                                                                                                                                       <  6:1063989/149‑1 (MQ=255)

                                                                                                                                                                   gCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACgg                                                                                                                                                                                                                                       <  7:1063792/149‑1 (MQ=255)

                                                                                                                                                                     aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgc                                                                                                                                                                                                                                     <  6:226645/149‑1 (MQ=255)

                                                                                                                                                                          tCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGaa                                                                                                                                                                                                                                <  5:1084386/149‑1 (MQ=255)

                                                                                                                                                                           cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAt                                                                                                                                                                                                                               <  8:1180305/149‑1 (MQ=255)

                                                                                                                                                                           cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAt                                                                                                                                                                                                                               <  7:640187/149‑1 (MQ=255)

                                                                                                                                                                              aCTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGAc                                                                                                                                                                                                                                         >  3:45336/1‑136 (MQ=255)

                                                                                                                                                                               cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTa                                                                                                                                                                                                                           <  7:1335128/149‑1 (MQ=255)

                                                                                                                                                                                tttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAc                                                                                                                                                                                                                          >  6:420559/1‑149 (MQ=255)

                                                                                                                                                                                  tGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgc                                                                                                                                                                                                                                     <  1:45336/136‑1 (MQ=255)

                                                                                                                                                                                   gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCgg                                                                                         >  2:126768/1‑275 (MQ=255)

                                                                                                                                                                                      ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGaa                                                                                                                                                                                                                    >  7:228686/1‑149 (MQ=255)

                                                                                                                                                                                      ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGaa                                                                                                                                                                                                                    >  7:103440/1‑149 (MQ=255)

                                                                                                                                                                                       tatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                                                                                                                                                                                                                   <  7:1099351/149‑1 (MQ=255)

                                                                                                                                                                                        atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATa                                                                                                                                                                                                                  <  7:88521/149‑1 (MQ=255)

                                                                                                                                                                                          tttCTGGAGGAGACTCGCAAAGTAGAATGCGCAGCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGc                                                                                                                                                                                                                <  6:664565/149‑1 (MQ=255)

                                                                                                                                                                                            tCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCtt                                                                                                                                                                                                              >  8:804262/1‑149 (MQ=255)

                                                                                                                                                                                             cTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCttt                                                                                                                                                                                                             <  6:346800/149‑1 (MQ=255)

                                                                                                                                                                                               ggAACAGACTCGCAAAGTAGAATGGAGAGTGTGGCAACCGTGTGGAGAAGGGATAAAAAAACGGGTAAGTCAGTGACCTGCCTGTTGATTGTCAGATAAGGCGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGt                                                                             <  4:126159/275‑1 (MQ=255)

                                                                                                                                                                                                gAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAt                                                                                                                                                                                                          >  6:736949/1‑149 (MQ=255)

                                                                                                                                                                                                 aGCAGACTCGCAAAGTAGAAGGCGAAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGTGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAtt                                                                                                                                                                                                         <  8:401999/149‑1 (MQ=255)

                                                                                                                                                                                                    aGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc                                                                                                                                                                                                      <  8:1242521/149‑1 (MQ=255)

                                                                                                                                                                                                     gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAACACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt                                                                                                                                                                                                     <  8:1429872/149‑1 (MQ=255)

                                                                                                                                                                                                     gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt                                                                                                                                                                                                     <  6:714946/149‑1 (MQ=255)

                                                                                                                                                                                                        tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                                                                                                                                                                                                  >  5:854602/1‑149 (MQ=255)

                                                                                                                                                                                                        tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                                                                                                                                                                                                  >  5:9566/1‑149 (MQ=255)

                                                                                                                                                                                                          gCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc                                                                                                                                                                                                >  5:646383/1‑149 (MQ=255)

                                                                                                                                                                                                           cAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGa                                                                                      <  4:130345/254‑1 (MQ=255)

                                                                                                                                                                                                                  gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCg                                                                                                                                                                                                                       <  8:550749/118‑1 (MQ=255)

                                                                                                                                                                                                                    aacgacAACGCGGCAACGGGTTCGGGAAGCGATAAAACAACTGGAAAGCCAGTGACCTGTCCGTTGGTTTTCAGTGAAGGGGAATCATTACAGCAGACGGGCTCGAATGGTACCCGGCATGGCTTTCAGTGCCTGCAGCGCTTTTTCGGCAACGTCTGCTTCTGCTTCAATATCAATAACCACATAGCCCATCTGGGCGGAAGTTTGCAGAGATTGCGCGGCGATGTTGACTCCCTGCTCTGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCg                                                        <  4:21182/270‑1 (MQ=255)

                                                                                                                                                                                                                     tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTc                                                                                                                                                                                        <  6:547827/146‑1 (MQ=255)

                                                                                                                                                                                                                        gcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAAc                                                                                                                                                                                  <  8:1270000/149‑1 (MQ=255)

                                                                                                                                                                                                                           aCGCGGCAACGGTGTGGAGAAGAGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTGGGCAACGTc                                                                                                                                                                               <  5:1270568/149‑1 (MQ=255)

                                                                                                                                                                                                                              cgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtc                                                                                                                                                                            <  6:646383/149‑1 (MQ=255)

                                                                                                                                                                                                                              cgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtc                                                                                                                                                                            <  6:747803/149‑1 (MQ=255)

                                                                                                                                                                                                                               ggCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCg                                             >  4:248073/1‑275 (MQ=255)

                                                                                                                                                                                                                                 cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAAACAAATAACCCATCTGGGCGGAAGGTTGCAGATATTGCGCGGCGATGTTGCCGCACTGCTCGGCGAAGCTTGTGTTCAGCGCCGTTAGAACTCACGGACGGTGTTCGTg                                           >  2:286045/1‑275 (MQ=255)

                                                                                                                                                                                                                                  aaCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgtcg                                                                                                                                                                        <  6:1084623/149‑1 (MQ=255)

                                                                                                                                                                                                                                     ggTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCt                                                                                                                                                                     <  5:736949/149‑1 (MQ=255)

                                                                                                                                                                                                                                     ggTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCt                                                                                                                                                                     <  7:1292929/149‑1 (MQ=255)

                                                                                                                                                                                                                                             gaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatc                                                                                                                                                             <  6:883682/149‑1 (MQ=255)

                                                                                                                                                                                                                                               aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaa                                                                                                                                                           <  5:431247/149‑1 (MQ=255)

                                                                                                                                                                                                                                               aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaa                                                                                                                                                           <  8:284037/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 ggATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaata                                                                                                                                                         >  6:1149629/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 ggATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaata                                                                                                                                                         >  8:1388696/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 ggATAAAAAAACGGGCAAGGCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGAAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaata                                                                                                                                                         >  8:920590/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  gATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGGTTCAATATCAATaa                                                                                                                                                        >  7:54280/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aTAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAAc                                                                                                                                                       >  5:233245/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    tAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTa                                                                                                                                                                                                                           <  8:697876/80‑1 (MQ=255)

                                                                                                                                                                                                                                                    tAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTa                                                                                                                                                                                                                           >  7:697876/1‑80 (MQ=255)

                                                                                                                                                                                                                                                     aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                                                                                                                                                                                                  <  6:1026831/104‑1 (MQ=255)

                                                                                                                                                                                                                                                     aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACca                                                                                                                                                     >  6:1340039/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      aatcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACcac                                                                                                                                                    >  6:873480/5‑149 (MQ=255)

                                                                                                                                                                                                                                                       aaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACcaca                                                                                                                                                   >  5:927992/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              ggCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCa                                                                                                                                            >  6:193376/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               gCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGCACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCAt                                                                                                                                           >  5:1105360/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                cAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATc                                                                                                                                          >  7:825543/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  aGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTg                                                                                                                                        >  7:922167/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                          aCCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAg                                                                                                                                >  5:579601/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                          aCCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTGGATGTGCATCAGACGACGCCCACCGTGCAGTGGCAGCgaga  <  2:248073/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                            cTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGtt                                                                                                                              <  5:1340039/149‑1 (MQ=255)

ACTGACTGTTAATAAGCGCTGAAACTTATGAGTAACAGTACAATCAGTATGATGACAAGTCGCATCATAACCCTTCTCCTTCAAGCCCTCGCTTCGGTGAGGGCTTTACCGTTACAGCCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACATAACCCATCTGGGCGGAAGTTTGCAGATATTGCGCGGCGATGTTGACGCCCTGCTCGGCGAAGATTTTGTTCAGCGCAGTTAGCACGCCCGGACGGTTTTCGTGGATGTGCATCAGACGACGCCCACCGTGCAGTGGCAGCGAGA  >  NC_000913/3056882‑3057422

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