BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	207,483	A→C	17.9%	V786V (GTA→GTC)	dnaE →	DNA polymerase III alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	207,483	0	A	C	17.9%	355.7 / 76.1	162	V786V (GTA→GTC)	dnaE	DNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand): ref base A (55/77); new base C (10/19); total (65/97)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

Only 100 of 165 total aligned reads displayed.

TACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGG  >  NC_000913/207385‑207581

tACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTTGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTaa                                                                                                   <  1:15430/100‑1 (MQ=255)

tACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTaa                                                                                                   <  1:169934/100‑1 (MQ=255)

tACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTaa                                                                                                   <  2:3625167/100‑1 (MQ=255)

    gATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAc                                                                                               <  1:2172862/100‑1 (MQ=255)

        tAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAccgcc                                                                                           <  2:1818824/100‑1 (MQ=255)

        tAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAccgcc                                                                                           <  1:3536975/100‑1 (MQ=255)

          aCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGa                                                                                         <  1:1129905/100‑1 (MQ=255)

           cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                        <  1:3431771/100‑3 (MQ=255)

           cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGat                                                                                        >  1:3585802/1‑100 (MQ=255)

           cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGat                                                                                        <  2:3341388/100‑1 (MQ=255)

           cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGat                                                                                        <  2:799190/100‑1 (MQ=255)

           cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGat                                                                                        <  2:1656287/100‑1 (MQ=255)

           cAAATCGCACTCTGCGGCCTATACTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                        <  2:1750844/100‑3 (MQ=255)

            aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcca                                                                                       <  1:2642530/100‑4 (MQ=255)

            aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcca                                                                                       <  1:1821823/100‑4 (MQ=255)

            aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGata                                                                                       <  1:2425158/100‑1 (MQ=255)

            aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGata                                                                                       <  1:3392427/100‑1 (MQ=255)

             aaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                      <  1:2059735/100‑5 (MQ=255)

             aaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                      <  1:1344646/100‑5 (MQ=255)

             aaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGatat                                                                                      <  1:1492297/100‑1 (MQ=255)

              aTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccata                                                                                     <  1:611035/100‑6 (MQ=255)

              aTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccata                                                                                     <  1:2319768/100‑6 (MQ=255)

                 gCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGac                                                                                  >  2:665752/1‑100 (MQ=255)

                 gCACTCTGCGGCCTATGCTGTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTCTCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGCTATgg                                                                                    >  2:1272821/1‑98 (MQ=255)

                  naCTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                 >  1:2941530/2‑100 (MQ=255)

                  cACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                 <  2:2505886/100‑1 (MQ=255)

                    ctctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGacaac                                                                               >  1:2524927/1‑100 (MQ=255)

                    ctctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                 >  2:400641/1‑98 (MQ=255)

                     tctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGacaaca                                                                              >  1:3347329/1‑100 (MQ=255)

                     tctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGacaaca                                                                              >  1:1176847/1‑100 (MQ=255)

                      ctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                 <  2:2057690/96‑1 (MQ=255)

                      ctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAacac                                                                             >  1:933066/1‑100 (MQ=255)

                       cgCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACAcc                                                                            <  1:388830/99‑1 (MQ=255)

                        gCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTTATGACCGCCGATATGGACAACACCg                                                                           <  2:2315626/100‑1 (MQ=255)

                        gCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAacac                                                                             <  2:2107765/98‑1 (MQ=255)

                        gCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCg                                                                           >  2:2468809/1‑100 (MQ=255)

                          ggCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGCTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCgag                                                                         >  2:809744/1‑100 (MQ=255)

                             cTATGCGTTGGTGTAATATCAAACGTTATGGCTGAAAGCGCACTATCATGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCgagacg                                                                      >  2:1589376/1‑98 (MQ=255)

                              tATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGATTTTATGGCGGCGTTAATGGCCGCCGAGAGGTACAACACCGAGAAgg                                                                     >  1:3414865/1‑100 (MQ=255)

                               aTGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggt                                                                    <  2:308605/100‑1 (MQ=255)

                               aTGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggt                                                                    >  1:1011779/1‑100 (MQ=255)

                                   tttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGgtggtg                                                                <  2:3043861/100‑1 (MQ=255)

                                    ttGGTGTTATATCAAACGTTATGGCTGAAAGCGCCCTAGCCTGCGGAGTTTATGGCTGCGGTAAGGACCGCCGATATGGACAACACCGAGAAGGtggtgg                                                               <  2:1258593/100‑1 (MQ=255)

                                    ttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGtggtgg                                                               <  2:2023381/100‑1 (MQ=255)

                                     tGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTggg                                                              <  1:1222223/100‑1 (MQ=255)

                                      gatgtCATATCAAACGTTATGGCTGAAATCTCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTggg                                                              <  2:342817/97‑1 (MQ=255)

                                         gtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtg                                                          <  2:2962744/100‑1 (MQ=255)

                                         gtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCt                                                           <  2:1712704/99‑1 (MQ=255)

                                         gtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGGCtg                                                          >  1:1658220/1‑100 (MQ=255)

                                           cATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggt                                                        >  1:310238/1‑100 (MQ=255)

                                            atatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGtggcggcgctagtg                                                       >  2:3518079/1‑92 (MQ=255)

                                             tatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                      >  1:2523656/1‑100 (MQ=255)

                                                cAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATg                                                   >  1:1689802/1‑100 (MQ=255)

                                                cAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATg                                                   <  1:1859030/100‑1 (MQ=255)

                                                cAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATg                                                   <  1:225944/100‑1 (MQ=255)

                                                  aaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAg                                                 <  1:3216803/100‑1 (MQ=255)

                                                  aaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGCTCTGGTGGAagag                                                 >  2:116595/1‑100 (MQ=255)

                                                    cGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTg                                               >  1:1571366/1‑100 (MQ=255)

                                                        aTGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTgg                                           <  1:1786296/100‑1 (MQ=255)

                                                         tGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGc                                          <  1:2976284/100‑1 (MQ=255)

                                                         tGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGc                                          >  2:2582726/1‑100 (MQ=255)

                                                           gCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGACGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGAGGAGTGCTAGCgg                                        >  2:437035/1‑100 (MQ=255)

                                                            cTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGa                                       >  1:811243/1‑100 (MQ=255)

                                                             tGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGAt                                      <  1:83147/100‑1 (MQ=255)

                                                              gAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATg                                     <  2:2649288/100‑1 (MQ=255)

                                                               aaaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATgg                                    <  1:197987/100‑1 (MQ=255)

                                                               aaaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATgg                                    <  1:197902/100‑1 (MQ=255)

                                                                aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATggg                                   <  1:2789985/100‑1 (MQ=255)

                                                                aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATggg                                   <  2:2318559/100‑1 (MQ=255)

                                                                aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATggg                                   <  2:3347304/100‑1 (MQ=255)

                                                                  gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATggggn                                 <  2:3372314/100‑2 (MQ=255)

                                                                  gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGc                                 <  2:462276/100‑1 (MQ=255)

                                                                  gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGc                                 <  2:518822/100‑1 (MQ=255)

                                                                  gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGGGGGTCTGGTGGATAAGTGc                                              >  2:2376654/1‑87 (MQ=255)

                                                                    gcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTg                               >  1:1862049/1‑100 (MQ=255)

                                                                    gcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTg                               <  1:3373176/100‑1 (MQ=255)

                                                                     cACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGa                              <  2:38079/100‑1 (MQ=255)

                                                                      aCTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGTCTGAGTGCCGGCGGATGGGGCTGa                              >  2:2254592/1‑99 (MQ=255)

                                                                      aCTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGGTGAGTGCTGGCGGATGGGGCTGaa                             <  1:2719741/100‑1 (MQ=255)

                                                                      aCTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGaa                             <  2:324596/100‑1 (MQ=255)

                                                                       cTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGaaa                            <  1:634087/100‑1 (MQ=255)

                                                                           ccTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATcc                        >  1:436219/1‑100 (MQ=255)

                                                                           ccTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATcc                        <  1:3544866/100‑1 (MQ=255)

                                                                             tGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTg                      <  1:1091553/100‑1 (MQ=255)

                                                                             tGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTg                      <  1:472762/100‑1 (MQ=255)

                                                                              gCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGc                     <  1:3078169/100‑1 (MQ=255)

                                                                                 gAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccac                  <  2:883066/100‑1 (MQ=255)

                                                                                   gtccatatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacca                >  1:2377543/9‑100 (MQ=38)

                                                                                    tccatatcggcggTCATTACCGCCGCTATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAg               >  1:2833948/8‑100 (MQ=38)

                                                                                    tccatatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGACTGGTGGATGAGTGCTAGCGGATGGGGCTGAAAATCCTGGCACCAg               >  2:1447322/8‑100 (MQ=255)

                                                                                     ccatatcggcggTCATTACCGCCGATATGGACAACACGGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCGGCCACCAg               >  2:1948798/7‑99 (MQ=255)

                                                                                       atatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGATCTGGCGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCCGata            >  1:879694/5‑100 (MQ=255)

                                                                                            ggcggTCAGTACCGCCGATATGGACAACACCGAGAAGGTGGTGCGTCTGGTGGTTGAGTCCTGGCGGCTGGGGCTGAAGATTCTGCCGCCGGTTATaaa        >  2:3314490/1‑99 (MQ=255)

                                                                                            ggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAAc       <  2:2621254/100‑1 (MQ=255)

                                                                                            ggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAAc       >  1:594531/1‑100 (MQ=255)

                                                                                              cggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTc     >  1:302814/1‑100 (MQ=255)

                                                                                               ggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTcc    <  2:1862040/100‑1 (MQ=255)

                                                                                                 tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg  <  2:2473633/100‑1 (MQ=255)

                                                                                                 tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg  <  1:731126/100‑1 (MQ=255)

                                                                                                 tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg  <  1:1365592/100‑1 (MQ=255)

TACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGG  >  NC_000913/207385‑207581

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