BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,857,315	C→T	Q734* (CAG→TAG)	mutS →	MutS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,857,315	0	C	T	99.1%	877.8 / ‑3.3	233	Q734* (CAG→TAG)	mutS	MutS
Reads supporting (aligned to +/- strand): ref base C (0/2); new base T (138/93); total (138/95)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 236 total aligned reads displayed.

GGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAG  >  NC_000913/2857225‑2857404

  tCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGtt                                                                                       <  1:4538877/93‑1 (MQ=255)

    tGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTAc                                                                                     <  1:6103546/93‑1 (MQ=255)

    tGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTAc                                                                                     <  2:4487767/93‑1 (MQ=255)

     gTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTAcc                                                                                    <  2:5940755/93‑1 (MQ=255)

         cTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGgag                                                                                >  2:2613467/1‑93 (MQ=255)

         cTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGgag                                                                                >  2:5871333/1‑93 (MQ=255)

            gcgtgggcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGaaa                                                                             <  1:51459/93‑1 (MQ=255)

              gtgggcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGGTACCCACTATTTCGAGCTGTCCTAGCTACCGGAGaaaag                                                                           >  2:399396/1‑92 (MQ=255)

              gtgggcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAAt                                                                           <  1:1440479/93‑1 (MQ=255)

                gggcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATgg                                                                         <  2:4841826/93‑1 (MQ=255)

                gggcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATgg                                                                         >  2:4622558/1‑93 (MQ=255)

                 ggcgtgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGa                                                                        <  2:4370525/93‑1 (MQ=255)

                     tgCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGc                                                                    >  2:367047/1‑93 (MQ=255)

                          ggAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGc                                                               <  1:175143/93‑1 (MQ=255)

                             aaaTCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTac                                                            >  2:1345787/1‑92 (MQ=255)

                             aaaTCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTaa                                                            <  1:4449951/93‑1 (MQ=255)

                               aTCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACg                                                          >  1:6139508/1‑93 (MQ=255)

                               aTCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACg                                                          >  1:2845313/1‑93 (MQ=255)

                               aTCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACg                                                          <  2:704666/93‑1 (MQ=255)

                                tCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGt                                                         <  2:2077707/93‑1 (MQ=255)

                                   ggCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTCACGTGCa                                                      >  2:392820/1‑93 (MQ=255)

                                   ggCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCa                                                      >  1:2581450/1‑93 (MQ=255)

                                      gAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCAtct                                                   >  1:126353/1‑93 (MQ=255)

                                      gAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCAtct                                                   >  2:1766877/1‑93 (MQ=255)

                                        aTAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCg                                                 >  1:891342/1‑93 (MQ=255)

                                        aTAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCg                                                 >  1:5902523/1‑93 (MQ=255)

                                        aTAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCg                                                 >  1:1658206/1‑93 (MQ=255)

                                         tAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGa                                                >  1:3446421/1‑93 (MQ=255)

                                         tAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGa                                                >  2:1250421/1‑93 (MQ=255)

                                         tAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGa                                                <  2:3286657/93‑1 (MQ=255)

                                         tAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGa                                                <  2:3099827/93‑1 (MQ=255)

                                         tAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGa                                                >  2:5183787/1‑93 (MQ=255)

                                           aGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATg                                              <  2:551837/93‑1 (MQ=255)

                                           aGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATg                                              >  2:4129142/1‑93 (MQ=255)

                                           aGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATg                                              >  2:5025772/1‑93 (MQ=255)

                                            gATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGc                                             >  2:1714425/1‑93 (MQ=255)

                                             aTTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCa                                            >  1:4244902/1‑93 (MQ=255)

                                             aTTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCa                                            >  1:4137599/1‑93 (MQ=255)

                                              ttAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCAc                                           >  1:3041423/1‑93 (MQ=255)

                                                  ggCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGa                                       <  1:1100423/93‑1 (MQ=255)

                                                  ggCATTGACGTTATTTGCTACCCACTATTCCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGa                                       <  1:4674259/93‑1 (MQ=255)

                                                   gCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAg                                      <  2:3334136/93‑1 (MQ=255)

                                                    cATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGc                                     >  1:1950779/1‑93 (MQ=255)

                                                     aTTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGGATCTCGATGCACTGGAGCa                                    >  1:5672637/1‑93 (MQ=255)

                                                     aTTGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCa                                    >  1:5025426/1‑93 (MQ=255)

                                                      ttGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTTGAGCAc                                   >  1:5470482/1‑93 (MQ=255)

                                                      ttGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCAc                                   <  2:5045340/93‑1 (MQ=255)

                                                      ttGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCAc                                   <  1:970743/93‑1 (MQ=255)

                                                      ttGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCAc                                   >  1:2717767/1‑93 (MQ=255)

                                                      ttGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCAc                                   >  1:4500902/1‑93 (MQ=255)

                                                      ttGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCAc                                   >  1:3204728/1‑93 (MQ=255)

                                                      ttGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCAc                                   >  1:4600103/1‑93 (MQ=255)

                                                       tGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACg                                  >  1:1008600/1‑93 (MQ=255)

                                                       tGACGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACg                                  >  2:2929992/1‑93 (MQ=255)

                                                         aCGTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGc                                >  2:5759638/1‑93 (MQ=255)

                                                           gTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGa                              >  1:3421855/1‑93 (MQ=255)

                                                           gTTATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGa                              >  2:1967906/1‑93 (MQ=255)

                                                            ttATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGac                             >  1:723595/1‑93 (MQ=255)

                                                            ttATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGac                             >  2:2810817/1‑93 (MQ=255)

                                                             tATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGaca                            >  2:1019891/1‑93 (MQ=255)

                                                             tATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGaca                            >  1:3572418/1‑93 (MQ=255)

                                                             tATTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGaca                            >  2:3777209/1‑93 (MQ=255)

                                                              aTTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGacac                           >  1:4092019/1‑93 (MQ=255)

                                                              aTTTGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGacac                           >  1:4958475/1‑93 (MQ=255)

                                                                ttGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCa                         >  2:725783/1‑93 (MQ=255)

                                                                 tGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCAt                        >  2:218146/1‑93 (MQ=255)

                                                                 tGCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCAt                        <  2:5368689/93‑1 (MQ=255)

                                                                  gCTACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCAtt                       <  2:5285823/93‑1 (MQ=255)

                                                                    tACCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGc                     >  2:3115874/1‑93 (MQ=255)

                                                                     aCCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGcc                    <  2:1899338/93‑1 (MQ=255)

                                                                     aCCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGcc                    <  1:3903175/93‑1 (MQ=255)

                                                                     aCCCACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCAGTGcc                    >  2:1275211/1‑93 (MQ=255)

                                                                        cACTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCttt                 >  1:2304678/1‑93 (MQ=255)

                                                                         ccTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTa                <  2:1082455/92‑1 (MQ=255)

                                                                         aCTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTa                >  1:3538137/1‑93 (MQ=255)

                                                                          cTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTAt               >  2:2375822/1‑93 (MQ=255)

                                                                          cTATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTAt               <  1:1584072/93‑1 (MQ=255)

                                                                           tATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATg              >  2:3435927/1‑93 (MQ=255)

                                                                           tATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATg              <  2:421471/93‑1 (MQ=255)

                                                                           tATTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATg              >  1:5947191/1‑93 (MQ=255)

                                                                            aTTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGc             >  2:1437636/1‑93 (MQ=255)

                                                                            aTTTCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGc             <  2:2159382/93‑1 (MQ=255)

                                                                             tttCGAGCTGACCTAGTTACCGGAGACAATGGAAGTCGTCGCTAACGTGCATCTCGATGCACTGGAGAACGGCGACACCATTGCCTTTATGca            >  1:3395769/1‑93 (MQ=255)

                                                                             tttCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGGGCACGGCGACACCATTGCCTTTATGca            >  2:5512559/1‑93 (MQ=255)

                                                                             tttCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGca            <  2:5676565/93‑1 (MQ=255)

                                                                              ttCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGcac           >  1:3858776/1‑93 (MQ=255)

                                                                              ttCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGcac           >  1:238873/1‑93 (MQ=255)

                                                                              ttCGAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGcac           >  1:3724936/1‑93 (MQ=255)

                                                                                 gAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGc        <  2:5950115/93‑1 (MQ=255)

                                                                                 gAGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGc        >  1:2179508/1‑93 (MQ=255)

                                                                                  aGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCg       <  1:5238502/93‑1 (MQ=255)

                                                                                  aGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCg       <  1:2799649/93‑1 (MQ=255)

                                                                                  aGCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCg       <  1:3274795/93‑1 (MQ=255)

                                                                                   gCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCgg      >  1:1359058/1‑92 (MQ=255)

                                                                                   gCTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGt      <  1:5318515/93‑1 (MQ=255)

                                                                                    gctgcccAATTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTg     <  2:5701599/89‑1 (MQ=255)

                                                                                    cTGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTg     <  1:603629/93‑1 (MQ=255)

                                                                                     tGACCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGc    >  2:2504315/1‑93 (MQ=255)

                                                                                       aCCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAg  >  1:2900504/1‑93 (MQ=255)

                                                                                       aCCTAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAg  <  2:2726199/93‑1 (MQ=255)

GGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTGACCCAGTTACCGGAGAAAATGGAAGGCGTCGCTAACGTGCATCTCGATGCACTGGAGCACGGCGACACCATTGCCTTTATGCACAGCGTGCAG  >  NC_000913/2857225‑2857404

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