BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,561,617	0	A	G	34.4%	526.5 / inf	392	S353S (AGT→AGC)	yjiJ	putative transporter YjiJ
Reads supporting (aligned to +/- strand): ref base A (115/142); new base G (58/77); total (173/219)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-03
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 396 total aligned reads displayed.

GCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGTTGGCGGGCGATAG  >  NC_000913/4561517‑4561717

 cTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATAc                                                                                                     <  1:8866548/101‑1 (MQ=255)

   tAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTg                                                                                                   >  2:8867814/1‑101 (MQ=255)

    aaaTAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTgg                                                                                                  <  2:8867463/101‑1 (MQ=255)

    aaaTAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTgg                                                                                                  <  2:8867578/101‑1 (MQ=255)

     aataatTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGt                                                                                                 <  1:8867168/101‑1 (MQ=255)

      ataatTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTc                                                                                                <  2:8867432/101‑1 (MQ=255)

       taatTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCa                                                                                               >  1:8867820/1‑101 (MQ=255)

       taatTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCa                                                                                               >  2:8867819/1‑101 (MQ=255)

        aatTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAg                                                                                              <  1:8867441/101‑1 (MQ=255)

        aatTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAg                                                                                              <  2:8867501/101‑1 (MQ=255)

         atTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAgc                                                                                             <  2:8867309/101‑1 (MQ=255)

         atTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAgc                                                                                             <  1:8867279/101‑1 (MQ=255)

         atTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAgc                                                                                             >  1:8867825/1‑101 (MQ=255)

           tGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAgcgc                                                                                           <  1:8866935/101‑1 (MQ=255)

            gCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGcc                                                                                          >  2:8867830/1‑101 (MQ=255)

              aCTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCgg                                                                                        >  1:8867833/1‑101 (MQ=255)

                tgtgCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGc                                                                                      <  1:8867293/101‑1 (MQ=255)

                tgtgCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGc                                                                                      <  1:8867261/101‑1 (MQ=255)

                   gCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGa                                                                                   >  1:8867841/1‑101 (MQ=255)

                   gCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGa                                                                                   >  2:8867840/1‑101 (MQ=255)

                     tGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAgg                                                                                 <  2:8867237/101‑1 (MQ=255)

                       aTAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGAt                                                                               <  1:8867349/101‑1 (MQ=255)

                       aTAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGAt                                                                               >  1:8867846/1‑101 (MQ=255)

                        tAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGAtt                                                                              <  2:8867417/101‑1 (MQ=255)

                         aaGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGAttt                                                                             <  2:8867050/101‑1 (MQ=255)

                         aaGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGAttt                                                                             >  2:8867852/1‑101 (MQ=255)

                          aGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTg                                                                            <  1:8867230/101‑1 (MQ=255)

                          aGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTg                                                                            >  1:8867853/1‑101 (MQ=255)

                              cagcaATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCa                                                                        <  1:8867500/101‑1 (MQ=255)

                              cagcaATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCa                                                                        <  2:8867031/101‑1 (MQ=255)

                              cagcaATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCa                                                                        >  1:8867858/1‑101 (MQ=255)

                               agcaATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCaa                                                                       <  1:8867699/101‑1 (MQ=255)

                                  aaTAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCaa                                                                       >  2:8867864/1‑98 (MQ=255)

                                  aaTAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGc                                                                    <  1:8867381/101‑1 (MQ=255)

                                     aaaCAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGt                                                                 >  2:8867868/1‑101 (MQ=255)

                                     aaaCAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGcc                                                                   >  2:8867872/1‑99 (MQ=255)

                                      aacaacGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTa                                                                <  2:8867199/101‑1 (MQ=255)

                                      aacaacGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTa                                                                >  1:8867876/1‑101 (MQ=255)

                                       acaacGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTaa                                                               >  2:8867878/1‑101 (MQ=255)

                                         aacGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAAt                                                             <  2:8867296/101‑1 (MQ=255)

                                         aacGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAAt                                                             <  2:8867170/101‑1 (MQ=255)

                                          acGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGAttt                                                                             >  1:8867886/1‑84 (MQ=255)

                                          acGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATg                                                            >  2:8867885/1‑101 (MQ=255)

                                            gccgccgcGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAg                                                          <  2:8867403/101‑1 (MQ=255)

                                             ccgccgcGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGc                                                         <  2:8867796/101‑1 (MQ=255)

                                             ccgccgcGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGc                                                         <  2:8866978/101‑1 (MQ=255)

                                             ccgccgcGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGc                                                         <  1:8867138/101‑1 (MQ=255)

                                              cgccgcgccgCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCg                                                        >  1:8867892/1‑101 (MQ=255)

                                               gccgcgccgcAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGt                                                       >  2:8867897/1‑101 (MQ=255)

                                               gccgcgccgcAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGt                                                       >  1:8867896/1‑101 (MQ=255)

                                                 cgcgccgcAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTca                                                     >  1:8867902/1‑101 (MQ=255)

                                                   cgccgcAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTcaca                                                   >  2:8867905/1‑101 (MQ=255)

                                                     ccgcAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACaaa                                                 >  2:8867911/1‑101 (MQ=255)

                                                     ccgcAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACaaa                                                 <  1:8867336/101‑1 (MQ=255)

                                                     ccgcAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACaaa                                                 <  1:8867045/101‑1 (MQ=255)

                                                      cgcAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAg                                                <  2:8867485/101‑1 (MQ=255)

                                                       gcAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGc                                               >  1:8867919/1‑101 (MQ=255)

                                                       gcAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGc                                               >  2:8867916/1‑101 (MQ=255)

                                                         agagagTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCt                                             >  1:8867926/1‑101 (MQ=255)

                                                         agagagTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCt                                             >  2:8867922/1‑101 (MQ=255)

                                                          gagagTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTa                                            <  1:8867549/101‑1 (MQ=255)

                                                           agagTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTaa                                           <  1:8867477/101‑1 (MQ=255)

                                                            gagTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTaaa                                          <  1:8867118/101‑1 (MQ=255)

                                                            gagTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTaaa                                          <  1:8867635/101‑1 (MQ=255)

                                                            gagTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTaaa                                          <  2:8867268/101‑1 (MQ=255)

                                                            gagTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTaaa                                          <  1:8867523/101‑1 (MQ=255)

                                                             agTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTaaaa                                         <  1:8867357/101‑1 (MQ=255)

                                                              gTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAg                                        <  1:8867201/101‑1 (MQ=255)

                                                               tGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGa                                       <  1:8867425/101‑1 (MQ=255)

                                                                ggCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCt                                                                                                    <  2:8867142/39‑1 (MQ=255)

                                                                ggCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGAt                                      <  1:8867342/101‑1 (MQ=255)

                                                                ggCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGAt                                      <  2:8867181/101‑1 (MQ=255)

                                                                 gcgcTGGCAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGAtt                                     >  1:8867933/1‑101 (MQ=255)

                                                                 gcgcTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGAtt                                     >  2:8867932/1‑101 (MQ=255)

                                                                 gcgcTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGAtt                                     <  2:8867517/101‑1 (MQ=255)

                                                                     tGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTgaga                                 >  2:8867939/1‑101 (MQ=255)

                                                                        cAAGCGCCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGAttt                              >  1:8867944/1‑101 (MQ=255)

                                                                        cAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGccc                                                                                     >  2:8867945/1‑46 (MQ=255)

                                                                         aaGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTc                             <  2:8867263/101‑1 (MQ=255)

                                                                           gcgcCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTccc                           >  2:8867951/1‑101 (MQ=255)

                                                                            cgcCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCg                          <  2:8867561/101‑1 (MQ=255)

                                                                            cgcCGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCg                          <  1:8867318/101‑1 (MQ=255)

                                                                             gcCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCgg                         >  2:8867956/1‑101 (MQ=255)

                                                                              ccGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGt                        >  1:8867958/1‑101 (MQ=255)

                                                                              ccGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGGCACAAAGCCTAAAAGATTGAGATTTCCCGGt                        <  2:8867346/101‑1 (MQ=255)

                                                                               cGATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTa                       >  2:8867963/1‑101 (MQ=255)

                                                                                gATGTTCCATTGCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTAc                      >  1:8867966/1‑101 (MQ=255)

                                                                                gATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTAc                      <  2:8867453/101‑1 (MQ=255)

                                                                                    ttCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTc                  <  1:8867451/101‑1 (MQ=255)

                                                                                    ttCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTc                  >  1:8867972/1‑101 (MQ=255)

                                                                                     tCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCa                 >  2:8867973/1‑101 (MQ=255)

                                                                                       cATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGt               <  1:8867599/101‑1 (MQ=255)

                                                                                         ttGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGTTg             >  2:8867979/1‑101 (MQ=255)

                                                                                           gCTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGTTggc           <  1:8867366/101‑1 (MQ=255)

                                                                                           gCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGTTggc           <  1:8867219/101‑1 (MQ=255)

                                                                                           gCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGTTggc           >  1:8867981/1‑101 (MQ=255)

                                                                                            cTGAGCATGCTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGTTggcg          <  2:8867444/101‑1 (MQ=255)

                                                                                             tGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGTTggcgg         <  2:8867647/101‑1 (MQ=255)

                                                                                             tGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGTTggcgg         <  2:8867539/101‑1 (MQ=255)

                                                                                                  aTACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGTTGGCGGGCGAt    <  1:8867394/101‑1 (MQ=255)

GCTTAAATAATTGCACTGTGCTGATAAGCGCAGCAATAAACAACGCCGCCGCGCCGCAGAGAGTGGCGCTGGCAAGCGCCGATGTTCCATTGCTGAGCATACTGGTCAGCGCCGGGCCGAGGATTTGCCCAATGCCGTAAATGAGCGTCACAAAGCCTAAAAGATTGAGATTTCCCGGTACGCTCAGTTGGCGGGCGATAG  >  NC_000913/4561517‑4561717

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: `‑`