BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,286,311	C→G	5.9%	intergenic (+24/+79)	proL → / ← yejO	tRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,286,311	0	C	G	5.9%	769.2 / 26.6	222	intergenic (+24/+79)	proL/yejO	tRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand): ref base C (102/105); new base G (5/8); total (109/113)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 228 total aligned reads displayed.

CGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCAC  >  NC_000913/2286050‑2286572

cGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCAAAGAAAAAACCAACCCTTACGGTTGGttt                                                                                                                                                                                                                                                           >  2:609495/1‑275 (MQ=255)

                     gAAAAACCACAGTTACGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACAAACCTTTACGGTTGGTTTTTATATTTcgtcatttattct                                                                                                                                                                                                                                      >  2:291222/1‑263 (MQ=255)

                             gcaGTTAATTTTTTTATGTAAAGTGCTGAAAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGCAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                                                                                                                                                                                              <  4:483735/274‑1 (MQ=255)

                                                           aGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCCACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTAtttt                                                                                                                                                                                               <  2:190534/275‑1 (MQ=255)

                                                              cGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTACCGTGCCGCCCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATATGCACTTAATTCGATCAACA‑ACCCTGTACCCATCACATCCTGTTTAtttttt                                                                                                                                                                                             >  4:190125/1‑272 (MQ=255)

                                                                   aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑acac                                                                                                                                                                                                                  >  2:373665/1‑248 (MQ=255)

                                                                      gTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATAT‑CTTACTTGAAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGa                                                                                                                                                                                    <  2:624416/275‑1 (MQ=255)

                                                                      gTTTTAACGTTAACTAATTAATATGAAATCAATTA‑AGAAAAATGTTTACTTGGAATCGGGTTGTAAAACGTTAGTATTAGCAGCCACGGTGCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGa                                                                                                                                                                                    <  4:627062/275‑1 (MQ=255)

                                                                       ttttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATca                                                                                                                                                                                                              <  4:373661/248‑1 (MQ=255)

                                                                                      attaattaTAAATCAGTTAGCGAAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCa                                                                                                                                                                                                                                              >  4:264354/1‑202 (MQ=255)

                                                                                          attaTAAATCAGTTAGCGAAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatta                                                                                                                                                                                                                                          <  2:264357/202‑1 (MQ=255)

                                                                                                                            aaTCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                                                                                                                                                                                    <  6:1393361/143‑1 (MQ=255)

                                                                                                                            aaTCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                                                                                                                                                                                    >  5:1393361/1‑143 (MQ=255)

                                                                                                                              tCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggtt                                                                                                                                                                                                                                                            <  5:1239993/149‑1 (MQ=255)

                                                                                                                                ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                                                                                                                                                                                          <  7:1437792/149‑1 (MQ=255)

                                                                                                                                 gtgtGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                                                                                                                                         <  7:1651611/149‑1 (MQ=255)

                                                                                                                                 gtgtGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                                                                                                                                         <  7:895122/149‑1 (MQ=255)

                                                                                                                                 gtgtGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                                                                                                                                         <  8:1350831/149‑1 (MQ=255)

                                                                                                                                 gtatgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                                                                                                                                         <  7:527147/146‑1 (MQ=255)

                                                                                                                                  tgtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                                                                                                                                        <  8:267128/149‑1 (MQ=255)

                                                                                                                                   gtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                                                                                                                                                                                       >  5:1104495/1‑149 (MQ=255)

                                                                                                                                    tgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                                                                                                                                                                                                      >  5:144303/1‑149 (MQ=255)

                                                                                                                                     ggAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttctt                                                                                                                                                                                                                                                     <  5:517801/149‑4 (MQ=255)

                                                                                                                                         aaCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCt                                                                                                                                                                                                                                                 >  8:444176/1‑149 (MQ=255)

                                                                                                                                            ggTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttt                                                                                                                                                                                                                                                           >  8:1260863/1‑136 (MQ=255)

                                                                                                                                            ggTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttt                                                                                                                                                                                                                                                           <  7:1260863/136‑1 (MQ=255)

                                                                                                                                            ggTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggt                                                                                                                                                                                                                                                                 <  1:3741/130‑1 (MQ=255)

                                                                                                                                             gtagtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                                                                                                                                        >  5:815260/1‑138 (MQ=255)

                                                                                                                                              tagtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTAtttt                                                                                                                                                                                               <  1:232545/193‑1 (MQ=255)

                                                                                                                                              tagtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCtgactgac                                                                                                             >  2:226113/1‑275 (MQ=255)

                                                                                                                                                        atccACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAACTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCATACCACACCCGCCTGAGTATTCCACGGCGATTCTCACCCTGCGCCATTGCCATAGCTGACTGACACATGAcctc                                                                                                   >  2:144200/3‑273 (MQ=255)

                                                                                                                                                          ccACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                                                                                                                                         >  5:1130397/1‑124 (MQ=255)

                                                                                                                                                           cACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGATCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATACCTGCAATTAATTCGCTAAATAGACCGTT‑CCACAACAAAGCCTGCTTATTTTCGGTTATCAGAACGTCCAGACCAAACCCAACCGAATATTACACGGCAATTCAACCCCTGCACCNNNTCCAAAGCCGACTGAAAAATAAcaaaacaaaaa                                                                                            >  3:263518/1‑269 (MQ=255)

                                                                                                                                                           cACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTGTATCTGCAATTAATTCGATAAAc                                                                                                                                                                                                                               <  6:415370/149‑1 (MQ=255)

                                                                                                                                                           cACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTNTATATCTGCAATTAATTCGATAANCAGANCCTG‑ACACATCAAAGCCTGTTTATTTTCTGTTATCAGCACGTCCAGACCacac                                                                                                                                                                     >  4:703317/1‑206 (MQ=255)

                                                                                                                                                            aCGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                                                                                                                                                                                               >  6:1652866/1‑148 (MQ=255)

                                                                                                                                                            aCGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                                                                                                                                                                                               <  5:1652866/148‑1 (MQ=255)

                                                                                                                                                              gAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGa                                                                                                                                                                                                                            <  7:841294/149‑1 (MQ=255)

                                                                                                                                                              gAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGa                                                                                                                                                                                                                            >  8:841294/1‑149 (MQ=255)

                                                                                                                                                               aGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc                                                                                                                                                                               <  2:703322/192‑1 (MQ=255)

                                                                                                                                                                  cGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                                                                                                                                     <  5:777189/136‑1 (MQ=255)

                                                                                                                                                                       cGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑acac                                                                                                                                                                                                                  >  8:1621526/1‑149 (MQ=255)

                                                                                                                                                                         tAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACAt                                                                                                                                                                                                                <  7:1621526/149‑1 (MQ=255)

                                                                                                                                                                            cgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGAt                                                                               <  2:430984/275‑1 (MQ=255)

                                                                                                                                                                                   tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGTCCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGt                                                                                                                                                                                                      <  6:482179/149‑1 (MQ=255)

                                                                                                                                                                                   tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGt                                                                                                                                                                                                      <  6:831464/149‑1 (MQ=255)

                                                                                                                                                                                       aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTAt                                                                                                                                                                                                  >  7:1562667/1‑149 (MQ=255)

                                                                                                                                                                                        gcgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTAtt                                                                                                                                                                                                 <  6:914056/149‑1 (MQ=255)

                                                                                                                                                                                            aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑a                                                                                                                                                                                                                     <  8:1662435/125‑1 (MQ=255)

                                                                                                                                                                                            aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCt                                                                                                                                                                                             <  6:104017/149‑1 (MQ=255)

                                                                                                                                                                                               gTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGc                                                                                                                                                                 <  1:102882/174‑1 (MQ=255)

                                                                                                                                                                                               gTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGc                                                                                                                                                                 >  3:102882/1‑174 (MQ=255)

                                                                                                                                                                                                tCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTa                                                                                                                                                                                         <  6:1380083/149‑1 (MQ=255)

                                                                                                                                                                                                tCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTa                                                                                                                                                                                         >  6:1616870/1‑149 (MQ=255)

                                                                                                                                                                                                 cATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTAt                                                                                                                                                                                        <  7:2912/149‑1 (MQ=255)

                                                                                                                                                                                                  aTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggtt                                                                                                                                                                                                                                                            <  6:1204537/81‑1 (MQ=38)

                                                                                                                                                                                                  aTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCt                                                                                                                                                                                             <  2:480634/143‑1 (MQ=255)

                                                                                                                                                                                                  aTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTCCTTTTTTc                                                                                                                                                                                       >  6:1254645/1‑149 (MQ=255)

                                                                                                                                                                                                          tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAACAACCAACCCt                                                                                                                                                                                                                                                                       >  4:461641/1‑62 (MQ=255)

                                                                                                                                                                                                          tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc                                                                                                                                                                               <  5:1549705/149‑1 (MQ=255)

                                                                                                                                                                                                          tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACGCCAATg                                                 >  2:58170/1‑275 (MQ=255)

                                                                                                                                                                                                           cGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc                                                                                                                                                                              <  5:1196327/149‑1 (MQ=255)

                                                                                                                                                                                                            gggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCg                                                                                                                                                                                                                         >  4:315339/1‑106 (MQ=255)

                                                                                                                                                                                                            gggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCa                                                                                                                                                                             >  6:784506/1‑149 (MQ=255)

                                                                                                                                                                                                            gggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCa                                                                                                                                                                             <  5:48885/149‑1 (MQ=255)

                                                                                                                                                                                                                      ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAccc                                                                                                                                                                   <  8:336605/149‑1 (MQ=255)

                                                                                                                                                                                                                       gTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCaaacttg                                                                                                                                                                  >  7:1338444/1‑146 (MQ=255)

                                                                                                                                                                                                                        ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGc                                                                                                                                                                 >  7:20212/1‑149 (MQ=255)

                                                                                                                                                                                                                        ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGc                                                                                                                                                                 >  6:305855/1‑149 (MQ=255)

                                                                                                                                                                                                                           aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                                                                                                                                                                                    <  6:957415/64‑1 (MQ=255)

                                                                                                                                                                                                                            aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGa                                                                                                                                                             >  5:455482/1‑149 (MQ=255)

                                                                                                                                                                                                                            aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGa                                                                                                                                                             <  8:1085271/149‑1 (MQ=255)

                                                                                                                                                                                                                                ctctCGTGCCGACCAAAAATACCAAGAAAACACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTAt                                                                                                                                                         <  6:1116318/149‑1 (MQ=255)

                                                                                                                                                                                                                                 tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGAc                                                                                                                                                                                                                           >  5:241289/1‑83 (MQ=255)

                                                                                                                                                                                                                                 tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTAtt                                                                                                                                                        >  7:1195332/1‑149 (MQ=255)

                                                                                                                                                                                                                                 tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTAtt                                                                                                                                                        >  8:1226906/1‑149 (MQ=255)

                                                                                                                                                                                                                                 tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTAtt                                                                                                                                                        >  6:1155044/1‑149 (MQ=255)

                                                                                                                                                                                                                                   tcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTcc                                                                                                                                                      >  5:1352001/1‑149 (MQ=255)

                                                                                                                                                                                                                                     gTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCAc                                                                                                                                                    >  8:58075/1‑149 (MQ=255)

                                                                                                                                                                                                                                            ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGAtt                                                                                                                                             >  5:957061/1‑149 (MQ=255)

                                                                                                                                                                                                                                               aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGAtt                                                                                                                                             >  7:416511/1‑146 (MQ=255)

                                                                                                                                                                                                                                                 aaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCt                                                                                                                                                               <  6:1424233/126‑1 (MQ=255)

                                                                                                                                                                                                                                                     ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTg                                                                                                                                    >  8:1593500/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                                                                                                                                                                                  <  6:531844/102‑1 (MQ=255)

                                                                                                                                                                                                                                                      cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCa                                                                                                                                                     >  7:455848/1‑131 (MQ=255)

                                                                                                                                                                                                                                                      cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgc                                                                                                                                   <  5:942831/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       aaGAAAAAACCAACCCTTATGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgcg                                                                                                                                  >  7:992203/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          aaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGCG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa                                                                                                                               >  8:707703/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          aaaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa                                                                                                                               >  5:779136/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          aaaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa                                                                                                                               >  5:180534/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            aaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAtt                                                                                                                             >  6:1542736/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            aaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAtt                                                                                                                             <  5:1026508/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              aaccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTTCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGc                                                                                                                           >  5:1371286/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGcc                                                                                                                          <  5:305855/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGcc                                                                                                                          <  7:1064994/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 caaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCAt                                                                                                                        >  6:1718741/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  aaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATa                                                                                                                       >  8:31671/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    cccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACg                                                       >  2:560747/1‑211 (MQ=255)

                                                                                                                                                                                                                                                                     ccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGct                                                                                                                    >  6:686772/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      cTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGctg                                                                                                                   >  5:657922/1‑149 (MQ=255)

CGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATAT‑CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTG‑ACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGACAAATGACCGCTTAACGTCGGGGTAAACGATGACCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCAC  >  NC_000913/2286050‑2286572

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