BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	224,920	A→T	noncoding (1150/1542 nt)	rrsH →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	224,920	0	A	T	90.4%	58.9 / 73.6	301	noncoding (1150/1542 nt)	rrsH	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (17/12); new base T (107/165); total (124/177)
Fisher's exact test for biased strand distribution p-value = 4.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

Only 100 of 307 total aligned reads displayed.

GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATA  >  NC_000913/224789‑225018

gATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGAc                                                                                            >  1:1588311/1‑140 (MQ=40)

 aTTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACt                                                                                           >  1:1752143/1‑140 (MQ=39)

                   ttcaagaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCggag                                                                         >  1:1142337/5‑140 (MQ=2)

                    gtgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggagg                                                                        <  2:1778473/140‑1 (MQ=32)

                     tgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggagga                                                                       <  1:1704560/140‑1 (MQ=32)

                     tgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTggagga                                                                       >  2:1651632/1‑140 (MQ=32)

                     caagaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCggag                                                                         >  1:728820/3‑138 (MQ=2)

                        gaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAgg                                                                    <  2:616344/140‑1 (MQ=2)

                             gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTggg                                                                >  1:789876/5‑140 (MQ=2)

                             gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:1211099/5‑140 (MQ=2)

                             gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  2:951339/5‑140 (MQ=2)

                             gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCggagg                                                                        >  2:87623/5‑131 (MQ=2)

                             gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:1184577/5‑140 (MQ=2)

                             gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:199144/5‑140 (MQ=2)

                             gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTCGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:607671/5‑140 (MQ=2)

                             gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCggagga                                                                       >  1:864143/1‑133 (MQ=2)

                             gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  2:765329/1‑140 (MQ=2)

                             gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               <  2:1205636/140‑1 (MQ=2)

                             gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               <  1:283631/140‑1 (MQ=2)

                             gtggtgCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  1:178076/1‑140 (MQ=2)

                             gTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTgggg                                                               >  2:1497198/1‑140 (MQ=32)

                              tggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGa                                                              <  2:934162/137‑1 (MQ=2)

                               ggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGAt                                                             <  1:110111/138‑1 (MQ=2)

                               ggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGAt                                                             <  2:1170816/138‑1 (MQ=2)

                                gtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGAt                                                             <  1:516475/139‑1 (MQ=2)

                                gtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATg                                                            >  1:960082/2‑140 (MQ=2)

                                  gcATGGTTGTCGTCAGCTCGTGTCGTAAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGAc                                                          <  2:219871/140‑1 (MQ=2)

                                   cATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACg                                                         <  2:570152/140‑1 (MQ=2)

                                       gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                     >  1:778466/3‑140 (MQ=2)

                                       gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                     <  2:129289/138‑1 (MQ=2)

                                        ttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaaa                                                    <  2:1184577/139‑2 (MQ=11)

                                        cTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAg                                                    >  2:1721801/1‑140 (MQ=32)

                                           tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGa                                                              >  1:850552/1‑128 (MQ=2)

                                           tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtc                                                  >  1:902450/1‑140 (MQ=2)

                                           tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                 >  2:806782/1‑140 (MQ=2)

                                           tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                 >  1:644917/1‑140 (MQ=2)

                                           tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                 >  2:422996/1‑140 (MQ=2)

                                           tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGa                                                              >  2:477411/1‑127 (MQ=2)

                                           tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                 >  1:1166517/1‑140 (MQ=2)

                                           tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                 >  1:1231169/1‑140 (MQ=2)

                                            cgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtcat                                                <  2:676606/140‑1 (MQ=2)

                                              tcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtcatca                                              <  1:1012359/140‑1 (MQ=2)

                                                aGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATcatcat                                             >  1:583785/1‑140 (MQ=2)

                                                aGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                                            <  1:422996/140‑1 (MQ=2)

                                                aGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                                            <  1:1219998/140‑1 (MQ=2)

                                                 gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           <  2:318631/140‑2 (MQ=2)

                                                 gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           <  2:1142337/140‑2 (MQ=2)

                                                 gCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           <  1:695428/140‑2 (MQ=2)

                                                  cTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           <  1:451475/140‑2 (MQ=2)

                                                      tttcgtgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtcatca                                              <  2:1129335/128‑1 (MQ=2)

                                                      tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  2:700541/140‑1 (MQ=2)

                                                      tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  1:34586/140‑1 (MQ=2)

                                                      tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  2:657762/140‑1 (MQ=2)

                                                      tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  1:476361/140‑1 (MQ=2)

                                                      tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  1:747891/140‑1 (MQ=2)

                                                      tgtcgtGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  1:537028/140‑1 (MQ=2)

                                                       gtcgtgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTa                                     >  1:665025/4‑140 (MQ=2)

                                                          gtgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaaat                                                   <  1:257526/123‑3 (MQ=2)

                                                           tgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  1:1036868/140‑1 (MQ=2)

                                                           tgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  1:1201321/140‑1 (MQ=2)

                                                           tgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  2:737499/140‑1 (MQ=2)

                                                           tgAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAc                                 <  1:1271751/140‑1 (MQ=32)

                                                            gagaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  2:710837/137‑1 (MQ=2)

                                                            gagaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           >  2:1024654/4‑128 (MQ=2)

                                                            gagaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcc                                <  1:15439/137‑1 (MQ=2)

                                                             agaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcct                               <  2:458/138‑2 (MQ=2)

                                                             agaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcct                               <  2:24334/138‑2 (MQ=2)

                                                              gaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcctg                              <  2:533615/139‑3 (MQ=2)

                                                              gaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcctg                              <  2:553322/139‑3 (MQ=2)

                                                                  ttGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAAACGGCGGGAGGTGGGGAAGACGTCAAATCATCATGCCCCTTATGACCTGGGCt                          >  2:592137/1‑140 (MQ=2)

                                                                   tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCtt                                      <  1:970694/128‑1 (MQ=2)

                                                                   tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                         <  1:831289/140‑1 (MQ=2)

                                                                    gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:1231169/140‑1 (MQ=2)

                                                                     ggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:832369/140‑1 (MQ=2)

                                                                      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  2:1097507/140‑1 (MQ=2)

                                                                      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCTGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:502462/140‑1 (MQ=2)

                                                                      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  2:14476/140‑1 (MQ=2)

                                                                      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:404948/140‑1 (MQ=2)

                                                                      gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  2:177081/140‑1 (MQ=2)

                                                                       ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:778466/140‑1 (MQ=2)

                                                                       ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:644917/140‑1 (MQ=2)

                                                                       ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:515424/140‑1 (MQ=2)

                                                                        tAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                     <  2:1004608/140‑1 (MQ=2)

                                                                         aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcc                                <  1:139174/127‑1 (MQ=2)

                                                                         aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAcc                                >  2:139174/1‑127 (MQ=2)

                                                                          aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:356619/135‑1 (MQ=2)

                                                                          aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:1048760/134‑1 (MQ=2)

                                                                           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 <  2:910073/140‑2 (MQ=2)

                                                                           gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 >  2:277854/1‑139 (MQ=2)

                                                                            tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                 >  2:374426/1‑139 (MQ=2)

                                                                            tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 >  2:1076668/1‑123 (MQ=2)

                                                                            tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  1:1076668/123‑1 (MQ=2)

                                                                               cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:1001282/131‑1 (MQ=2)

                                                                                gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  1:372530/140‑1 (MQ=2)

                                                                                    cGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATgg        <  1:1497198/140‑1 (MQ=32)

                                                                                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:869305/140‑1 (MQ=2)

                                                                                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:977036/140‑1 (MQ=2)

                                                                                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:466811/140‑1 (MQ=2)

                                                                                     gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  1:1131151/140‑1 (MQ=2)

                                                                                      aGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  2:401745/122‑1 (MQ=2)

GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATA  >  NC_000913/224789‑225018

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