BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	224,893	T→A	noncoding (1123/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,893	0	T	A	89.1%	42.7 / 85.9	304	noncoding (1123/1542 nt)	rrsH	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (18/15); new base A (107/164); total (125/179)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 311 total aligned reads displayed.

TTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATA  >  NC_000913/224761‑225018

ttGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGACCGGCAACGAGCGCAACCCTTATCTTTTGTTGCCa                                                                                                                         <  2:1322585/140‑1 (MQ=255)

             aaCTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTccggccgg                                                                                                           <  2:1500765/140‑1 (MQ=255)

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

                                 gTGCCTTCGGGAACTGCGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCa                                                                                       >  2:1492743/1‑140 (MQ=21)

                                    ccTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTg                                                                                    <  1:1757734/140‑1 (MQ=34)

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

                                               ttcaagaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCggag                                                                         <  1:564328/136‑1 (MQ=2)

                                                tcaagaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCggagg                                                                        <  1:835615/137‑1 (MQ=2)

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

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

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

                                                 caagaCAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCggagga                                                                       <  1:226366/138‑1 (MQ=2)

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

                                                       aGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTg                                                                  <  1:586664/140‑1 (MQ=2)

                                                         gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAAc                                                                              >  2:66936/5‑125 (MQ=2)

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

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

                                                         gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  2:1219998/5‑140 (MQ=2)

                                                         gtggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTgggg                                                               >  2:514722/5‑140 (MQ=2)

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

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

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

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

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

                                                          tggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCGGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGa                                                              >  2:368521/4‑140 (MQ=2)

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

                                                          tggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGa                                                              >  2:476361/4‑140 (MQ=2)

                                                           ggtgcATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGAt                                                             <  2:17865/138‑1 (MQ=2)

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

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

                                                               cATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGAc                                                          <  2:1083427/140‑1 (MQ=2)

                                                                aTGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACg                                                         <  2:537190/140‑1 (MQ=2)

                                                                  ggttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCa                                                      <  1:1220273/137‑1 (MQ=2)

                                                                   gttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                     <  1:87241/138‑1 (MQ=11)

                                                                    ttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaa                                                     >  2:733351/2‑140 (MQ=2)

                                                                    ttGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaaa                                                    <  2:1069943/139‑2 (MQ=2)

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

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

                                                                     tGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCaaa                                                    <  1:647495/140‑2 (MQ=2)

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

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

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

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

                                                                       tcgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtca                                                 >  2:1022787/1‑140 (MQ=11)

                                                                        cgtcAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCAGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAAtcat                                                <  1:669422/140‑1 (MQ=2)

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

                                                                            aGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGTGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                                            >  1:520935/1‑140 (MQ=2)

                                                                            aGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATg                                            >  2:568283/1‑140 (MQ=2)

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

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

                                                                              cTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgc                                           >  1:86000/1‑139 (MQ=2)

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

                                                                              cTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATgcc                                          >  1:458/1‑138 (MQ=2)

                                                                                 gtgtCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCt                                       >  1:386029/1‑140 (MQ=2)

                                                                                 gtgtCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCt                                       >  2:404948/1‑140 (MQ=2)

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

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

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

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

                                                                                       tgAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGAc                                 <  2:705732/140‑1 (MQ=2)

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

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

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

                                                                                         agaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAGTCATCATGCCCCTTATGAcct                               <  1:410015/138‑2 (MQ=2)

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

                                                                                              ttGGGTTAAGTCCCCCAACGAGCGCACCCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGc                           <  2:1165549/140‑1 (MQ=2)

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

                                                                                                gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGTGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  1:321697/140‑1 (MQ=2)

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

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

                                                                                                gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                        <  1:914195/140‑1 (MQ=2)

                                                                                                 ggTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTaca                       >  2:1602805/1‑140 (MQ=32)

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

                                                                                                  gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                      <  1:788965/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‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCGAATCATCATGCCCCTTATGACCTGGGCTacaca                     >  1:430534/1‑140 (MQ=2)

                                                                                                     aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCATATCATCATGCCCCTTATGACCTGGGCTACACACg                   >  2:53356/1‑140 (MQ=2)

                                                                                                     aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACg                   >  1:196622/1‑140 (MQ=2)

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

                                                                                                      aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGc                           <  1:626431/131‑1 (MQ=2)

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

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

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

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

                                                                                                        tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                         >  2:1116951/1‑132 (MQ=2)

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

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

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

                                                                                                            gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  2:342908/140‑1 (MQ=2)

                                                                                                            gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  1:71439/140‑1 (MQ=2)

                                                                                                             cAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACaa           >  2:394804/1‑140 (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‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  1:1131151/140‑1 (MQ=2)

                                                                                                                  ggcgaAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCGAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGgcg      <  2:1243256/139‑1 (MQ=12)

                                                                                                                  aGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:294938/140‑1 (MQ=2)

                                                                                                                      cAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATa  <  2:1410256/140‑1 (MQ=32)

TTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATA  >  NC_000913/224761‑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: `‑`