BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,425,975	A→G	5.0%	intergenic (‑59/+31)	hisQ ← / ← hisJ	histidine ABC transporter permease/histidine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,425,975	0	A	G	5.0%	696.2 / 16.8	210	intergenic (‑59/+31)	hisQ/hisJ	histidine ABC transporter permease/histidine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand): ref base A (110/88); new base G (8/3); total (118/92)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-03

Only 100 of 215 total aligned reads displayed.

CCGGCACGCCACGAATCAGCGTGGTGTACCCTTCGAAAATCAGCCCCGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACC  >  NC_000913/2425736‑2426237

                                           ccccGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTCGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTCGCCACCATAAACATCAAAATCGAAGTACCTTTTCGCTACTTtct                                                                                                                                                                                           >  1:482/1‑274 (MQ=255)

                                                                                             tgatTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAg                                                                                                                                                                                                                                                                      <  8:961659/149‑1 (MQ=255)

                                                                                                  aCAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGc                                                                                                                                                                                                                                                                 <  8:776769/149‑1 (MQ=255)

                                                                                                     gCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCggg                                                                                                                                                                                                                                                              <  6:1044339/149‑1 (MQ=255)

                                                                                                           aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                                                                                                                                                        <  6:413471/149‑1 (MQ=255)

                                                                                                           aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                                                                                                                                                        <  5:17549/149‑1 (MQ=255)

                                                                                                           aCTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttct                                                                                                                                                                                                                                                        <  7:511808/149‑1 (MQ=255)

                                                                                                            cTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctttctt                                                                                                                                                                                                                                                       >  6:1480910/1‑149 (MQ=255)

                                                                                                                agagCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCAc                                                                                                                                                                                                                                                   >  5:797078/1‑149 (MQ=255)

                                                                                                                   gCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGa                                                                                                                                                                                                                                                >  8:344773/1‑149 (MQ=255)

                                                                                                                     tGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                                                                                                                                                              <  5:1382921/149‑1 (MQ=255)

                                                                                                                     tGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                                                                                                                                                              <  8:1147275/149‑1 (MQ=255)

                                                                                                                       aTAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAACTCGAAGTAATTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCCCATTTCGTCACAGGCTTTGTTCAGTGCTTCGCACCGTTCGTTCTCTTCTTTACGCAg                                                                                                               >  4:81960/1‑274 (MQ=255)

                                                                                                                       aTAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTNCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAANCCGAATTTCGGCAACGGCTTTGTTCAGTGCNNTTCGCAGTTCGTTATCTTCTTTACTCa                                                                                                                >  3:48309/1‑273 (MQ=255)

                                                                                                                       aTAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGATTTCTTTCACCGATTGCCACTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTCCCGTCAGCGCGCATTTCGGCAACGGCTTTGTTCATTGCTTCGAGCAGTTCTTTATCTTCTTTACGCAg                                                                                                               >  2:283109/1‑274 (MQ=255)

                                                                                                                         aGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGa                                                                                                                                                                                                                                          <  7:763397/149‑1 (MQ=255)

                                                                                                                           ccAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGAtt                                                                                                                                                                                                                                        >  8:497947/1‑149 (MQ=255)

                                                                                                                            cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                                                                                                                                                       >  7:554925/1‑149 (MQ=255)

                                                                                                                            cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                                                                                                                                                       >  5:1230335/1‑149 (MQ=255)

                                                                                                                            cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTa                                                                                                                                                                                                                                       <  6:304009/149‑1 (MQ=255)

                                                                                                                                tCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGcca                                                                                                                                                                                                                                   >  5:758241/1‑149 (MQ=255)

                                                                                                                                   aGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGccacca                                                                                                                                                                                                                                <  5:836074/149‑1 (MQ=255)

                                                                                                                                         aCGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGCATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGTAGGGGGCTTTCTTTCACCAATGGCGGTTAGCCACCATAAACATCAAAATCGAAGTACTCTTTCGCTAATTTCCCGTAAGTACCGTAGACGCGCATTTCGGCACAGGCTTTGTTCAGTCCTTCGCGCCGTTCGTTATACTCTTTACCGAGTCCCCTGCCGCTCCCTAc                                                                                             >  4:222341/1‑274 (MQ=255)

                                                                                                                                         aCGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACa                                                                                                                                                                                                                          >  5:597279/1‑149 (MQ=255)

                                                                                                                                           gAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATc                                                                                                                                                                                                                        <  5:240916/149‑1 (MQ=255)

                                                                                                                                            aGCGCACCCTCTAAAATAAACCCTGAAAACCCTTACAACTGGCAGCCTGGCCTGT‑GAGTGTTGGTGAGTTCCCGCCTCTTGTCTCCCTCCGGGTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATACACATCCAAATCGAAGTACTTTGTCGCTAATTTCTCGTAAGTACCGTCAGCGCACATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCa                                                                                         <  4:207897/275‑1 (MQ=255)

                                                                                                                                                   ccTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGa                                                                                                                                                                                                                >  6:750083/1‑149 (MQ=255)

                                                                                                                                                    cTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                               <  6:77981/149‑1 (MQ=255)

                                                                                                                                                     tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctt                                                                                                                                                                                                                                                           >  7:959139/1‑104 (MQ=255)

                                                                                                                                                     tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGctt                                                                                                                                                                                                                                                           <  8:959139/104‑1 (MQ=255)

                                                                                                                                                     tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAg                                                                                                                                                                                                              >  5:586489/1‑149 (MQ=255)

                                                                                                                                                       tAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGACGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCANAATCGAAGTACTTTTTCGCTCATTTCCCGTAAGTACCTTCAGCGCGCATTTCGGCAAAGGCTTTTTCTAGTGCTNCCCCCATTTCTTTTTTTTCTTTACGCAGGCNNATGCCCGTCACTACCCCAACCAGTTTTTc                                                                              >  3:47882/1‑275 (MQ=255)

                                                                                                                                                        aaaaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTAc                                                                                                                                                                                                           <  8:144473/149‑1 (MQ=255)

                                                                                                                                                         aaaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTAc                                                                                                                                                                                                           >  7:1374288/1‑148 (MQ=255)

                                                                                                                                                          aaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACtt                                                                                                                                                                                                         <  7:534173/149‑1 (MQ=255)

                                                                                                                                                            tAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACtttt                                                                                                                                                                                                       >  7:612986/1‑149 (MQ=255)

                                                                                                                                                                acCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGc                                                                                                                                                                                                   <  6:1081809/149‑1 (MQ=255)

                                                                                                                                                                 ccTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGAtt                                                                                                                                                                                                                                        <  6:662315/111‑1 (MQ=255)

                                                                                                                                                                  cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAggg                                                                                                                                                                                                                                                                    <  6:262356/82‑1 (MQ=255)

                                                                                                                                                                  cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                               >  5:1293783/1‑135 (MQ=255)

                                                                                                                                                                  cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                               <  6:1293783/135‑1 (MQ=255)

                                                                                                                                                                   tGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTaa                                                                                                                                                                                                >  7:977406/1‑149 (MQ=255)

                                                                                                                                                                        aCCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                                                                                                                                                              >  8:58211/1‑98 (MQ=255)

                                                                                                                                                                        aCCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATg                                                                                                                                                                                                                                              <  7:58211/98‑1 (MQ=255)

                                                                                                                                                                        aCCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtct                                                                                                                                                                                           >  8:1213635/1‑149 (MQ=255)

                                                                                                                                                                          ccATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGAt                                                                                                                                                                                                                                               <  6:933062/95‑1 (MQ=255)

                                                                                                                                                                          ccATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGAt                                                                                                                                                                                                                                               >  5:933062/1‑95 (MQ=255)

                                                                                                                                                                          ccATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAAc                                                                                                                                                                                                                           <  6:850506/115‑1 (MQ=255)

                                                                                                                                                                          ccATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                                                                                                                                                         >  8:1004632/1‑149 (MQ=255)

                                                                                                                                                                           cATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                               >  5:472356/1‑126 (MQ=255)

                                                                                                                                                                           cATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                                                                                                                                               <  6:472356/126‑1 (MQ=255)

                                                                                                                                                                             tACAACCTGCTGCCTGTCCTGTCGTGTG‑TGGTAAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCg                                                        <  1:34126/275‑1 (MQ=255)

                                                                                                                                                                             tACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTaa                                                                                                                                                                                      >  7:473804/1‑149 (MQ=255)

                                                                                                                                                                                       tgcCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCg                                                                                                     >  4:237075/1‑220 (MQ=255)

                                                                                                                                                                                        gcCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTTTCCCTCCGGTTACCCGAAGGTCGCGGGCTTTCTTTCACCGATGGTGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgc                                                                                                                                                                           >  7:1398930/1‑149 (MQ=255)

                                                                                                                                                                                        gcCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgc                                                                                                                                                                           <  8:554925/149‑1 (MQ=255)

                                                                                                                                                                                        gcCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTGTCATCTTTACCGAACGGGCCCCCGGATTTggactc                                             >  2:305253/1‑270 (MQ=255)

                                                                                                                                                                                        gcCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATC‑AAATCGaa                                                                                                                                                                                                               >  7:1154481/1‑112 (MQ=255)

                                                                                                                                                                                         cctgtcctgtCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTcg                                                                                                                                         >  2:71326/1‑182 (MQ=255)

                                                                                                                                                                                            gtcctgtcGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgcgc                                                                                                                                                                       <  8:1357297/149‑1 (MQ=255)

                                                                                                                                                                                             tcctgtcGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCa                                                                                                                                                                      >  8:70163/1‑149 (MQ=255)

                                                                                                                                                                                             tcctgtcGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCAGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCa                                                                                                                                                                      >  5:1165006/1‑149 (MQ=255)

                                                                                                                                                                                               ctgtcGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCAtt                                                                                                                                                                    >  5:90065/1‑149 (MQ=255)

                                                                                                                                                                                                 gtcGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTc                                                                                                                                                                  <  8:281896/149‑1 (MQ=255)

                                                                                                                                                                                                 gtcGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTc                                                                                                                                                                  <  7:210911/149‑1 (MQ=255)

                                                                                                                                                                                                   cGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCgg                                                                                                                                                                <  5:992438/149‑1 (MQ=255)

                                                                                                                                                                                                      gtg‑tgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCaa                                                                                                                                                             >  5:261981/1‑149 (MQ=255)

                                                                                                                                                                                                       tg‑tgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCaaa                                                                                                                                                            <  6:597279/149‑1 (MQ=255)

                                                                                                                                                                                                       tg‑tgGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCAAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCCAAAGgtt                               >  4:302405/1‑274 (MQ=255)

                                                                                                                                                                                                          tggtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAgg                                                                                                                                                          >  5:778538/1‑149 (MQ=255)

                                                                                                                                                                                                            gtgAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAg                                                                                                                                                           >  7:1139083/1‑146 (MQ=255)

                                                                                                                                                                                                              gAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCttt                                                                                                                                                      <  5:662694/149‑1 (MQ=255)

                                                                                                                                                                                                               aaTTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCgg                                                                                                                                                                >  6:647611/1‑138 (MQ=255)

                                                                                                                                                                                                               aaTTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCgg                                                                                                                                                                <  5:647611/138‑1 (MQ=255)

                                                                                                                                                                                                                  tACCGCCTCTTGGCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACAATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGGACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTGAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGaaa                    <  2:219041/275‑1 (MQ=255)

                                                                                                                                                                                                                    ccGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAg                                                                                                                                                <  8:670435/149‑1 (MQ=255)

                                                                                                                                                                                                                    ccGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAg                                                                                                                                                <  6:1357732/149‑1 (MQ=255)

                                                                                                                                                                                                                     cGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                                                                                                                                                          >  7:1049388/1‑106 (MQ=255)

                                                                                                                                                                                                                     cGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGt                                                                                                                                               >  6:1230534/1‑149 (MQ=255)

                                                                                                                                                                                                                        ctctTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCt                                                                                                                                            >  7:1010720/1‑149 (MQ=255)

                                                                                                                                                                                                                          ctTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTc                                                                                                                                          >  6:812818/1‑149 (MQ=255)

                                                                                                                                                                                                                            ccgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                                                                                                                                                         <  7:567796/97‑1 (MQ=38)

                                                                                                                                                                                                                             cgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAAtt                                                                                                                                                                                              >  5:1326271/3‑94 (MQ=38)

                                                                                                                                                                                                                             cgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAAtt                                                                                                                                                                                              <  6:1326271/92‑1 (MQ=38)

                                                                                                                                                                                                                             cgctcccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTcgcg                                                                                                                                       >  7:1203352/3‑149 (MQ=255)

                                                                                                                                                                                                                                tccctccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCtt                                                                                                                                           >  8:109516/1‑142 (MQ=255)

                                                                                                                                                                                                                                  ccttcGGGTAACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                                                                                                                                                          >  5:271023/1‑93 (MQ=255)

                                                                                                                                                                                                                                    tccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGAcc  <  2:109772/275‑1 (MQ=255)

                                                                                                                                                                                                                                    tccGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGAcc  <  4:62736/275‑1 (MQ=255)

                                                                                                                                                                                                                                       ggTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTa                                                                                                                             <  7:1433257/149‑1 (MQ=255)

                                                                                                                                                                                                                                        gTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCg                                                                                                                                >  4:138145/1‑145 (MQ=255)

                                                                                                                                                                                                                                           aCCCGAAGGGAGTGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtctt                                                                                                                         <  6:580450/149‑1 (MQ=255)

                                                                                                                                                                                                                                           aCCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtctt                                                                                                                         >  6:671792/1‑149 (MQ=255)

                                                                                                                                                                                                                                            cccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtcttc                                                                                                                        >  7:335371/1‑149 (MQ=255)

                                                                                                                                                                                                                                            cccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAt                                                                                                                            <  2:138145/145‑1 (MQ=255)

                                                                                                                                                                                                                                             ccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATGGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtcttct                                                                                                                       >  6:1209271/1‑149 (MQ=255)

                                                                                                                                                                                                                                             ccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCTCGCATTTCGGCGAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtcttct                                                                                                                       >  8:233518/1‑149 (MQ=255)

                                                                                                                                                                                                                                             ccGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTAtcttc                                                                                                                        <  5:658215/148‑1 (MQ=255)

                                                                                                                                                                                                                                              cGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATcttctt                                                                                                                      >  6:1127390/1‑149 (MQ=255)

                                                                                                                                                                                                                                              cGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATcttctt                                                                                                                      <  7:285676/149‑1 (MQ=255)

CCGGCACGCCACGAATCAGCGTGGTGTACCCTTCGAAAATCAGCCCCGAAAGCCGATTTTGCGAGAGCTTACCGCCAGCGCCAATTAAACCGATGATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTG‑TGGTGAATTACCGCCTCTTGTCTCCCTCCGGTTACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACC  >  NC_000913/2425736‑2426237

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