BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,942,018	G→T	noncoding (211/1542 nt)	rrsC →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,942,018	0	G	T	99.4%	359.5 / ‑3.9	155	noncoding (211/1542 nt)	rrsC	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (1/0); new base T (71/83); total (72/83)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 159 total aligned reads displayed.

GAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACG  >  NC_000913/3941883‑3942153

gAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGcc                                                                                                                                     >  2:1301939/1‑140 (MQ=17)

  aaCAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCctc                                                                                                                                   >  2:17896/1‑140 (MQ=17)

  aaCAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCctc                                                                                                                                   <  2:1977365/140‑1 (MQ=18)

  aaCAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCctc                                                                                                                                   <  1:1131822/140‑1 (MQ=18)

  aaCAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCctc                                                                                                                                   >  1:495142/1‑140 (MQ=17)

    cAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCtt                                                                                                                                 <  1:2184063/140‑1 (MQ=18)

    cAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCtt                                                                                                                                 <  1:1428471/140‑1 (MQ=18)

         tGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCAt                                                                                                                            <  1:1307121/140‑1 (MQ=18)

         tGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCAt                                                                                                                            <  1:570053/140‑1 (MQ=18)

         tGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCAt                                                                                                                            <  2:2218314/140‑1 (MQ=18)

          gCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATc                                                                                                                           >  2:1773478/1‑140 (MQ=17)

          gCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATc                                                                                                                           >  1:995299/1‑140 (MQ=17)

               ttCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGAtg                                                                                                                      >  1:330118/1‑140 (MQ=21)

               ttCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGAtg                                                                                                                      >  1:911107/1‑140 (MQ=21)

                 cGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGAtgtg                                                                                                                    <  2:66737/140‑1 (MQ=18)

                      aCGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAg                                                                                                               >  1:1283113/1‑140 (MQ=17)

                      aCGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAg                                                                                                               <  2:1626255/140‑1 (MQ=18)

                       cGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGTCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGa                                                                                                              <  1:1726371/140‑1 (MQ=16)

                       cGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTTCTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGa                                                                                                              <  1:1160338/140‑1 (MQ=12)

                       cGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGa                                                                                                              <  2:1142431/140‑1 (MQ=18)

                       cGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGa                                                                                                              <  2:1283113/140‑1 (MQ=17)

                             gCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGAt                                                                                                        <  2:468085/140‑1 (MQ=18)

                             gCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGAt                                                                                                        <  1:848641/140‑1 (MQ=18)

                                gACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAg                                                                                                     >  2:304572/1‑140 (MQ=17)

                                  cGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCt                                                                                                   <  1:1037210/140‑1 (MQ=18)

                                  cGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCt                                                                                                   >  1:1112207/1‑140 (MQ=17)

                                  cGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCt                                                                                                   >  1:1846236/1‑140 (MQ=17)

                                     gTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAg                                                                                                               <  1:897088/125‑1 (MQ=18)

                                        aGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAgg                                                                                             >  1:818662/1‑140 (MQ=17)

                                        aGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAgg                                                                                             >  2:1864048/1‑140 (MQ=17)

                                        aGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAgg                                                                                             <  2:2076330/140‑1 (MQ=17)

                                              gTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGt                                                                                       >  2:1876242/1‑140 (MQ=17)

                                              gTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGt                                                                                       >  2:788746/1‑140 (MQ=17)

                                              gTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGt                                                                                       <  2:2053265/140‑1 (MQ=18)

                                               tCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTa                                                                                      <  2:538182/140‑1 (MQ=18)

                                               tCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTa                                                                                      >  1:1874104/1‑140 (MQ=17)

                                               tCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTa                                                                                      >  1:521921/1‑140 (MQ=17)

                                                  gggAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCAT‑ACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACgg                                                                                  >  1:479951/1‑140 (MQ=12)

                                                        cTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGTTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCAc                                                                             >  1:1236213/1‑140 (MQ=17)

                                                         tGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCAcc                                                                            <  1:89049/140‑1 (MQ=17)

                                                         tGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCAcc                                                                            <  1:1940299/140‑1 (MQ=17)

                                                         tGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCAcc                                                                            <  1:1941217/140‑1 (MQ=17)

                                                                tGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcga                                                                     <  2:1924830/140‑1 (MQ=17)

                                                                tGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcga                                                                     <  1:699139/140‑1 (MQ=17)

                                                                tGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcga                                                                     <  2:151356/140‑1 (MQ=17)

                                                                tGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcga                                                                     <  1:304572/140‑1 (MQ=17)

                                                                tGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcga                                                                     <  1:27487/140‑1 (MQ=17)

                                                                tGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCCAGGcga                                                                     <  2:661877/140‑1 (MQ=17)

                                                                 ggAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcgac                                                                    <  1:85813/140‑1 (MQ=17)

                                                                  gAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcgacg                                                                   <  1:1876242/140‑1 (MQ=17)

                                                                   aGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcgacga                                                                  >  2:74314/1‑140 (MQ=17)

                                                                     ggggATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATc                                                                >  2:1066297/1‑140 (MQ=17)

                                                                     ggggATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATc                                                                >  1:1063459/1‑140 (MQ=17)

                                                                     ggggATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATc                                                                >  1:1552760/1‑140 (MQ=17)

                                                                         aTAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTa                                                            >  1:796131/1‑140 (MQ=17)

                                                                          tAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAg                                                           >  2:113739/1‑140 (MQ=17)

                                                                           aaCTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAga                                                          <  2:716457/140‑2 (MQ=17)

                                                                             ctactGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTg                                                        >  1:2005518/1‑140 (MQ=17)

                                                                             ctactGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTg                                                        >  1:649513/1‑140 (MQ=17)

                                                                                 tGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCt                                                    >  2:599506/1‑140 (MQ=17)

                                                                                 tGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCt                                                    <  1:559065/140‑1 (MQ=17)

                                                                                  ggAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcgac                                                                    <  2:881900/123‑1 (MQ=17)

                                                                                  ggAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTg                                                   <  1:223543/140‑1 (MQ=17)

                                                                                    aaaCGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTgag                                                 <  1:1326666/140‑1 (MQ=17)

                                                                                        ggTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGa                                             >  1:1675330/1‑140 (MQ=17)

                                                                                          tAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATg                                           <  2:22368/140‑1 (MQ=17)

                                                                                          tAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATg                                           >  2:1527889/1‑140 (MQ=17)

                                                                                                aTACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                     <  2:1216339/140‑1 (MQ=17)

                                                                                                aTACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                     <  1:542730/140‑1 (MQ=17)

                                                                                                aTACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                     <  1:1593510/140‑1 (MQ=17)

                                                                                                      cATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTTGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACt                               >  1:2094181/1‑140 (MQ=17)

                                                                                                      cATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                     <  1:79979/134‑1 (MQ=17)

                                                                                                      cATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACt                               >  2:826029/1‑140 (MQ=17)

                                                                                                      cATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACt                               >  2:2125145/1‑140 (MQ=17)

                                                                                                      cATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACt                               >  1:1197313/1‑140 (MQ=17)

                                                                                                      cATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACt                               >  1:1827188/1‑140 (MQ=17)

                                                                                                      cATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACt                               >  2:916615/1‑140 (MQ=17)

                                                                                                       aTAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTg                              <  2:1552760/140‑1 (MQ=17)

                                                                                                          aCGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGaa                           >  2:132743/1‑140 (MQ=17)

                                                                                                           cGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAAc                          >  2:1110879/1‑140 (MQ=17)

                                                                                                           cGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAAc                          >  2:1357633/1‑140 (MQ=17)

                                                                                                           cGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGCGGGGCAACGGCTCACCTAGGCGACGATCCCTAGCTGATCTGAGAGGACGACCAGCCACACCGGGAc                          >  2:1220377/1‑140 (MQ=11)

                                                                                                             tCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCacac                                >  2:2024329/1‑132 (MQ=17)

                                                                                                             tCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCacac                                <  1:2024329/132‑1 (MQ=17)

                                                                                                                 aaGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGAGGGGATTAGCTAGTAGGTGGGGTAAGGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGac                    <  2:2228086/140‑1 (MQ=17)

                                                                                                                  aGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                   <  2:494460/140‑1 (MQ=17)

                                                                                                                  aGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                   <  1:1463421/140‑1 (MQ=17)

                                                                                                                  aGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                   <  1:1864048/140‑1 (MQ=17)

                                                                                                                  aGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                   <  2:1675330/140‑1 (MQ=17)

                                                                                                                  aGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                   <  1:423159/140‑1 (MQ=17)

                                                                                                                  aGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                   <  1:1728793/140‑1 (MQ=17)

                                                                                                                        aaGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTcc             >  2:1254377/1‑140 (MQ=18)

                                                                                                                        aaGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTcc             <  1:1277519/140‑1 (MQ=17)

                                                                                                                        aaGAGGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTcc             >  1:977988/1‑140 (MQ=17)

                                                                                                                         agagGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCa            <  2:557015/140‑1 (MQ=17)

                                                                                                                           agGGGGACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGACTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGa          <  2:303292/140‑1 (MQ=17)

                                                                                                                             ggggACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACt        >  1:1458177/1‑140 (MQ=18)

                                                                                                                             ggggACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGAATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACt        >  2:1474732/1‑140 (MQ=12)

                                                                                                                              gggACCTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTc       >  2:10376/1‑140 (MQ=17)

                                                                                                                                   cTTCTGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACg  <  2:1458177/140‑1 (MQ=17)

GAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACG  >  NC_000913/3941883‑3942153

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`