BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,361,697	C→G	7.5%	H174Q (CAC→CAG)	gltF →	periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,361,697	0	C	G	7.5%	369.8 / 17.3	107	H174Q (CAC→CAG)	gltF	periplasmic protein
Reads supporting (aligned to +/- strand): ref base C (54/45); new base G (3/5); total (57/50)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 107 total aligned reads displayed.

GGGAACGAAAAACGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCACCTTCCCGTTA  >  NC_000913/3361550‑3361844

gggAACGAAAAACGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACt                                                                                                                                                    >  2:611655/1‑149 (MQ=255)

 ggAACGAAAAACGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTg                                                                                                                                                   >  1:348904/1‑149 (MQ=255)

    aCGAAAAACGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTc                                                                                                                                                <  1:768217/149‑1 (MQ=255)

      gAAAAACGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTTGACCTGATTCAAACTGACCACTGGTcag                                                                                                                                              >  4:787746/1‑149 (MQ=255)

      gAAAAACGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTcag                                                                                                                                              <  2:348904/149‑1 (MQ=255)

        aaaaCGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCAGTGGTcagtt                                                                                                                                            <  4:732311/149‑3 (MQ=255)

        aaaaCGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCAGTGGTcagtt                                                                                                                                            <  4:640883/149‑3 (MQ=255)

        aaaaCGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCAGTGGTcagtt                                                                                                                                            <  3:260294/149‑3 (MQ=255)

        aaaaCGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTcagca                                                                                                                                            <  1:146992/149‑1 (MQ=255)

               cTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTca                                                                                                                                               >  3:175332/1‑139 (MQ=255)

               cTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTca                                                                                                                                               <  4:175332/139‑1 (MQ=255)

                tttGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGAcccc                                                                                                                                    >  3:602863/1‑149 (MQ=255)

                    gTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCcaca                                                                                                                                >  2:601609/1‑149 (MQ=255)

                      ttAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACGCCCACATg                                                                                                                              <  3:815057/149‑1 (MQ=255)

                        aGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGCCAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGa                                                                                                                            >  1:103749/1‑149 (MQ=255)

                          gTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGaaa                                                                                                                          <  1:268318/149‑1 (MQ=255)

                            tCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGc                                                                                                                        <  3:1882/149‑1 (MQ=255)

                             cTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGcc                                                                                                                       <  1:596558/149‑1 (MQ=255)

                                   cTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTcagcag                                                                                                                                           <  4:200431/123‑1 (MQ=255)

                                   cTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTcagcag                                                                                                                                           >  3:200431/1‑123 (MQ=255)

                                    tAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTa                                                                                                                >  3:321822/1‑149 (MQ=255)

                                    tAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTa                                                                                                                >  2:830687/1‑149 (MQ=255)

                                    tAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTa                                                                                                                >  1:114269/1‑149 (MQ=255)

                                      aTGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAAt                                                                                                              <  3:275234/149‑1 (MQ=255)

                                       tGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAAATCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATg                                                                                                             <  4:602863/149‑1 (MQ=255)

                                         gACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGgc                                                                                                           <  4:826024/149‑1 (MQ=255)

                                         gACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGgc                                                                                                           <  4:76199/149‑1 (MQ=255)

                                          aCAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGgcg                                                                                                          >  3:6532/1‑149 (MQ=255)

                                            aaaaaaTTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATc                                                                                                                  >  4:294142/1‑139 (MQ=255)

                                            aaaaaaTTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATc                                                                                                                  <  3:294142/139‑1 (MQ=255)

                                            aaaaaaTTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCa                                                                                                        >  2:745642/1‑149 (MQ=255)

                                            aaaaaaTTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCa                                                                                                        >  2:341659/1‑149 (MQ=255)

                                            aaaaaaTTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCa                                                                                                        >  1:220965/1‑149 (MQ=255)

                                                  ttGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGc                                                                                                  >  3:800179/1‑149 (MQ=255)

                                                    ggCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTc                                                                                                >  3:542502/1‑149 (MQ=255)

                                                     gcgcTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCg                                                                                               >  2:721183/1‑149 (MQ=255)

                                                       gcTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCt                                                                                             <  1:11956/149‑1 (MQ=255)

                                                         tGCATCCATTGGTATTGACAGGTCTAATGGTGCGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCCTCTAATTGCGCATTCTGCTCGCTGa                                                                                           <  2:103749/149‑1 (MQ=255)

                                                            aTCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACaa                                                                                        <  4:800179/149‑1 (MQ=255)

                                                            aTCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACaa                                                                                        >  1:885782/1‑149 (MQ=255)

                                                            aTCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACaa                                                                                        <  1:745642/149‑1 (MQ=255)

                                                                aTTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGc                                                                                                                        >  4:66744/1‑113 (MQ=255)

                                                                aTTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGc                                                                                                                        <  3:66744/113‑1 (MQ=255)

                                                                aTTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTg                                                                                    >  4:75680/1‑149 (MQ=255)

                                                                aTTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTg                                                                                    >  4:129397/1‑149 (MQ=255)

                                                                   ggTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTc                                                                                 <  1:867429/149‑1 (MQ=255)

                                                                    gTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCa                                                                                <  3:648005/149‑1 (MQ=255)

                                                                        tGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCaa                                                                            <  2:61565/149‑1 (MQ=255)

                                                                          aCAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGc                                                                                                                        <  1:191006/103‑1 (MQ=255)

                                                                          aCAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGc                                                                                                                        >  2:191006/1‑103 (MQ=255)

                                                                          aCAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATc                                                                          >  1:52651/1‑149 (MQ=255)

                                                                            aGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGa                                                                        <  3:813597/149‑1 (MQ=255)

                                                                              gTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCa                                                                                >  1:293198/1‑139 (MQ=255)

                                                                              gTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCa                                                                                <  2:293198/139‑1 (MQ=255)

                                                                              gTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAgc                                                                      >  2:719383/1‑149 (MQ=255)

                                                                                        gggATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTc                                                            <  1:764109/149‑1 (MQ=255)

                                                                                        gggATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTc                                                            >  4:434517/1‑149 (MQ=255)

                                                                                        gggATTCACGAGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGTAGCGACCCCCACATGGAAAGCCGCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTc                                                            >  2:194256/1‑149 (MQ=255)

                                                                                           aTTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGt                                                         <  1:506869/149‑1 (MQ=255)

                                                                                           aTTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCACACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGt                                                         <  1:741453/149‑1 (MQ=255)

                                                                                               aCGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCa                                                     >  4:265230/1‑149 (MQ=255)

                                                                                               aCGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCa                                                     >  2:623681/1‑149 (MQ=255)

                                                                                               aCGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCa                                                     >  1:311886/1‑149 (MQ=255)

                                                                                                cgcgGCGGACGATACAGGTGAAATTCCTGTTGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCaa                                                    <  3:787746/149‑1 (MQ=255)

                                                                                                    gcggACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAAc                                                <  4:542502/149‑1 (MQ=255)

                                                                                                       gACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAgg                                             >  1:809124/1‑149 (MQ=255)

                                                                                                       gACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAgg                                             <  4:704213/149‑1 (MQ=255)

                                                                                                        aCGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAggg                                            <  1:623681/149‑1 (MQ=255)

                                                                                                            tACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAAc                                        <  4:159500/149‑1 (MQ=255)

                                                                                                               aGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGa                                     <  2:114269/149‑1 (MQ=255)

                                                                                                                ggTGAACTTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCGCATGGTAAGCCTCATCTATTGGTGCATTCTTCTCGCTGACCAGTTGCTAAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGac                                    <  1:719383/149‑1 (MQ=255)

                                                                                                                 gTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGaca                                   <  3:853376/149‑1 (MQ=255)

                                                                                                                  tGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGacac                                  >  2:291182/1‑149 (MQ=255)

                                                                                                                  tGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCCCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGacac                                  >  4:533389/1‑149 (MQ=255)

                                                                                                                   gAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACAcc                                 >  1:120980/1‑149 (MQ=255)

                                                                                                                        tCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTg                            >  3:329526/1‑149 (MQ=255)

                                                                                                                           tGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTa                         <  1:286531/149‑1 (MQ=255)

                                                                                                                           tGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCAACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTa                         >  3:264393/1‑149 (MQ=255)

                                                                                                                             tCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTAtt                       <  3:731293/149‑1 (MQ=255)

                                                                                                                                 ccTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAg                   <  1:574579/149‑1 (MQ=255)

                                                                                                                                  cTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACa                                               <  1:231924/120‑1 (MQ=255)

                                                                                                                                  cTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACa                                               >  2:231924/1‑120 (MQ=255)

                                                                                                                                  cTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGc                  >  1:791926/1‑149 (MQ=255)

                                                                                                                                         ctgcTGACCAGTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTAt                        >  1:826377/4‑136 (MQ=255)

                                                                                                                                         ctgcTGACCAGTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTAt                        <  2:826377/133‑1 (MQ=255)

                                                                                                                                         ctgcTGACCAGTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCAcc           >  2:47484/4‑149 (MQ=255)

                                                                                                                                         aaaCTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCAcc           >  4:173030/1‑149 (MQ=255)

                                                                                                                                         aaaCTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCAcc           <  2:311886/149‑1 (MQ=255)

                                                                                                                                          tgcTGACCAGTGGTCAGTAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTa                                                               >  1:287780/3‑96 (MQ=255)

                                                                                                                                          tgcTGACCAGTGGTCAGTAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTa                                                               <  2:287780/94‑1 (MQ=255)

                                                                                                                                          aaCTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCACCt          >  3:497012/1‑149 (MQ=255)

                                                                                                                                           aCTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCACCtt         >  2:636516/1‑149 (MQ=255)

                                                                                                                                              gACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCACCTTccc      >  2:809717/1‑149 (MQ=255)

                                                                                                                                               aCCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCACCTTCCCg     >  4:786284/1‑149 (MQ=255)

                                                                                                                                                 cACTGGTCAGCAGCGACCCCCGCATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGa                                     >  1:641531/1‑115 (MQ=255)

                                                                                                                                                 cACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTAt                        <  4:157385/128‑1 (MQ=255)

                                                                                                                                                 cACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTAt                        >  3:157385/1‑128 (MQ=255)

                                                                                                                                                 cACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCACCTTCCCGtt   >  1:209202/1‑149 (MQ=255)

                                                                                                                                                 cACTGGTCAGCAGCGACCCCCACATGAAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGa                                     <  2:641531/115‑1 (MQ=255)

                                                                                                                                                  aCTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCACCTTCCCGTTa  >  3:328354/1‑149 (MQ=255)

GGGAACGAAAAACGCCTTTGGTTTAGGTTCTGCGCCTAATGGACAAAAAATTGGCGCTGCATCCATTGGTATTGACAGGTCTAATGGTGGGATTCACGCGGCGGACGATACAGGTGAAATTCCTGTCGACCTGATTCAAACTGACCACTGGTCAGCAGCGACCCCCACATGGAAAGCCTCATCTAATGGCGCATTCTGCTCGCTGACAAGTTGCTCAGCAATCGAGCGCGGTTATTCTGTTGCAAAAACAGGGGAACTGACACCTGTTGCTATTACAGCGGTCACCTTCCCGTTA  >  NC_000913/3361550‑3361844

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