BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,846,935	G→A	P9L (CCG→CTG)	topB ←	DNA topoisomerase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,846,935	0	G	A	99.6%	1148.8 / ‑4.4	267	P9L (CCG→CTG)	topB	DNA topoisomerase III
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (110/156); total (111/156)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 267 total aligned reads displayed.

TACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAA  >  NC_000913/1846837‑1847032

tACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGt                                                                                                 <  2:3326009/101‑1 (MQ=255)

tACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGt                                                                                                 >  1:3326437/1‑101 (MQ=255)

 acacCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGtt                                                                                                <  1:3326126/101‑1 (MQ=255)

   acCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGtttt                                                                                              <  1:3325836/101‑1 (MQ=255)

   acCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGtttt                                                                                              >  1:3326442/1‑101 (MQ=255)

    ccAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGttttt                                                                                             >  2:3326445/1‑101 (MQ=255)

      aGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCg                                                                                           <  2:3325925/101‑1 (MQ=255)

          caccacCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCaa                                                                                       <  2:3325961/101‑1 (MQ=255)

          caccacCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCaa                                                                                       >  2:3326449/1‑101 (MQ=255)

           accaccTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAAt                                                                                      <  1:3325891/101‑1 (MQ=255)

           accaccTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAAt                                                                                      >  1:3326450/1‑101 (MQ=255)

              accTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATaaa                                                                                   <  1:3325941/101‑1 (MQ=255)

               ccTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAaac                                                                                  <  2:3326229/101‑1 (MQ=255)

                cTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATaa                                                                                    <  2:3325856/98‑1 (MQ=255)

                 tGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAaacaa                                                                                >  1:3326454/1‑101 (MQ=255)

                  gACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAaacaac                                                                               <  2:3325831/101‑1 (MQ=255)

                  gACCATTACCGCACTCGATAAAGCCATCGCCTTCCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAaacaac                                                                               <  2:3325766/101‑1 (MQ=255)

                   aCCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAAcc                                                                              <  2:3325902/101‑1 (MQ=255)

                    ccATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCg                                                                             >  2:3326458/1‑101 (MQ=255)

                     cATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGc                                                                            >  2:3326461/1‑101 (MQ=255)

                      aTTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCa                                                                           >  1:3326463/1‑101 (MQ=255)

                       ttACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCAt                                                                          >  2:3326465/1‑101 (MQ=255)

                        tACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCAtt                                                                         <  1:3325895/101‑1 (MQ=255)

                          ccGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTg                                                                        <  2:3326248/100‑1 (MQ=255)

                          ccGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGa                                                                       <  2:3326074/101‑1 (MQ=255)

                           cGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGaa                                                                      <  2:3325995/101‑1 (MQ=255)

                            gCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAAt                                                                     <  1:3326004/101‑1 (MQ=255)

                             cACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAAtt                                                                    >  2:3326470/1‑101 (MQ=255)

                             cACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAAtt                                                                    <  2:3325870/101‑1 (MQ=255)

                             cACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAAtt                                                                    >  1:3326469/1‑101 (MQ=255)

                              aCTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTa                                                                   <  1:3326103/101‑1 (MQ=255)

                              aCTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTa                                                                   <  1:3326082/101‑1 (MQ=255)

                               cTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTaa                                                                  <  2:3325906/101‑1 (MQ=255)

                                tCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAAc                                                                 <  2:3325731/101‑1 (MQ=255)

                                 cGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTaa                                                                  <  1:3326041/99‑1 (MQ=255)

                                 cGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACg                                                                <  1:3325980/101‑1 (MQ=255)

                                  gATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGa                                                               <  1:3326023/101‑1 (MQ=255)

                                  gATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGa                                                               >  2:3326474/1‑101 (MQ=255)

                                  gATAAAGCCATCGCCTTTCCGGAGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGa                                                               <  1:3325679/101‑1 (MQ=255)

                                   aTAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGaa                                                              <  2:3325942/101‑1 (MQ=255)

                                   aTAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGaa                                                              <  1:3325998/101‑1 (MQ=255)

                                   aTAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGaa                                                              >  1:3326477/1‑101 (MQ=255)

                                   aTAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGaa                                                              <  1:3326054/101‑1 (MQ=255)

                                   aTAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCACTTTTTCGGCAATAAACAACCGCATTGAATTAACGaa                                                              <  2:3326211/101‑1 (MQ=255)

                                    tAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAAt                                                             <  1:3325956/101‑1 (MQ=255)

                                    tAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAAt                                                             <  2:3326130/101‑1 (MQ=255)

                                     aaaGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAAtc                                                            <  2:3326016/101‑1 (MQ=255)

                                      aaGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAAtct                                                           >  2:3326480/1‑101 (MQ=255)

                                       aGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAAtctc                                                          <  1:3326032/101‑1 (MQ=255)

                                        gCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCa                                                         >  1:3326481/1‑101 (MQ=255)

                                        gCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCa                                                         <  2:3326061/101‑1 (MQ=255)

                                         ccATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCaa                                                        <  2:3326044/101‑1 (MQ=255)

                                          cATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAAc                                                       <  1:3325971/101‑1 (MQ=255)

                                           aTCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAAcc                                                      >  1:3326483/1‑101 (MQ=255)

                                             cGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCAt                                                    <  1:3325983/101‑1 (MQ=255)

                                              gCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATg                                                   >  1:3326486/1‑101 (MQ=255)

                                              gCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATg                                                   >  2:3326489/1‑101 (MQ=255)

                                              gCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATg                                                   >  2:3326491/1‑101 (MQ=255)

                                                 tttCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCa                                                <  1:3326111/101‑1 (MQ=255)

                                                 tttCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCa                                                >  2:3326494/1‑101 (MQ=255)

                                                 tttCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCa                                                <  1:3325931/101‑1 (MQ=255)

                                                  ttCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCAc                                               <  2:3326350/101‑1 (MQ=255)

                                                   tCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACg                                              >  2:3326498/1‑101 (MQ=255)

                                                     cGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGa                                             <  1:3326172/100‑1 (MQ=255)

                                                     cGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGAc                                            >  2:3326500/1‑101 (MQ=255)

                                                      ggTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGAcc                                           <  2:3325967/101‑1 (MQ=255)

                                                        tGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGAccgc                                         >  1:3326503/1‑101 (MQ=255)

                                                         gCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGAccgcc                                        >  1:3326505/1‑101 (MQ=255)

                                                           ggTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACcgccgc                                      >  2:3326510/1‑101 (MQ=255)

                                                           ggTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACcgccgc                                      >  1:3326509/1‑101 (MQ=255)

                                                                gggCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACcgccgc                                      <  1:3325945/96‑1 (MQ=255)

                                                                 ggCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTg                                <  2:3325976/101‑1 (MQ=255)

                                                                  gCAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTgg                               >  2:3326514/1‑101 (MQ=255)

                                                                   cAGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGc                              >  1:3326517/1‑101 (MQ=255)

                                                                    aGGACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCa                             >  2:3326519/1‑101 (MQ=255)

                                                                      gACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCAcc                           >  2:3326522/1‑101 (MQ=255)

                                                                      gACATCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCAcc                           >  1:3326523/1‑101 (MQ=255)

                                                                         aTCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCAc                                               >  2:3326530/1‑78 (MQ=255)

                                                                         aTCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGt                        <  2:3326100/101‑1 (MQ=255)

                                                                         aTCAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGt                        >  2:3326527/1‑101 (MQ=255)

                                                                           cAGCAATGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTc                      <  1:3325604/101‑1 (MQ=255)

                                                                               aaTGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCa                  <  1:3326266/101‑1 (MQ=255)

                                                                               aaTGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCa                  >  1:3326532/1‑101 (MQ=255)

                                                                                aTGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCaa                 <  2:3326089/101‑1 (MQ=255)

                                                                                aTGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCaa                 <  2:3326030/101‑1 (MQ=255)

                                                                                 tGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAAt                <  2:3326186/101‑1 (MQ=255)

                                                                                 tGGCGCGCGCCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAAt                <  1:3326446/101‑1 (MQ=255)

                                                                                     gcgcgcCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGcc            >  2:3326536/1‑101 (MQ=255)

                                                                                      cgcgcCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCg           <  1:3326189/101‑1 (MQ=255)

                                                                                        cgcCATACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACg                                              <  2:3326067/64‑1 (MQ=255)

                                                                                        cgcCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTc         <  2:3326164/101‑1 (MQ=255)

                                                                                        cgcCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTc         >  1:3326537/1‑101 (MQ=255)

                                                                                         gcCAGCCTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCa        <  1:3326122/101‑1 (MQ=255)

                                                                                         gcCAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCg           >  2:3326540/1‑98 (MQ=255)

                                                                                           cAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGt      <  2:3326157/101‑1 (MQ=255)

                                                                                           cAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGt      <  1:3326150/101‑1 (MQ=255)

                                                                                           cAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGt      <  1:3326331/101‑1 (MQ=255)

                                                                                           cAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGt      >  1:3326545/1‑101 (MQ=255)

                                                                                           cAGACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGt      >  1:3326542/1‑101 (MQ=255)

                                                                                             gACTCAGTTTTTCGGCAATAAACAACCGCATTGAATTaa                                                                  <  2:3325681/39‑1 (MQ=255)

TACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAA  >  NC_000913/1846837‑1847032

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