BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	158,711	C→T	R128H (CGC→CAC)	pcnB ←	poly(A) polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	158,711	0	C	T	99.7%	1166.0 / ‑6.0	325	R128H (CGC→CAC)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (164/160); total (164/161)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 339 total aligned reads displayed.

TACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGA  >  NC_000913/158621‑158800

tACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCgtgt                                                                                         >  2:3445726/1‑93 (MQ=255)

tACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCgtgt                                                                                         >  1:3886309/1‑93 (MQ=255)

 aCAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGtt                                                                                        >  1:4893470/1‑93 (MQ=255)

   aGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTgg                                                                                      <  2:3527345/93‑1 (MQ=255)

     gCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGa                                                                                    <  2:3896205/93‑1 (MQ=255)

     gCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGa                                                                                    <  2:8485300/93‑1 (MQ=255)

      gtGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAg                                                                                   <  1:1288749/92‑1 (MQ=255)

      cTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAg                                                                                   <  1:5434564/93‑1 (MQ=255)

        gTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGt                                                                                 >  1:464066/1‑93 (MQ=255)

        gTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGt                                                                                 >  1:5604307/1‑93 (MQ=255)

          tGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCg                                                                               <  1:3147352/93‑1 (MQ=255)

            aTAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTg                                                                             >  1:9004958/1‑93 (MQ=255)

              aGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCg                                                                           <  1:5044461/93‑1 (MQ=255)

               gTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCgg                                                                          >  2:1475106/1‑93 (MQ=255)

                tGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGt                                                                         >  2:9170122/1‑93 (MQ=255)

                tGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGt                                                                         <  1:9432213/93‑1 (MQ=255)

                 gAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTc                                                                        <  2:5459134/93‑1 (MQ=255)

                  aaaTCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCg                                                                       >  1:8280807/1‑93 (MQ=255)

                  aaaTCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCg                                                                       >  2:1654210/1‑93 (MQ=255)

                  aaaTCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCg                                                                       >  1:910972/1‑93 (MQ=255)

                   aaTCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGc                                                                      >  1:2262356/1‑93 (MQ=255)

                       gcggcgCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAGCATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGAc                                                                  <  2:4995467/93‑1 (MQ=255)

                         ggcgCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGt                                                                <  1:9255792/93‑1 (MQ=255)

                          gcgcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGtt                                                               >  1:751710/1‑93 (MQ=255)

                          gcgcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGtt                                                               >  1:7931615/1‑93 (MQ=255)

                          gcgcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGtt                                                               >  2:7592440/1‑93 (MQ=255)

                          gcgcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGtt                                                               <  1:4291972/93‑1 (MQ=255)

                           cgcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCTGTTTTGCCGGTGTTGGGGGGGGGGGCGGTGGCTGACGTTa                                                              >  2:7178953/1‑93 (MQ=255)

                            gcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTAc                                                             >  1:605855/1‑93 (MQ=255)

                             cTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTAcc                                                            <  2:7003688/93‑1 (MQ=255)

                              tGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCt                                                           <  2:8799222/93‑1 (MQ=255)

                               gggCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCtt                                                          <  1:9849454/93‑1 (MQ=255)

                                ggCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTc                                                         >  2:2020927/1‑93 (MQ=255)

                                 gCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCg                                                        >  1:2892617/1‑93 (MQ=255)

                                    tcttctTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTCCCTTCgtgg                                                     >  1:2344844/1‑93 (MQ=255)

                                    tcttctTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtgg                                                     >  2:8177283/1‑93 (MQ=255)

                                    tcttctTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtgg                                                     >  1:6903393/1‑93 (MQ=255)

                                    tcttctTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtgg                                                     <  1:9892461/93‑1 (MQ=255)

                                      ttcttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtggtg                                                   <  1:7228261/93‑1 (MQ=255)

                                      ttcttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtgctg                                                   >  1:3444542/1‑93 (MQ=255)

                                       tcttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGgtgt                                                  >  1:6019410/1‑93 (MQ=255)

                                       tcttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGgtgt                                                  <  1:1959593/93‑1 (MQ=255)

                                        cttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTc                                                 >  1:184781/1‑93 (MQ=255)

                                         ttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGAGGTTACCTTCGTGGTGTcc                                                >  1:4123519/1‑93 (MQ=255)

                                          tcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCa                                               <  1:2685366/93‑1 (MQ=255)

                                            gATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACg                                             >  1:819096/1‑93 (MQ=255)

                                            gATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACg                                             <  2:1850422/93‑1 (MQ=255)

                                             aTGGAGCCAAAAATGTTGTCGCGCAGCAACATGTGGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACgg                                            <  1:9701467/93‑1 (MQ=255)

                                              tGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGa                                           >  2:3058045/1‑93 (MQ=255)

                                               ggAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGaa                                          <  2:4130774/93‑1 (MQ=255)

                                               ggAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGaa                                          >  1:5582979/1‑93 (MQ=255)

                                               ggAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGaa                                          <  2:117209/93‑1 (MQ=255)

                                                gAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAg                                         <  1:10007085/93‑1 (MQ=255)

                                                gAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAg                                         >  1:1637140/1‑93 (MQ=255)

                                                 aGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAgg                                        <  1:9960843/93‑1 (MQ=255)

                                                    cGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCg                                     >  2:4087592/1‑93 (MQ=255)

                                                     gAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGc                                    >  2:1709402/1‑93 (MQ=255)

                                                     gAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGc                                    >  1:8390613/1‑93 (MQ=255)

                                                        aaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAAc                                 <  2:4730079/93‑1 (MQ=255)

                                                        aaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAAc                                 >  2:3092118/1‑93 (MQ=255)

                                                        aaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCATGGTGTCCACGGAAGGTCGCAAc                                 >  1:7411446/1‑93 (MQ=255)

                                                         aTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACt                                <  2:2157676/93‑1 (MQ=255)

                                                          tgttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACtt                               <  1:3468974/93‑1 (MQ=255)

                                                          tgttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACtt                               <  2:9303082/93‑1 (MQ=255)

                                                          tgttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACtt                               <  2:5684637/93‑1 (MQ=255)

                                                            ttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCg                             >  1:1211554/1‑93 (MQ=255)

                                                             tgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGa                            >  2:8385138/1‑93 (MQ=255)

                                                             tgtAGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGa                            <  1:6407902/93‑1 (MQ=255)

                                                              gtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGAt                           >  1:2996101/1‑93 (MQ=255)

                                                              gtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGAt                           <  2:6043909/93‑1 (MQ=255)

                                                               tCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATa                          >  2:2391867/1‑93 (MQ=255)

                                                               tCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATa                          <  2:8238937/93‑1 (MQ=255)

                                                               tCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATa                          <  2:6405667/93‑1 (MQ=255)

                                                               tCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATa                          <  1:7739800/93‑1 (MQ=255)

                                                                 gcgcAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAt                        <  1:1662325/93‑1 (MQ=255)

                                                                 gcgcAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAt                        >  1:9164465/1‑93 (MQ=255)

                                                                 gcgcAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAt                        <  2:5042002/93‑1 (MQ=255)

                                                                 gcgcAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGACAAt                        <  2:9839146/93‑1 (MQ=255)

                                                                   gcagcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAtct                      <  1:8641973/93‑1 (MQ=255)

                                                                      gcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCg                   <  2:9713477/93‑1 (MQ=255)

                                                                      gcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCg                   >  1:5501084/1‑93 (MQ=255)

                                                                       caACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCgg                  <  2:6745027/93‑1 (MQ=255)

                                                                        aaCATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCggg                 >  2:5941709/1‑93 (MQ=255)

                                                                        aaCATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCggg                 >  2:2888298/1‑93 (MQ=255)

                                                                         aCATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGc                >  2:6244101/1‑93 (MQ=255)

                                                                          cATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGcc               <  1:4474594/93‑1 (MQ=255)

                                                                           aTGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGGAGGTCGCAACTTCGATAATCTCGGGGCCa              >  1:3702414/1‑93 (MQ=255)

                                                                            tGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCaa             >  2:4348477/1‑93 (MQ=255)

                                                                            tGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCaa             >  2:2319245/1‑93 (MQ=255)

                                                                            tGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCaa             >  1:2147851/1‑93 (MQ=255)

                                                                            tGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCaa             >  2:2618611/1‑93 (MQ=255)

                                                                             gCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCaaa            >  2:796720/1‑93 (MQ=255)

                                                                               cGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACa          <  2:7680432/93‑1 (MQ=255)

                                                                                 ttttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACAtt        >  2:4229867/1‑93 (MQ=255)

                                                                                 ttttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACAtt        <  2:898643/93‑1 (MQ=255)

                                                                                     gCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACAt    >  2:6898873/1‑93 (MQ=255)

                                                                                      cccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATg   <  2:7767500/93‑1 (MQ=255)

                                                                                      cccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATg   >  1:7899920/1‑93 (MQ=255)

                                                                                       ccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGa  <  1:5216102/93‑1 (MQ=255)

                                                                                       ccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGa  >  2:8033027/1‑93 (MQ=255)

TACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGA  >  NC_000913/158621‑158800

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