BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,965,820	G→A	56.5%	A20A (GCG→GCA)	trxA →	thioredoxin 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,965,820	0	G	A	56.5%	69.6 / 229.0	184	A20A (GCG→GCA)	trxA	thioredoxin 1
Reads supporting (aligned to +/- strand): ref base G (36/44); new base A (62/42); total (98/86)
Fisher's exact test for biased strand distribution p-value = 5.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 186 total aligned reads displayed.

CGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAA  >  NC_000913/3965722‑3965919

   aaCCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACg                                                                                                 <  1:536486/100‑1 (MQ=255)

     ccAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACggg                                                                                               >  1:3762854/1‑100 (MQ=255)

     ccAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACggg                                                                                               >  1:1573849/1‑100 (MQ=255)

     ccAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACggg                                                                                               >  1:2484691/1‑100 (MQ=255)

     ccAACACGCCAGGCTGATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACggg                                                                                               <  1:4078792/100‑1 (MQ=255)

       aaCACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGc                                                                                             <  2:3946978/100‑1 (MQ=255)

       aaCACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGc                                                                                             <  2:1067529/100‑1 (MQ=255)

         ctcGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGa                                                                                           <  1:4136683/98‑1 (MQ=255)

         cacGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGATGACAGTTTTGACACGGATGTACTCAAAGCAGCCGGGGCGa                                                                                           >  2:3235725/1‑100 (MQ=255)

           cGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGaag                                                                                         >  2:4136656/1‑98 (MQ=255)

            gCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATcc                                                                                        <  2:1723376/100‑1 (MQ=255)

            gCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGAGATcc                                                                                        >  2:3060602/1‑100 (MQ=255)

             ccAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCt                                                                                       >  2:2884377/1‑100 (MQ=255)

                 gCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCtcgtc                                                                                   <  2:2929872/100‑1 (MQ=255)

                 gCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCtcgtc                                                                                   >  2:927514/1‑100 (MQ=255)

                  cTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCtcgtcg                                                                                  <  1:72319/100‑1 (MQ=255)

                   ttATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGa                                                                                 >  1:2613844/1‑100 (MQ=255)

                     aTTCCGGTGGAGTTATATATGAGCGATAAAAGTATTCCCCTGACTGACGACAGTTTTGACACGGATGGACTCAAAGCAGACGGGGCGATCCTCGTCGAtt                                                                               <  1:415729/100‑1 (MQ=255)

                      ttCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGAttt                                                                              >  1:1488458/1‑100 (MQ=255)

                      ttCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGAttt                                                                              <  2:2971381/100‑1 (MQ=255)

                        ccTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCt                                                                            >  1:4021854/1‑100 (MQ=255)

                         cTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTg                                                                           <  2:1604159/100‑1 (MQ=255)

                         cTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTg                                                                           >  2:1438003/1‑100 (MQ=255)

                         cTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTg                                                                           >  2:2841024/1‑100 (MQ=255)

                         cTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTg                                                                           <  2:1273826/100‑1 (MQ=255)

                         cTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGAATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTg                                                                           <  2:88017/100‑1 (MQ=255)

                            tgGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGc                                                                        <  1:2766923/100‑1 (MQ=255)

                             ggAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCa                                                                       <  1:2059677/100‑1 (MQ=255)

                                 ttATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAgt                                                                   <  1:3389321/100‑1 (MQ=255)

                                     atatGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAgtggtg                                                               >  1:3182426/1‑100 (MQ=255)

                                      tatGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGc                                                              <  1:3852541/100‑1 (MQ=255)

                                      tatGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGc                                                              <  1:3699025/100‑1 (MQ=255)

                                      tatGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGc                                                              >  2:704474/1‑100 (MQ=255)

                                       atGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCg                                                             >  1:860/1‑100 (MQ=255)

                                       atGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCg                                                             >  1:4109488/1‑100 (MQ=255)

                                       atGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCg                                                             >  1:3659580/1‑100 (MQ=255)

                                        tGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCgg                                                            >  2:1222317/1‑100 (MQ=255)

                                        tGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCgg                                                            >  1:2133147/1‑100 (MQ=255)

                                         gAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGt                                                           <  2:1614227/100‑1 (MQ=255)

                                         gAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGt                                                           >  1:1096343/1‑100 (MQ=255)

                                         gAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGt                                                           >  1:3747017/1‑100 (MQ=255)

                                            cGATAAAATTAATCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCg                                                        >  2:2158354/1‑100 (MQ=255)

                                             gATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGt                                                       <  1:2189264/100‑1 (MQ=255)

                                               tAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGc                                                     >  2:3165461/1‑100 (MQ=255)

                                               tAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGc                                                     >  1:1331268/1‑100 (MQ=255)

                                                 aaaTTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGt                                                       >  2:731969/1‑96 (MQ=255)

                                                 aaaTTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCaa                                                   >  1:3694508/1‑100 (MQ=255)

                                                 aaaTTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCaa                                                   >  1:1630712/1‑100 (MQ=255)

                                                   attattCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCaaaa                                                 <  1:1980484/100‑1 (MQ=255)

                                                    ttattCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAAt                                                <  1:2433162/100‑1 (MQ=255)

                                                     tattCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATg                                               >  1:3193609/1‑100 (MQ=255)

                                                       ttCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGAt                                             <  2:3640144/100‑1 (MQ=255)

                                                         cACCTGACTGACGACAGTTTTGGCACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCg                                           >  2:2394496/1‑100 (MQ=255)

                                                         cACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATc                                            <  2:1028458/99‑1 (MQ=255)

                                                         cACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCg                                           <  1:1267731/100‑1 (MQ=255)

                                                         cACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCg                                           >  1:1140310/1‑100 (MQ=255)

                                                           ccTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGcc                                         <  1:3556345/100‑1 (MQ=255)

                                                           ccTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGcc                                         >  2:617630/1‑100 (MQ=255)

                                                            ctgactgaCGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGccc                                        <  1:2139887/100‑1 (MQ=255)

                                                            ctgactgaCGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGccc                                        >  2:3681277/1‑100 (MQ=255)

                                                             tgactgacGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGcccc                                       <  1:3331266/100‑1 (MQ=255)

                                                             tgactgacGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGcccc                                       <  2:3416416/100‑1 (MQ=255)

                                                             tgactgacGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGcccc                                       <  1:2307272/100‑1 (MQ=255)

                                                             tgactgacGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGcccc                                       <  2:3959323/100‑1 (MQ=255)

                                                              gactgacGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCg                                      <  1:1752539/100‑1 (MQ=255)

                                                                ctgacGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGAt                                    >  1:2932894/1‑100 (MQ=255)

                                                                  gacgacAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTc                                  >  2:1480624/1‑100 (MQ=255)

                                                                   acgacAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCt                                 >  1:3226904/1‑100 (MQ=255)

                                                                   acgacAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCt                                 >  2:2349270/1‑100 (MQ=255)

                                                                    cgacAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTg                                >  1:3031818/1‑100 (MQ=255)

                                                                     gacAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTgg                               <  1:1891592/100‑1 (MQ=255)

                                                                      acAGTGTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGa                              <  2:573540/100‑1 (MQ=255)

                                                                       cAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGAt                             >  1:493775/1‑100 (MQ=255)

                                                                       cAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGAt                             <  1:2334197/100‑1 (MQ=255)

                                                                        aGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATg                            <  2:2027762/100‑1 (MQ=255)

                                                                        aGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATg                            <  1:2264202/100‑1 (MQ=255)

                                                                        aGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATg                            <  1:2841106/100‑1 (MQ=255)

                                                                        aGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTGCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATg                            <  2:2586009/100‑1 (MQ=255)

                                                                         gTTTTCCCACGGCTGTCCTCAACGCGGACGGGGCGATCCTCCTCACTTTCTGGGCCCGGCGGTGCGCTCCGCACAAAATGCTCGCCCCGCTTCTGGa                              >  2:1871224/1‑97 (MQ=255)

                                                                               acacGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGc                     <  2:3694484/100‑1 (MQ=255)

                                                                                cacGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCt                    >  1:1478839/1‑100 (MQ=255)

                                                                                cacGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCt                    >  2:2867837/1‑100 (MQ=255)

                                                                                   ggATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGAc                 <  2:227898/100‑1 (MQ=255)

                                                                                   ggATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGAc                 <  2:3989239/100‑1 (MQ=255)

                                                                                      tGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGaa              >  1:1400072/1‑100 (MQ=255)

                                                                                      tGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGaa              >  2:2755581/1‑100 (MQ=255)

                                                                                        tACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAata            >  1:3123820/1‑100 (MQ=255)

                                                                                         aCTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGCTGAAAGCGCTGACGAata            >  2:1028608/1‑99 (MQ=255)

                                                                                          cTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATc          >  1:2474735/1‑100 (MQ=255)

                                                                                          cTCAAAGCAGACGGGGCGATCCTCGTAGATTTCTCGGCAGAGTGGTGCGGTCCCTGCAACATGGTCGCCCCGATTCTTGATGAAACCGCTGACGaaga            >  2:2687727/1‑96 (MQ=255)

                                                                                           tCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCa         >  1:3173559/1‑100 (MQ=255)

                                                                                            cAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCa         <  2:3878155/99‑1 (MQ=255)

                                                                                            cAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAg        >  1:932405/1‑100 (MQ=255)

                                                                                            cAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTAGATGAAATCGCTGACGAATATc          >  1:2706314/1‑98 (MQ=255)

                                                                                            cAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAg        >  2:3608238/1‑100 (MQ=255)

                                                                                            cAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAg        >  1:1310849/1‑100 (MQ=255)

                                                                                            cAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAg        <  1:2867853/100‑1 (MQ=255)

                                                                                            cAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAg        <  1:2872853/100‑1 (MQ=255)

                                                                                             aaaGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAgg       <  1:2576131/100‑1 (MQ=255)

                                                                                                  ggACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCaaa  >  1:3050059/1‑100 (MQ=255)

CGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGGCAAA  >  NC_000913/3965722‑3965919

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