BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,079,672	G→A	57.8%	R123H (CGC→CAC)	putP →	proline:sodium symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,079,672	0	G	A	57.8%	72.0 / 176.1	149	R123H (CGC→CAC)	putP	proline:sodium symporter
Reads supporting (aligned to +/- strand): ref base G (29/34); new base A (43/43); total (72/77)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

GGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTGAAAGTAC  >  NC_000913/1079573‑1079768

ggAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCg                                                                                                  <  1:2272565/100‑1 (MQ=255)

  aaGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCa                                                                                                <  2:3558882/100‑1 (MQ=255)

  aaGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCa                                                                                                <  2:1329739/100‑1 (MQ=255)

    gCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACAtt                                                                                              <  2:3332931/100‑1 (MQ=255)

      tggtggCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACAtttt                                                                                            >  2:2695365/1‑100 (MQ=255)

       ngtggCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTg                                                                                           >  1:1302275/2‑100 (MQ=255)

       ggtggCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTg                                                                                           <  2:1679300/100‑1 (MQ=255)

        gtggCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTgc                                                                                          <  1:1778490/100‑1 (MQ=255)

         tggCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTgcg                                                                                         <  2:2856895/100‑1 (MQ=255)

         tggCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTgcg                                                                                         <  2:2275421/100‑1 (MQ=255)

         tggCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTgcg                                                                                         <  2:243228/100‑1 (MQ=255)

         tggCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTgcg                                                                                         <  2:2149449/100‑1 (MQ=255)

          ggCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTgcgc                                                                                        <  2:2919775/100‑1 (MQ=255)

          ggCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTgcg                                                                                         <  1:3653181/99‑1 (MQ=255)

            ccGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCAt                                                                                      >  2:2773282/1‑100 (MQ=255)

            ccGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCAt                                                                                      <  2:2004151/100‑1 (MQ=255)

            ccGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTCGAAGATAAAAGCCGCATTTTgcgc                                                                                        >  1:3964527/1‑98 (MQ=255)

             cGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCAtt                                                                                     <  2:1574582/100‑1 (MQ=255)

              gggCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTa                                                                                    <  2:866619/100‑1 (MQ=255)

               ggCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTAt                                                                                   >  1:1107079/1‑100 (MQ=255)

                 cGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTAtct                                                                                 >  2:2118645/1‑100 (MQ=255)

                  ggTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTAtctc                                                                                >  2:394821/1‑100 (MQ=255)

                  ggTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTAtctc                                                                                >  1:703183/1‑100 (MQ=255)

                   gTTGCGTGTGCATACCGAATACAACAATAACGCCTTTACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTAtctc                                                                                <  2:447785/99‑1 (MQ=255)

                        gtgtGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCt                                                                          <  1:4106493/100‑1 (MQ=255)

                           tgCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGTGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGt                                                                       <  1:2365221/100‑1 (MQ=255)

                           tgCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGt                                                                       >  2:2347352/1‑100 (MQ=255)

                           tgCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGt                                                                       <  2:4223076/100‑1 (MQ=255)

                             cATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGt                                                                       >  1:1063047/1‑98 (MQ=255)

                               tACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTAtt                                                                   >  2:2007701/1‑100 (MQ=255)

                                aCCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAGAGCCGGATTTTGCGAATTATCTCTGCGCTGGTTAttt                                                                  >  1:2881953/1‑100 (MQ=255)

                                aCCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTAttt                                                                  <  2:101808/100‑1 (MQ=255)

                                 ccGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTAtttt                                                                 <  2:2991832/100‑1 (MQ=255)

                                  cGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTg                                                                <  2:17019/100‑1 (MQ=255)

                                  cGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTg                                                                <  1:1287633/100‑1 (MQ=255)

                                  cGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTAtttt                                                                 >  1:4219356/1‑99 (MQ=255)

                                  cGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTg                                                                >  2:1254897/1‑100 (MQ=255)

                                   gAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGc                                                               <  2:1147746/100‑1 (MQ=255)

                                   gAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGc                                                               <  1:4130977/100‑1 (MQ=255)

                                    aaTACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCt                                                              <  2:2359789/100‑1 (MQ=255)

                                     aTACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTg                                                             <  1:2157587/100‑1 (MQ=255)

                                     aTACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTg                                                             >  2:2752394/1‑100 (MQ=255)

                                        caacaaTAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCCTTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGttc                                                          >  1:810016/1‑100 (MQ=255)

                                        caacaaTAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGttc                                                          >  1:2847992/1‑100 (MQ=255)

                                        caacaaTAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGttc                                                          >  1:1345393/1‑100 (MQ=255)

                                        caacaaTAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGttc                                                          >  1:3081189/1‑100 (MQ=255)

                                         aacaaTAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGttct                                                         <  2:227114/100‑1 (MQ=255)

                                          acaaTAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGttctt                                                        >  2:33878/1‑100 (MQ=255)

                                              tAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCAcc                                                    >  1:2079759/1‑100 (MQ=255)

                                              tAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCAcc                                                    >  1:1991960/1‑100 (MQ=255)

                                               aaCGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCa                                                   <  2:2852665/100‑1 (MQ=255)

                                                aCGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCAt                                                  <  2:2916164/100‑1 (MQ=255)

                                                 cGCCTTAACACTGCCGGATTATTTCACCGTGCGCTTTGAAGATAAATGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTACTTCACCAtt                                                 <  1:3738628/100‑1 (MQ=255)

                                                 cGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCAtt                                                 >  2:4221792/1‑100 (MQ=255)

                                                   ccTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTa                                               >  1:1396553/1‑100 (MQ=255)

                                                     ttAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTAtt                                             >  1:3826365/1‑100 (MQ=255)

                                                      tAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTg                                            >  1:4018624/1‑100 (MQ=255)

                                                       aaCACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTgc                                           <  2:270629/100‑1 (MQ=255)

                                                       aaCACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTgc                                           <  2:4023977/100‑1 (MQ=255)

                                                       aaCACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTgc                                           <  2:1645208/100‑1 (MQ=255)

                                                       aaCACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTgc                                           <  2:1231823/100‑1 (MQ=255)

                                                       aaCACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTgc                                           <  2:467112/100‑1 (MQ=255)

                                                          acTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCt                                        <  1:1137965/100‑1 (MQ=255)

                                                          acTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCt                                        <  2:3738195/100‑1 (MQ=255)

                                                           cTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCtt                                       >  2:1771723/1‑100 (MQ=255)

                                                             gCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCg                                     <  2:1719882/100‑1 (MQ=255)

                                                               cGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCggg                                   <  2:1087171/100‑1 (MQ=255)

                                                                     attTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATtgt                             >  2:1276884/1‑100 (MQ=255)

                                                                      tttCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATtgtg                            >  2:3154349/1‑100 (MQ=255)

                                                                      tttCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGGGCTTCGGGCATtgt                             >  1:1011549/1‑99 (MQ=255)

                                                                      tttCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATtgtg                            <  1:3491015/100‑1 (MQ=255)

                                                                      tgtCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATtgt                             <  1:2601173/97‑1 (MQ=255)

                                                                       ttCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTgg                           <  2:1074752/100‑1 (MQ=255)

                                                                       ttCACCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTgg                           <  2:3597817/100‑1 (MQ=255)

                                                                         cACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCa                         >  1:3647111/1‑100 (MQ=255)

                                                                          aCCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAg                        >  2:3165252/1‑100 (MQ=255)

                                                                          aCCGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTACTTTACCAATTATTGCGCGTCGGGCAtggtgg                           >  1:2100965/1‑97 (MQ=255)

                                                                           ccGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAgg                       >  1:103353/1‑100 (MQ=255)

                                                                           ccGGGCGCTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAgg                       >  2:3604049/1‑100 (MQ=255)

                                                                             gggCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGgcg                     <  1:1741837/100‑1 (MQ=255)

                                                                               gcgcTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGgcgcg                   <  2:3764478/100‑1 (MQ=255)

                                                                                  cTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGt                >  2:2327740/1‑100 (MQ=255)

                                                                                  cTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGt                >  1:163695/1‑100 (MQ=255)

                                                                                  cTTTGAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGt                >  1:2758650/1‑100 (MQ=255)

                                                                                      gAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGt            >  1:2702082/1‑100 (MQ=255)

                                                                                      gAAGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGt            <  1:4192646/100‑1 (MQ=255)

                                                                                       aaGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGtt           >  1:1674543/1‑100 (MQ=255)

                                                                                        aGATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTGATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGttt          <  1:941599/100‑1 (MQ=255)

                                                                                         gATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTg         <  2:885385/100‑1 (MQ=255)

                                                                                         gATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTg         >  2:355758/1‑100 (MQ=255)

                                                                                         gATAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTg         <  2:1041710/100‑1 (MQ=255)

                                                                                           tAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTGaa       >  1:2028146/1‑100 (MQ=255)

                                                                                           tAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTGaa       >  2:1196145/1‑100 (MQ=255)

                                                                                           tAAAAGCCACATTTTGCGCATTATCTCTGCGCTGGTCATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTGaa       >  1:418948/1‑100 (MQ=255)

                                                                                             aaaGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCTGGCATTGTGGCAGGCGCGCGTCTGTTTGAAAg     <  2:4051591/100‑1 (MQ=255)

                                                                                             aaaGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTGAAAg     <  1:303682/100‑1 (MQ=255)

                                                                                             aaaGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTGAAAg     <  1:3461854/100‑1 (MQ=255)

                                                                                             aaaGCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTGAAAg     <  2:3570617/100‑1 (MQ=255)

                                                                                               aGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTGAAAGTa   >  1:1146575/1‑100 (MQ=255)

                                                                                                gCCACATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTGAAAGTAc  >  2:3037846/1‑100 (MQ=255)

GGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACAATAACGCCTTAACACTGCCGGATTATTTCACCGGGCGCTTTGAAGATAAAAGCCGCATTTTGCGCATTATCTCTGCGCTGGTTATTTTGCTGTTCTTCACCATTTATTGCGCTTCGGGCATTGTGGCAGGCGCGCGTCTGTTTGAAAGTAC  >  NC_000913/1079573‑1079768

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