BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,796,375	T→A	15.8%	I138F (ATT→TTT)	ygaU ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,796,375	0	T	A	15.8%	423.9 / 63.5	165	I138F (ATT→TTT)	ygaU	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base T (67/72); new base A (4/22); total (71/94)
Fisher's exact test for biased strand distribution p-value = 2.12e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 170 total aligned reads displayed.

AATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGT  >  NC_000913/2796276‑2796473

aaTGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAt                                                                                                    <  2:3857961/100‑1 (MQ=255)

  tGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAttt                                                                                                  >  2:534319/1‑100 (MQ=255)

  tGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAttt                                                                                                  <  1:3857410/100‑1 (MQ=255)

   gcgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCACGGATAAAtttt                                                                                                 >  2:1116610/1‑100 (MQ=255)

    cgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGCTGGCTGGCCATTAACTACTCTTCAGGAATACGCAACACCTGCCCCGGACAAAttt                                                                                                  >  2:2594738/1‑98 (MQ=255)

    cgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTa                                                                                                <  1:1912420/100‑1 (MQ=255)

     gcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTAt                                                                                               >  2:2881036/1‑100 (MQ=255)

      ccAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATc                                                                                              >  1:914918/1‑100 (MQ=255)

         gCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCgg                                                                                           <  1:3935381/100‑1 (MQ=255)

          cATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCggg                                                                                          <  1:1092540/100‑1 (MQ=255)

           aTTTACTGACGCGAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCgggg                                                                                         <  1:1316302/100‑2 (MQ=255)

           aTTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGc                                                                                         >  1:3154237/1‑100 (MQ=255)

           aTTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCgggg                                                                                         <  1:2028685/100‑2 (MQ=255)

           aTTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCgggg                                                                                         <  1:1470177/100‑2 (MQ=255)

           aTTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCgggg                                                                                         <  1:2783385/100‑2 (MQ=255)

           aTTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTCGCCCCGGATAAAATTTATCCgggg                                                                                         <  2:452929/100‑2 (MQ=255)

            tttACTGACTCCAATTCCACTTATCCAGCAGGCGATGGCTGGAAATTAACTACTTTTCCGGAATACGCAACACCTGCCCCGGATAAAATTTATCCgggg                                                                                         <  2:551055/99‑2 (MQ=255)

            tttACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCt                                                                                        <  2:3058042/100‑1 (MQ=255)

            tttACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCt                                                                                        <  1:2360800/100‑1 (MQ=255)

               aCTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTG‑CCCGGATAAAATTTATCCggggcaagt                                                                                    <  2:2635242/100‑7 (MQ=38)

               aCTGACGCAAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCggggcaag                                                                                     <  2:225694/100‑6 (MQ=39)

                cTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCttttt                                                                                    <  1:1515382/100‑1 (MQ=255)

                cTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCttttt                                                                                    >  2:1144605/1‑100 (MQ=255)

                  gACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAg                                                                                  <  1:2154384/100‑1 (MQ=255)

                  gACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAg                                                                                  <  1:3168103/100‑1 (MQ=255)

                  gACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAg                                                                                  <  1:3262629/100‑1 (MQ=255)

                    cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCa                                                                                >  1:3410048/1‑100 (MQ=255)

                       cAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCg                                                                             <  2:3154216/100‑1 (MQ=255)

                       cAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCg                                                                             >  2:1869670/1‑100 (MQ=255)

                        aaTTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCgg                                                                            <  1:3766032/100‑1 (MQ=255)

                        aaTTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCgg                                                                            <  1:1948399/100‑1 (MQ=255)

                        aaTTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCgg                                                                            <  2:2476663/100‑1 (MQ=255)

                           tcacATTATCCAGCAGCCGATGGTTGGCACTTAATTACTCTTCATGAATACGCAATACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGttt                                                                         <  2:3578989/97‑1 (MQ=255)

                           tCCCATTATCCAGCAGGCGAGGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGttt                                                                         <  2:3680906/100‑1 (MQ=255)

                            cccATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTa                                                                        >  2:174300/1‑100 (MQ=255)

                            cccATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTa                                                                        >  1:1068914/1‑100 (MQ=255)

                             ccATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTAt                                                                       >  1:3685901/1‑100 (MQ=255)

                              cATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTAtt                                                                      <  2:1718095/100‑1 (MQ=255)

                               aTTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTc                                                                     >  2:1630300/1‑100 (MQ=255)

                               aTTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTc                                                                     <  2:2641496/100‑1 (MQ=255)

                               aTTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTc                                                                     <  1:3948723/100‑1 (MQ=255)

                                ttATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCg                                                                    <  1:355841/100‑1 (MQ=255)

                                 tATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGc                                                                   <  1:843561/100‑1 (MQ=255)

                                  aTCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCt                                                                  <  1:422124/100‑1 (MQ=255)

                                   tCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCtt                                                                 >  1:2623076/1‑100 (MQ=255)

                                     cagcagGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCg                                                               >  2:3150827/1‑100 (MQ=255)

                                      agcagGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCg                                                               <  2:616490/99‑1 (MQ=255)

                                      agcagGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGa                                                              >  1:131784/1‑100 (MQ=255)

                                       gcagGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGaa                                                             <  2:1378394/100‑1 (MQ=255)

                                          ggCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCgagatcggaag                                                          >  1:3125448/1‑90 (MQ=255)

                                          ggCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGAt                                                          >  1:603502/1‑100 (MQ=255)

                                          ggCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTCATTCGCTTCGAAGAt                                                          >  1:3883431/1‑100 (MQ=255)

                                           gCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGAtt                                                         <  1:2757736/100‑1 (MQ=255)

                                               tggctggcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTAt                                                     <  1:3052909/100‑1 (MQ=255)

                                               tggctggcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTAt                                                     <  2:3402347/100‑1 (MQ=255)

                                               tggctggcAATTAACTACTCTTCCGGAATACGCAACACTCGCCCCGGATAAATTTTATCCGGGCTTTTTCGCATCCGTTTATTCGCTTTGAAGATTTTAt                                                     <  2:217859/100‑1 (MQ=255)

                                                 gctggcAATTAACTACTCTTCCGGCATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTg                                                   <  1:2961403/100‑1 (MQ=255)

                                                  ctggcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGt                                                  <  1:3635708/100‑1 (MQ=255)

                                                     gcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACa                                               <  1:3359905/100‑1 (MQ=255)

                                                      cAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAg                                              <  1:2437024/100‑1 (MQ=255)

                                                        aTTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGAt                                            >  1:1756056/1‑100 (MQ=255)

                                                          tAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTTTCTTTTTTTTTTTCTCTTCGAACATTTTATTTTACAGAt                                            >  2:3209784/1‑98 (MQ=255)

                                                             ctactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCg                                       <  1:1191832/100‑1 (MQ=255)

                                                             ctactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCg                                       <  2:3883401/100‑1 (MQ=255)

                                                             ctactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTAGAAGATTTTATTGTACAGATTTGc                                        >  2:3055540/1‑99 (MQ=255)

                                                              tactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGt                                      <  1:1630309/100‑1 (MQ=255)

                                                                ctctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGtt                                     >  2:2223508/1‑99 (MQ=255)

                                                                ctctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTa                                    >  2:440379/1‑100 (MQ=255)

                                                                 tctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTAc                                   <  1:495001/100‑1 (MQ=255)

                                                                 tctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTAc                                   <  2:131784/100‑1 (MQ=255)

                                                                  ctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGtt                                     <  1:2516792/97‑1 (MQ=255)

                                                                  ctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGCTTATTCGCATCGAAGATTTTATTGTACAGCTTAACATTAcc                                  >  2:235415/1‑100 (MQ=255)

                                                                     ccGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGTGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTa                               <  2:3454425/100‑1 (MQ=255)

                                                                       ggAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGa                             >  2:2099761/1‑100 (MQ=255)

                                                                       ggAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCATTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGa                             >  1:3349166/1‑100 (MQ=255)

                                                                        gAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                            <  2:1324261/100‑1 (MQ=255)

                                                                         aaTACGCAACACTTGCCCCGGATAAATTTTATCCTGTCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTATCGTAGAc                            <  2:3379600/99‑1 (MQ=255)

                                                                         aaTACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAcc                           <  2:3928482/100‑1 (MQ=255)

                                                                         aaTACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAcc                           <  1:2138520/100‑1 (MQ=255)

                                                                             cGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGtt                       >  2:1880982/1‑100 (MQ=255)

                                                                              gCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGttt                      <  2:897715/100‑1 (MQ=255)

                                                                               cAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTg                     <  2:3391121/100‑1 (MQ=255)

                                                                               cAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTg                     <  1:535398/100‑1 (MQ=255)

                                                                                 ctaaaaagcccGGATAAAATTTATCCGGGCTTTTTCGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGCTTAGCGTTACCGTAGACCTGTTTGGa                   >  1:2293467/9‑100 (MQ=38)

                                                                                    cTTGCCCCGGATAAGTTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAAt                >  2:1937027/1‑100 (MQ=255)

                                                                                    cTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCTAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAAt                >  2:843801/1‑100 (MQ=255)

                                                                                    cTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAAt                <  1:2923051/100‑1 (MQ=255)

                                                                                       tccccGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGc             >  2:1796882/2‑100 (MQ=255)

                                                                                       gCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGc             <  2:1566636/100‑1 (MQ=255)

                                                                                       gCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGc             >  1:3436676/1‑100 (MQ=255)

                                                                                        gcccGGATAAAATTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCa            >  2:2282802/2‑100 (MQ=255)

                                                                                         cccGGATAAATTTTATCCGGGCTTTTTCGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCAc           <  2:2448783/100‑1 (MQ=255)

                                                                                         cccGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCAc           >  1:3669631/1‑100 (MQ=255)

                                                                                         cccGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCAc           >  2:2890122/1‑100 (MQ=255)

                                                                                            ggATAAATTTTATCCGGGCTTTTTAGCATCGGTTTACTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCTTAGACCTGTTTGGAAATGGCACt          >  2:3326974/1‑98 (MQ=255)

                                                                                               tAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCagag     >  2:2970789/1‑100 (MQ=255)

                                                                                               tAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCagag     <  1:2950215/100‑1 (MQ=255)

                                                                                                aaaTTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAgt    >  1:3388429/1‑100 (MQ=255)

                                                                                                 aaTTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAgtg   <  1:1477407/100‑1 (MQ=255)

                                                                                                  aTTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAgtgt  >  2:244948/1‑100 (MQ=255)

AATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGT  >  NC_000913/2796276‑2796473

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