BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	338,334	T→A	54.0%	S4T (TCA→ACA)	yahG →	DUF1116 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	338,334	0	T	A	54.0%	34.2 / 213.7	153	S4T (TCA→ACA)	yahG	DUF1116 family protein
Reads supporting (aligned to +/- strand): ref base T (42/27); new base A (43/39); total (85/68)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

Only 100 of 156 total aligned reads displayed.

TGCTTGATGCCGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCAACTCGACTG  >  NC_000913/338237‑338431

tgctTGATGCCGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGGAGCCTGATGAGCCAGACa                                                                                                 >  2:325018/1‑100 (MQ=255)

 gctTGATGCCGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGacac                                                                                                <  2:1310306/100‑1 (MQ=255)

  ctTGATGCCGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACt                                                                                               >  2:3138671/1‑100 (MQ=255)

   ttGATGCCGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTg                                                                                              >  1:241476/1‑100 (MQ=255)

    tGATGCCGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGt                                                                                             >  1:611281/1‑100 (MQ=255)

    tGATGCCGGAGTGATTCTGGAGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGAAAGCCTGATGAGCCAGTCACTGt                                                                                             >  2:2621071/1‑100 (MQ=255)

      aTGCCGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGttt                                                                                           >  2:3288972/1‑100 (MQ=255)

        gCCGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAg                                                                                         <  1:2923861/100‑1 (MQ=255)

         ccGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGc                                                                                        <  2:3343241/100‑1 (MQ=255)

         ccGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGc                                                                                        >  1:2756397/1‑100 (MQ=255)

          cGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGc                                                                                        <  2:1664732/99‑1 (MQ=255)

          cGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGcc                                                                                       <  1:4048041/100‑1 (MQ=255)

            gAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCaa                                                                                     >  1:3456366/1‑100 (MQ=255)

            gAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCaa                                                                                     >  1:3040555/1‑100 (MQ=255)

             aGTGATTCTGGCGAGCAGCAGCACCAATACCGGATGGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAAc                                                                                    <  2:245011/100‑1 (MQ=255)

              gTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAAcc                                                                                   <  1:76315/100‑1 (MQ=255)

              gTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAAcc                                                                                   <  2:2618412/100‑1 (MQ=255)

                gATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCAt                                                                                 >  1:1041472/1‑100 (MQ=255)

                 actCTGGCGAGCAACAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCAtt                                                                                >  1:2807241/3‑100 (MQ=255)

                 aTTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCAtt                                                                                >  1:4067795/1‑100 (MQ=255)

                 aTTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCAtt                                                                                >  2:197654/1‑100 (MQ=255)

                   tgtGGCGAGCAGCAGCACCAATGCTGGATTGCTGGCGCGTGGATCTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCAtt                                                                                <  2:904885/96‑1 (MQ=255)

                    cTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGaa                                                                             >  1:2571650/1‑100 (MQ=255)

                     tGGCGAGCAGCAGCACCAATACCGGAGTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCGTGATGAGACAGACACTGTTTAGCCAACCATTg                                                                               <  2:627254/97‑1 (MQ=255)

                       gCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACg                                                                           >  2:690591/1‑99 (MQ=255)

                        cGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGtt                                                                         >  2:3461393/1‑100 (MQ=255)

                         gAGCAGCAGCAGCAATGCCGGATTGCTGGCGCGTGAAGTGATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGGCACTGTTTAGCCAACCATTGAACGtta                                                                        <  2:3897523/100‑1 (MQ=255)

                         gAGCAACAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGtta                                                                        <  2:2608136/100‑1 (MQ=255)

                          agcagcagcACCAATACCGGATTGCTGGCGCTTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGttat                                                                       <  1:1787501/100‑1 (MQ=255)

                          agcagcagcACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGttat                                                                       <  2:354933/100‑1 (MQ=255)

                           gcagcagcaCCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGttatt                                                                      >  1:951124/1‑100 (MQ=255)

                            cagcagcaCCAATACCGGATTGCTGGGGCGTGAATTTATCTGCAAAGTGGAGGAAGCCTGCTGAGCCAGACACTGTTTAGCCAACCaaggaa                                                                             >  2:3629026/1‑87 (MQ=255)

                            cagcagcaCCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGttatta                                                                     <  2:321555/100‑1 (MQ=255)

                              gcagcaCCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAAc                                                                   >  1:3716060/1‑100 (MQ=255)

                                   aCCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCgg                                                              >  1:3453147/1‑100 (MQ=255)

                                   aCCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCgg                                                              <  2:2269892/100‑1 (MQ=255)

                                    ccAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGc                                                             >  1:4108637/1‑100 (MQ=255)

                                    ccAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGc                                                             >  1:1519017/1‑100 (MQ=255)

                                    ccAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGc                                                             <  1:3088157/100‑1 (MQ=255)

                                         aCCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGc                                                        >  2:303323/1‑100 (MQ=255)

                                          ccGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGcc                                                       >  1:2666061/1‑100 (MQ=255)

                                          ccGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGcc                                                       >  1:2584939/1‑100 (MQ=255)

                                          ccGGATTGCTGGCGCATGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCAGCGcc                                                       >  1:2402884/1‑100 (MQ=255)

                                                 gCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTa                                                >  1:3525052/1‑100 (MQ=255)

                                                  cTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAg                                               <  2:3571316/100‑1 (MQ=255)

                                                  cTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACTTCGGCATCTCCATGTTTa                                                >  1:608505/1‑99 (MQ=255)

                                                  cTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAg                                               >  1:3663036/1‑100 (MQ=255)

                                                  cTGGCGCGTGAATGTATCTGCAAAGCGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACCGTATTAACGTCGGCATCGCCACGTTTAg                                               >  2:820584/1‑100 (MQ=255)

                                                   tGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCATCCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGc                                              <  2:2774369/100‑1 (MQ=255)

                                                   tGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGc                                              >  2:3378014/1‑100 (MQ=255)

                                                    ggCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCg                                             >  2:2980216/1‑100 (MQ=255)

                                                     gcgcGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCTACCATTGAACGTTATTAACGTCGGCATCGCCATGTGTAGCGa                                            <  2:3202630/100‑1 (MQ=255)

                                                      cgcgTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATGGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGAt                                           >  1:3110997/1‑100 (MQ=255)

                                                        cgTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGa                                         >  1:799103/1‑100 (MQ=255)

                                                         gTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGAc                                        >  1:1228804/1‑100 (MQ=255)

                                                         gTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGAc                                        <  2:1098034/100‑1 (MQ=255)

                                                         gTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGAc                                        <  1:1940212/100‑1 (MQ=255)

                                                          tGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGAcc                                       >  1:2323268/1‑100 (MQ=255)

                                                           gAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCt                                      <  2:629674/100‑1 (MQ=255)

                                                            aaTTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTg                                     >  1:3845309/1‑100 (MQ=255)

                                                             aTTTATCTGCAAAGGGGATGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATGAACGTCGGCATCGCCATGTTTAGTGATGACCt                                      >  2:2029971/1‑98 (MQ=255)

                                                              tttATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGaa                                   >  1:3014488/1‑100 (MQ=255)

                                                               ttATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGaaa                                  <  1:1996499/100‑1 (MQ=255)

                                                                  tCTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAag                               <  1:85308/100‑1 (MQ=255)

                                                                   cTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAagc                              <  2:3845284/100‑1 (MQ=255)

                                                                   cTGCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAagc                              >  2:1699506/1‑100 (MQ=255)

                                                                    tGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAagca                             >  1:31297/1‑100 (MQ=255)

                                                                    tGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAagca                             >  1:492686/1‑100 (MQ=255)

                                                                     gCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAagcag                            <  2:793665/100‑1 (MQ=255)

                                                                     gCAAAGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAagcag                            <  1:17197/100‑1 (MQ=255)

                                                                      cAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAagcagc                           <  2:1625318/100‑1 (MQ=255)

                                                                      cAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGACATGATTAGCGATGACCTGAAAAagcagc                           >  2:237013/1‑100 (MQ=255)

                                                                      cAAAGGGGAGGAAGCCTGATGAGCCAGATACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAagcag                            <  2:2482856/99‑1 (MQ=255)

                                                                       aaaGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAgcagca                          <  1:2363865/100‑1 (MQ=255)

                                                                       aaaGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAgcagca                          <  1:3263521/100‑1 (MQ=255)

                                                                        aaGGGGAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCAt                         <  2:423222/100‑1 (MQ=255)

                                                                          ggggAGGAAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGACGGCATCGCCATATTTAACGATGACCTGAAAAAGCAGCATGt                       >  2:2901386/1‑100 (MQ=255)

                                                                             gaggaAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGa                    >  2:1275298/1‑100 (MQ=255)

                                                                             gaggaAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGa                    <  2:1497661/100‑1 (MQ=255)

                                                                               ggaAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAGGCAGCATGTAGAAg                  >  2:2817094/1‑100 (MQ=255)

                                                                               ggaAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAg                  <  2:2135980/100‑1 (MQ=255)

                                                                               ggaAGCCTGATGAGCCAGTCACTGTTTAGCCAACCAGTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAg                  >  1:2916041/1‑100 (MQ=255)

                                                                               ggaAGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAg                  <  1:3138690/100‑1 (MQ=255)

                                                                                 aaGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTg                <  2:3378582/100‑1 (MQ=255)

                                                                                  aGCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGa               >  1:3969396/1‑100 (MQ=255)

                                                                                   gCCTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGAc              <  1:727612/100‑1 (MQ=255)

                                                                                    ccTGATGAGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACt             <  1:3172239/100‑1 (MQ=255)

                                                                                     cTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTc            >  1:835185/1‑100 (MQ=255)

                                                                                     cTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTc            <  1:3316956/100‑1 (MQ=255)

                                                                                      tgatgaGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCa           >  1:2114010/1‑100 (MQ=255)

                                                                                      tgatgaGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCa           >  2:592284/1‑100 (MQ=255)

                                                                                      tgatgaGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCa           >  2:685272/1‑100 (MQ=255)

                                                                                       gatgaGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCaa          <  2:437353/100‑1 (MQ=255)

                                                                                        atgaGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCAAc         <  2:3716036/100‑1 (MQ=255)

                                                                                        atgaGCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCAAc         <  2:1392057/100‑1 (MQ=255)

                                                                                         tgaGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCAAc         >  2:21825/1‑99 (MQ=255)

                                                                                          tccgatctCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGCTGACCTGAAAAAGCAGCATGTAGAAGTGACTCAACTc       <  2:2218310/93‑1 (MQ=255)

                                                                                          gaGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCAACTc       >  1:2021910/1‑100 (MQ=255)

                                                                                            gCCAGACACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCAACTCGa     >  2:1717308/1‑100 (MQ=255)

                                                                                               aGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCAACTCGACTg  <  1:718986/100‑1 (MQ=255)

TGCTTGATGCCGGAGTGATTCTGGCGAGCAGCAGCACCAATACCGGATTGCTGGCGCGTGAATTTATCTGCAAAGGGGAGGAAGCCTGATGAGCCAGTCACTGTTTAGCCAACCATTGAACGTTATTAACGTCGGCATCGCCATGTTTAGCGATGACCTGAAAAAGCAGCATGTAGAAGTGACTCAACTCGACTG  >  NC_000913/338237‑338431

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