BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,211,614	C→T	6.7%	intergenic (+37/+66)	icdC → / ← iraM	pseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion/RpoS stabilzer during Mg starvation, anti‑RssB factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,211,614	0	C	T	6.7%	480.4 / 17.1	149	intergenic (+37/+66)	icdC/iraM	pseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion/RpoS stabilzer during Mg starvation, anti‑RssB factor
Reads supporting (aligned to +/- strand): ref base C (73/66); new base T (6/4); total (79/70)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-02

Only 100 of 154 total aligned reads displayed.

TACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCTTTTGAGTAAACTCCATAAATCCTTGTTGAATGGTGCGATGTGATAAATAGTAATAGGATATTCTTTATCCTT  >  NC_000913/1211392‑1211847

tACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATattatt                                                                                                                                                                                                           >  2:560466/1‑255 (MQ=255)

                                        ggACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGTCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGa                                                                                                                                                                                >  2:575940/1‑242 (MQ=255)

                                                   gCTGACCTGATTGTTAAAGGTAGGGAAGGCGAAATCAATGAAAAGATCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATg                                                                                                                                                                                 <  1:92376/230‑1 (MQ=255)

                                                                  aaaGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                                                                                                                                                                       >  1:35437/1‑177 (MQ=255)

                                                                              ggCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgtt                                                                                                                                                                                                                               <  4:48785/157‑1 (MQ=255)

                                                                                cgcAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGccc                                                                                                                                                                                                                                     <  5:159654/149‑1 (MQ=255)

                                                                                 gcAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAGTATTGAAACGGGCGTATAACAGGACCg                                                                                                                                                                                                                                    >  8:933817/1‑149 (MQ=255)

                                                                                 gcAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCg                                                                                                                                                                                                                                    >  5:888/1‑149 (MQ=255)

                                                                                 gcAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCg                                                                                                                                                                                                                                    >  8:597535/1‑149 (MQ=255)

                                                                                   aatcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGc                                                                                                                                                                                                                                       >  8:270090/1‑144 (MQ=255)

                                                                                   aatcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGc                                                                                                                                                                                                                                       <  7:270090/144‑1 (MQ=255)

                                                                                   aatcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgtt                                                                                                                                                                                                                                  >  5:1032354/1‑149 (MQ=255)

                                                                                    atcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTGTTATACGCCCgttt                                                                                                                                                                                                                                 >  8:51162/1‑148 (MQ=255)

                                                                                    atcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttg                                                                                                                                                                                                                                 >  7:789838/1‑149 (MQ=255)

                                                                                     tcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgt                                                                                                                                                                                                                                <  6:341943/149‑1 (MQ=255)

                                                                                     tcaatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCgttgt                                                                                                                                                                                                                                >  6:177139/1‑149 (MQ=255)

                                                                                       aatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttt                                                                                                                                                                                                                              <  7:864204/149‑1 (MQ=255)

                                                                                       aatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGCGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttt                                                                                                                                                                                                                              <  6:855906/149‑1 (MQ=255)

                                                                                       aatGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTCCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttt                                                                                                                                                                                                                              <  6:1065631/149‑1 (MQ=255)

                                                                                          gCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttat                                                                                                                                                                                                                           <  6:649987/149‑1 (MQ=255)

                                                                                            cAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattt                                                                                                                                                                                                                         >  7:432666/1‑149 (MQ=255)

                                                                                             aaGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTtttattta                                                                                                                                                                                                                        >  5:118732/1‑149 (MQ=255)

                                                                                              aGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTttatttat                                                                                                                                                                                                                       <  6:593125/149‑1 (MQ=255)

                                                                                                aCCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGACATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgt                                                                                                                                                                                                                     <  6:921176/149‑1 (MQ=255)

                                                                                                 ccGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTAtgtg                                                                                                                                                                                                                    <  8:10737/149‑1 (MQ=255)

                                                                                                        ttATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtatta                                                                                                                                                                                                             <  6:888/149‑1 (MQ=255)

                                                                                                          aTGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGAtatt                                                                                                                                                                                                              >  4:653522/1‑146 (MQ=255)

                                                                                                            gACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTaa                                                                                                                                                                                                         >  6:507983/1‑149 (MQ=255)

                                                                                                            gACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTaa                                                                                                                                                                                                         <  5:909784/149‑1 (MQ=255)

                                                                                                              cTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATAttatta                                                                                                                                                                                                          <  2:653522/146‑1 (MQ=255)

                                                                                                               ttCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATTTTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAg                                                                                                                                                                                                      >  5:131543/1‑149 (MQ=255)

                                                                                                               ttCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAg                                                                                                                                                                                                      >  7:126586/1‑149 (MQ=255)

                                                                                                               ttCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAg                                                                                                                                                                                                      >  7:479440/1‑149 (MQ=255)

                                                                                                                 cGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCatat                                                                                                                                                                                                 >  4:190503/1‑152 (MQ=255)

                                                                                                                        cTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAg                                                                                                                                                                                             >  5:829391/1‑149 (MQ=255)

                                                                                                                         tGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGc                                                                                                                                                                                            >  6:47125/1‑149 (MQ=255)

                                                                                                                           aTGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCat                                                                                                                                                                                          >  7:974156/1‑149 (MQ=255)

                                                                                                                            tGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCata                                                                                                                                                                                         <  6:699322/149‑1 (MQ=255)

                                                                                                                            tGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCata                                                                                                                                                                                         >  7:782500/1‑149 (MQ=255)

                                                                                                                                 ggCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTtat                                                                                                                                                                                    >  6:164909/1‑149 (MQ=255)

                                                                                                                                   cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTatat                                                                                                                                                                                  <  5:177139/149‑1 (MQ=255)

                                                                                                                                           tgctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCa                                                                                                                                                                          >  8:19164/1‑149 (MQ=255)

                                                                                                                                             ctgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTtata                                                                                                                                           >  2:455785/1‑178 (MQ=255)

                                                                                                                                                 aaTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGc                                                                                                                                       <  4:455785/178‑1 (MQ=255)

                                                                                                                                                 aaTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACAGGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACtt                                                                                                                                                                    <  8:907960/149‑1 (MQ=255)

                                                                                                                                                  aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTg                                                                                                                                                                   >  6:802145/1‑149 (MQ=255)

                                                                                                                                                        agagTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCa                                                                                                                                                             <  7:1062232/149‑1 (MQ=255)

                                                                                                                                                         gagTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCat                                                                                                                                                            >  8:236882/1‑149 (MQ=255)

                                                                                                                                                               ggTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGc                                                                                                                                                      <  5:945164/149‑1 (MQ=255)

                                                                                                                                                                  gAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCata                                                                                                                                                   >  8:657380/1‑149 (MQ=255)

                                                                                                                                                                   aaGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCatat                                                                                                                                                  >  5:594616/1‑149 (MQ=255)

                                                                                                                                                                     gCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCatat                                                                                                                                                          <  7:56394/139‑1 (MQ=255)

                                                                                                                                                                       gATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCTATTACTTGAGCCCATATGGGCATAttttt                                                                                                                                              >  8:1014748/1‑149 (MQ=255)

                                                                                                                                                                        atgatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGATCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta                                                                                                                                             >  8:724106/4‑149 (MQ=255)

                                                                                                                                                                        atcatcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta                                                                                                                                             >  7:39389/1‑149 (MQ=255)

                                                                                                                                                                           atcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATa                                                                                                                                                                                                       >  8:976621/1‑88 (MQ=255)

                                                                                                                                                                           atcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATa                                                                                                                                                                                                       <  7:976621/88‑1 (MQ=255)

                                                                                                                                                                            tcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAAt                                                                                                                                         <  6:826430/149‑1 (MQ=255)

                                                                                                                                                                            tcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAAt                                                                                                                                         <  8:852754/149‑1 (MQ=255)

                                                                                                                                                                            tcGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAAt                                                                                                                                         >  5:558726/1‑149 (MQ=255)

                                                                                                                                                                             cGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATg                                                                                                                                        <  8:782500/149‑1 (MQ=255)

                                                                                                                                                                             cGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATg                                                                                                                                        >  7:608566/1‑149 (MQ=255)

                                                                                                                                                                               aaaaCATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCTTTTGAGTAAACTCAATAAATCCTTGTTGAATGGTGGTATGTGATAAAAAGTAATAGGAAAtttttt        >  4:635094/1‑275 (MQ=255)

                                                                                                                                                                                     tGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACtatt                                                                                                                                >  5:313097/1‑149 (MQ=255)

                                                                                                                                                                                      gTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACtatta                                                                                                                               >  8:209516/1‑149 (MQ=255)

                                                                                                                                                                                       tAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACtattat                                                                                                                              <  6:725921/149‑1 (MQ=255)

                                                                                                                                                                                         aTCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGt                                                                                                                            >  5:845612/1‑149 (MQ=255)

                                                                                                                                                                                            tctcCATGTTTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTaaa                                                                                                                         <  8:432666/149‑1 (MQ=255)

                                                                                                                                                                                               ccATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACat                                                                                                                      <  5:394876/149‑1 (MQ=255)

                                                                                                                                                                                                cATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTtata                                                                                                                                                                                   >  5:730011/1‑87 (MQ=255)

                                                                                                                                                                                                 aTGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACattt                                                                                                                    >  5:606812/1‑149 (MQ=255)

                                                                                                                                                                                                 aTGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACattt                                                                                                                    >  5:780453/1‑149 (MQ=255)

                                                                                                                                                                                                  tgtgTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACattta                                                                                                                   <  8:68463/149‑1 (MQ=255)

                                                                                                                                                                                                    tgtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACatttatt                                                                                                                 >  5:503704/1‑149 (MQ=255)

                                                                                                                                                                                                    tgtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACatttatt                                                                                                                 <  5:228777/149‑1 (MQ=255)

                                                                                                                                                                                                     gtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACatttattt                                                                                                                <  5:1017264/149‑1 (MQ=255)

                                                                                                                                                                                                     gtTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACatttattt                                                                                                                <  7:1014748/149‑1 (MQ=255)

                                                                                                                                                                                                      ttAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTTTTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTg                                                                                                               <  7:933817/149‑1 (MQ=255)

                                                                                                                                                                                                      ttAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCata                                                                                                                                                                                                  >  1:205673/1‑66 (MQ=255)

                                                                                                                                                                                                      ttAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTg                                                                                                               >  8:968673/1‑149 (MQ=255)

                                                                                                                                                                                                        aaaTATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGtt                                                                                                             >  7:605097/1‑149 (MQ=255)

                                                                                                                                                                                                          atatTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGa                                                                                                                                                                                              <  3:205673/66‑1 (MQ=255)

                                                                                                                                                                                                           tatTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACatttat                                                                                                                  >  5:968897/1‑141 (MQ=255)

                                                                                                                                                                                                           tatTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACatttat                                                                                                                  <  6:968897/141‑1 (MQ=255)

                                                                                                                                                                                                           tatTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTAtt                                                                                                          >  6:284655/1‑149 (MQ=255)

                                                                                                                                                                                                            atTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTAttt                                                                                                         <  6:200010/149‑1 (MQ=255)

                                                                                                                                                                                                            atTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTAttt                                                                                                         <  6:135953/149‑1 (MQ=255)

                                                                                                                                                                                                             ttGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTAtttt                                                                                                        >  6:1049851/1‑149 (MQ=255)

                                                                                                                                                                                                              tGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCatat                                                                                                                                                                                                 >  2:121833/1‑59 (MQ=255)

                                                                                                                                                                                                              aacaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTg                                                                                                       >  6:537839/4‑149 (MQ=255)

                                                                                                                                                                                                               gAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGc                                                                                                      <  5:520497/149‑1 (MQ=255)

                                                                                                                                                                                                               acaaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGc                                                                                                      >  6:1021487/3‑149 (MQ=255)

                                                                                                                                                                                                                caaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCt                                                                                                     >  8:723345/2‑149 (MQ=255)

                                                                                                                                                                                                                 aaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttt                                                                                                                                              <  8:634772/107‑1 (MQ=255)

                                                                                                                                                                                                                 aaCGGGCGTGTTATACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATAttttt                                                                                                                                              >  7:634772/1‑107 (MQ=255)

                                                                                                                                                                                                                         tataACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTta                                                                                                                                             <  2:458516/100‑1 (MQ=255)

                                                                                                                                                                                                                            aaCACGCCCGTTGTTTTATTTATGTGGATATTATTACTAGCATATCGAGCATATTTATATGAAGCCGATTACTTGAGCCCATATGGGTATATTTTTATGATGCAACTAATATGTAAACATTTATTTGTTATTTTGCTTTCCCCTggagg                                                                                         >  7:673121/1‑149 (MQ=255)

                                                                                                                                                                                                                             acacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTggagga                                                                                        >  7:89652/1‑149 (MQ=255)

                                                                                                                                                                                                                             acacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTggagga                                                                                        >  5:165418/1‑149 (MQ=255)

                                                                                                                                                                                                                              cacGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGac                                                                                       >  7:277949/1‑149 (MQ=255)

TACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCTTTTGAGTAAACTCCATAAATCCTTGTTGAATGGTGCGATGTGATAAATAGTAATAGGATATTCTTTATCCTT  >  NC_000913/1211392‑1211847

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