BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,253,830	C→G	6.3%	S708T (AGT→ACT)	ycgV ←	putative adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,253,830	0	C	G	6.3%	1083.0 / 49.6	303	S708T (AGT→ACT)	ycgV	putative adhesin
Reads supporting (aligned to +/- strand): ref base C (148/135); new base G (11/8); total (159/144)
Fisher's exact test for biased strand distribution p-value = 8.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-04

Only 100 of 304 total aligned reads displayed.

ACTGTCCAGTACGTGGAAACTATTTTTCTGGCGCGATGCTTTTATAACGAGATCGCTGTAAAAACCGTTTTGTGCCATGTAACTGGCGTACATTCCCATGTAGTCTGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGATCGGGGTTGACTATGGGAGGCTGAGCGGGAGCTGGTGTGGGTTCAGGATTAGGAGTCGGTTCCGGAAC  >  NC_000913/1253571‑1254095

acgttCCATTACGCGGAAACTATTTTTCTGGCGCAATTCGTTTATAACGTGATCGCTGTAAAAACCGTTTTGTGCCATGTAACTGGCGTACATTCCCATGTAGTCTGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAATCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCt                                                                                                                                                                                                                                                        <  1:32156/275‑1 (MQ=255)

  tGTCCAGTACGTGGAAACTATTTTTCTGGCGCGATGCTTTTATAACGAGATCGCTGTAAAAACCGTTTTGTGCCATGTAACTGGCGTACATTCCCATGTAGTCTGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAgg                                                                                                                                                                                                                                                          <  2:161403/275‑1 (MQ=255)

                           cTGGCGCGATGCTTTTATAACCAGATCGCTGTAAAAACCGTGTTATGCCATGTAACTGGCGTACATTCCCATGTAGTCTGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCTGCGATGCATGTGTTTAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGa                                                                                                                                                                                                                                 <  2:701627/275‑1 (MQ=255)

                                                                  gTTTTGTGCCATGTAACTGGCGTACATTCCCATGTAGTCTGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCCTGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTc                                                                                                                                                                                                                                                              >  1:291740/1‑207 (MQ=255)

                                                                                  cttgCGTACATTCCCATGTAGTCTGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                                                                                                                                          <  2:611693/272‑1 (MQ=255)

                                                                                               ccATGTAGTCTGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCATATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                                                                                                                                              >  2:67286/1‑274 (MQ=255)

                                                                                               ccATGTAGTCTGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAg                                                                                                                                                             >  2:439449/1‑275 (MQ=255)

                                                                                                         tGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCa                                                                                                                                                                                                                                                            >  4:409109/1‑170 (MQ=255)

                                                                                                         tGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGCTTGCCACTGGCAAAGGAGTCCAGGCt                                                                                                                                                                                                                                                        >  1:190607/1‑174 (MQ=255)

                                                                                                             cGTGCGGTACCGTCGCCTCCGCTATAGTCCGACGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAcc                                                                                                                                               <  2:293212/275‑1 (MQ=255)

                                                                                                                      ccGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                                                                                                                                                    >  8:1269642/1‑149 (MQ=255)

                                                                                                                      ccGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                                                                                                                                                    <  8:951211/149‑1 (MQ=255)

                                                                                                                          cGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTTACCCATGTCAAAGCCGCTCAGTTTGCCAGTGGCaaactga                                                                                                                                                                                                                                                                <  5:499097/149‑5 (MQ=255)

                                                                                                                              tCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCa                                                                                                                                                                                                                                                            <  7:791202/149‑1 (MQ=255)

                                                                                                                              tCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCa                                                                                                                                                                                                                                                            >  5:917811/1‑149 (MQ=255)

                                                                                                                                   tataGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACAACAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCtt                                                                                                                                                                                                                                                       >  5:1134108/1‑149 (MQ=255)

                                                                                                                                        tCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTccccc                                                                                                                                                                                                                                                  >  8:107933/1‑149 (MQ=255)

                                                                                                                                         ccGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGGGTCCAGGCTTCCCCCa                                                                                                                                                                                                                                                 >  7:241393/1‑149 (MQ=255)

                                                                                                                                         ccGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCa                                                                                                                                                                                                                                                 >  5:185643/1‑149 (MQ=255)

                                                                                                                                                tGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                                                                                                                                          <  8:1007047/149‑1 (MQ=255)

                                                                                                                                                tGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                                                                                                                                          <  6:1586067/149‑1 (MQ=255)

                                                                                                                                                tGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                                                                                                                                          <  6:1651338/149‑1 (MQ=255)

                                                                                                                                                 gCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgcg                                                                                                                                                                                                                                         >  7:1008795/1‑149 (MQ=255)

                                                                                                                                                    tgtgTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCaa                                                                                                                                                                                                                                      >  8:887560/1‑149 (MQ=255)

                                                                                                                                                    tgtgTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCaa                                                                                                                                                                                                                                      >  7:1531801/1‑149 (MQ=255)

                                                                                                                                                       gtTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCa                                                                                                                                                                                                                                   >  5:626654/1‑149 (MQ=255)

                                                                                                                                                        ttGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTACAGGCTTCCCCCATAACTgc                                                                                                                                                                                                                                          >  7:405757/1‑141 (MQ=255)

                                                                                                                                                           aGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGata                                                                                                                                                                                                                               <  7:953775/149‑1 (MQ=255)

                                                                                                                                                           aGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTAc                                                                                                                                                                           >  4:747382/1‑201 (MQ=255)

                                                                                                                                                            gCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGatat                                                                                                                                                                                                                              >  5:700748/1‑149 (MQ=255)

                                                                                                                                                              cAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCa                                                                                                                                                                                                                                       >  8:1332331/1‑138 (MQ=255)

                                                                                                                                                               aaTATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTAc                                                                                                                                                                                                                           <  8:477458/149‑1 (MQ=255)

                                                                                                                                                                atatACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAGTGGCaaactg                                                                                                                                                                                                                                                                 <  8:986503/110‑4 (MQ=255)

                                                                                                                                                                atatACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAGTGGCaaactg                                                                                                                                                                                                                                                                 >  7:986503/1‑107 (MQ=255)

                                                                                                                                                                    aCAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCt                                                                                                                                                                                                                      >  5:1583879/1‑149 (MQ=255)

                                                                                                                                                                      aGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCttt                                                                                                                                                                                                                    >  6:31887/1‑149 (MQ=255)

                                                                                                                                                                       gACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTa                                                                                                                                                                                                                   <  6:131692/149‑1 (MQ=255)

                                                                                                                                                                       gACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCg                                                                                                                                                                              >  2:749930/1‑186 (MQ=255)

                                                                                                                                                                       gACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGGTTAGACGTGGgcgct                                                                                     >  2:311013/1‑274 (MQ=255)

                                                                                                                                                                           gACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCt                                                                                                                                                                                                                                                        <  4:62452/108‑1 (MQ=255)

                                                                                                                                                                            aCATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                                                                                                                                          >  8:1350786/1‑121 (MQ=255)

                                                                                                                                                                            aCATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                                                                                                                                          <  7:1350786/121‑1 (MQ=255)

                                                                                                                                                                             cATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTgg                                                                                                                                                                                                                                                                        >  2:144221/1‑90 (MQ=255)

                                                                                                                                                                             cATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGa                                                                                                                                                                                                             <  7:887560/149‑1 (MQ=255)

                                                                                                                                                                             cATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACa                                                                                                                                                                <  3:295329/194‑1 (MQ=255)

                                                                                                                                                                               tACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAg                                                                                                                                                                                                     <  4:150829/155‑1 (MQ=255)

                                                                                                                                                                                aCAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAttt                                                                                                                                                                                                          >  7:32630/1‑149 (MQ=255)

                                                                                                                                                                                   gnggCATTAAATCAGAGAGACGTTTCTCCCCACNNAACTGGATACCGCTGTAACCCATGTCAAAGNCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCNGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGAt                                                                     <  3:113458/277‑1 (MQ=255)

                                                                                                                                                                                    cGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCa                                                                                                                                                                                                      <  7:218173/149‑1 (MQ=255)

                                                                                                                                                                                     ggCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                                                                                                                                                    >  7:1342008/1‑86 (MQ=255)

                                                                                                                                                                                     ggCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAg                                                                                                                                                                                                     >  5:737953/1‑149 (MQ=255)

                                                                                                                                                                                           aCATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACAt                                                                                                                                                               <  4:305311/181‑1 (MQ=255)

                                                                                                                                                                                            cATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCAccc                                                                                                                                                                                              <  8:328064/149‑1 (MQ=255)

                                                                                                                                                                                             aTCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCa                                                                                                                                                                                             >  7:1156472/1‑149 (MQ=255)

                                                                                                                                                                                               cAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGAt                                                                     >  1:113458/1‑267 (MQ=255)

                                                                                                                                                                                               cAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATGGAGATACTTGCCGCCAGCATCTGCGGTCGTGGTTGGTTTAGGCGTGGCGTCGGTGCCAGCCGTAGGACGGGGGTTg                                                             >  4:476005/1‑275 (MQ=255)

                                                                                                                                                                                                  agagaCGTTTATGCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGt                                                                                                                                                                                        <  5:1358305/149‑1 (MQ=255)

                                                                                                                                                                                                  agagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGt                                                                                                                                                                                        >  6:659578/1‑149 (MQ=255)

                                                                                                                                                                                                   gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTAcc                                                                                                                                                                                                                          >  5:1455200/1‑114 (MQ=255)

                                                                                                                                                                                                   gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACGCATCCGtt                                                                                                                                                                                       <  7:1553348/149‑1 (MQ=255)

                                                                                                                                                                                                      aCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAATGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCa                                                                                                                                                                                    >  5:1209349/1‑149 (MQ=255)

                                                                                                                                                                                                      aCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCa                                                                                                                                                                                    >  8:605866/1‑149 (MQ=255)

                                                                                                                                                                                                       cGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCTCTGCCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCAt                                                                                                                                                                                   <  6:1227474/149‑1 (MQ=255)

                                                                                                                                                                                                            aTCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTa                                                                                                                                                                                                                            >  4:116268/1‑103 (MQ=255)

                                                                                                                                                                                                                ccaccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                                                                                                                                          <  5:1673552/149‑1 (MQ=255)

                                                                                                                                                                                                                  accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGt                                                                                                                                                                        >  7:1263132/1‑149 (MQ=255)

                                                                                                                                                                                                                  accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATAGTACCGTCTGTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGt                                                                                                                                                                        >  7:1454584/1‑149 (MQ=255)

                                                                                                                                                                                                                   ccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGtt                                                                                                                                                                       >  5:77816/1‑149 (MQ=255)

                                                                                                                                                                                                                     aaaCTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGtttt                                                                                                                                                                     >  6:103988/1‑149 (MQ=255)

                                                                                                                                                                                                                         tGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAAc                                                                                                                                                                 >  6:1491402/1‑149 (MQ=255)

                                                                                                                                                                                                                           gATACCGCTGTAACCAATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGc                                                                                                                                                                                     >  5:1255270/1‑127 (MQ=255)

                                                                                                                                                                                                                            aTACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                                                                                                                                              <  7:1378752/149‑1 (MQ=255)

                                                                                                                                                                                                                             tACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAg                                                                                                                                                             <  7:107933/149‑1 (MQ=255)

                                                                                                                                                                                                                                cGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCtt                                                                                                                                                                                                                                                       <  3:243949/56‑1 (MQ=255)

                                                                                                                                                                                                                                cGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCtt                                                                                                                                                                                                                                                       >  1:243949/1‑56 (MQ=255)

                                                                                                                                                                                                                                cGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTc                                                                                                                                                          <  7:761351/149‑1 (MQ=255)

                                                                                                                                                                                                                                   tGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAAt                                                                                                                                                       <  6:74440/149‑1 (MQ=255)

                                                                                                                                                                                                                                   tGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAAt                                                                                                                                                       <  5:461169/149‑1 (MQ=255)

                                                                                                                                                                                                                                    gTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                                                                                                                                              >  8:1181298/1‑141 (MQ=255)

                                                                                                                                                                                                                                       aCCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGt                                                                                                                                                   <  8:789338/149‑1 (MQ=255)

                                                                                                                                                                                                                                         ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTaa                                                                                                                                                 >  7:965767/1‑149 (MQ=255)

                                                                                                                                                                                                                                           aTGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAcc                                                                                                                                               <  8:169972/149‑1 (MQ=255)

                                                                                                                                                                                                                                           aTGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAcc                                                                                                                                               >  7:1309086/1‑149 (MQ=255)

                                                                                                                                                                                                                                              tCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGAc                                                                                                                                            >  7:521689/1‑149 (MQ=255)

                                                                                                                                                                                                                                               cacaGCCGCTCAGTTTGCCACTGGCTAAGGAGTCCAGGCTTCCCCCATAACTGCGCAGCCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGATCGGGGTTGACTATGGGAGGCTGAGCGGGAGCTGGTGTGGGTTCAGGATTAGGAGTc             <  2:476009/272‑1 (MQ=255)

                                                                                                                                                                                                                                                aaaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAt                                                                                                                                                                                                            <  6:1417267/83‑1 (MQ=255)

                                                                                                                                                                                                                                                aaaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCA‑CCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAc                                                                                                                                                >  8:1545273/1‑142 (MQ=255)

                                                                                                                                                                                                                                                 aaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACAtt                                                                                                                                         <  7:1620644/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  aGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAc                                                                                                                                                <  8:1365753/141‑1 (MQ=255)

                                                                                                                                                                                                                                                   gCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTga                                                                                                                                       <  8:1602650/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      gatCAGTGTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGATCGGGGTTGACTATGGGAGGCTGAGCGGGAGCTGGTGTGGGTTCAGGATTAGGAGTCGGTTCCGGAAc  <  1:84266/277‑1 (MQ=255)

                                                                                                                                                                                                                                                       cTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATa                                                                                                                                   >  6:878582/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        tCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATaa                                                                                                                                  <  5:276723/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        ctcttttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCTTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACTGACATTGAGATaa                                                                                                                                  >  8:1437607/5‑149 (MQ=255)

                                                                                                                                                                                                                                                         tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAg                                                                                                                                                                                                     >  6:198686/4‑81 (MQ=255)

                                                                                                                                                                                                                                                         tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGc                                                                                                                                                                                     >  6:495238/4‑97 (MQ=255)

                                                                                                                                                                                                                                                         tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAg                                                                                                                                                    >  6:408492/4‑130 (MQ=255)

                                                                                                                                                                                                                                                         tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAt                                                                                                                                 >  6:941816/4‑149 (MQ=255)

                                                                                                                                                                                                                                                           gTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTa                                                                                                                                                                            <  4:779281/104‑1 (MQ=255)

                                                                                                                                                                                                                                                           gTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGATCgggg                                                                <  2:569646/212‑1 (MQ=255)

ACTGTCCAGTACGTGGAAACTATTTTTCTGGCGCGATGCTTTTATAACGAGATCGCTGTAAAAACCGTTTTGTGCCATGTAACTGGCGTACATTCCCATGTAGTCTGAACGTGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGATCGGGGTTGACTATGGGAGGCTGAGCGGGAGCTGGTGTGGGTTCAGGATTAGGAGTCGGTTCCGGAAC  >  NC_000913/1253571‑1254095

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