BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,888,030	A→T	7.2%	intergenic (‑39/+31)	rnd ← / ← fadD	ribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,888,030	0	A	T	7.2%	577.9 / 27.1	170	intergenic (‑39/+31)	rnd/fadD	ribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)
Reads supporting (aligned to +/- strand): ref base A (81/75); new base T (7/5); total (89/81)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

Only 100 of 177 total aligned reads displayed.

TAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTAT  >  NC_000913/1887882‑1888178

taaGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGa                                                                                                                                                       <  4:163218/149‑1 (MQ=255)

taaGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGa                                                                                                                                                       <  4:1557956/149‑1 (MQ=255)

taaGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGa                                                                                                                                                       >  3:550395/1‑149 (MQ=255)

    tGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTaa                                                                                                                                                   <  3:467403/149‑1 (MQ=255)

       gcgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCCAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAACACGCCGGATTAA‑CCg                                                                                                                                               <  1:1158158/149‑1 (MQ=255)

        cgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCgg                                                                                                                                              <  2:793511/149‑1 (MQ=255)

           aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAAT‑CCGGCGt                                                                                                                                           <  2:781391/149‑1 (MQ=255)

           aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGt                                                                                                                                           <  2:8518/149‑1 (MQ=255)

           aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGt                                                                                                                                           <  1:1320995/149‑1 (MQ=255)

             gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCt                                                                                                                                         >  1:429727/1‑149 (MQ=255)

             gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCt                                                                                                                                         >  3:1332082/1‑149 (MQ=255)

             gAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCt                                                                                                                                         >  2:1106419/1‑149 (MQ=255)

              aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAAT‑CCGGCGtttt                                                                                                                                        <  2:1239439/149‑4 (MQ=255)

              aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAAT‑CCGGCGttt                                                                                                                                         >  2:662541/1‑146 (MQ=255)

              aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTg                                                                                                                                        <  2:1305018/149‑1 (MQ=255)

              aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTg                                                                                                                                        >  2:567652/1‑149 (MQ=255)

               aCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGGAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTga                                                                                                                                       <  3:349139/149‑1 (MQ=255)

                 aaaTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTgacg                                                                                                                                     >  2:119685/1‑149 (MQ=255)

                  aaTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGG‑TTAATCCGGCGtttttttt                                                                                                                                    <  2:1629486/149‑8 (MQ=255)

                      cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgactga                                                                                                                                <  3:965289/149‑1 (MQ=255)

                      cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgactga                                                                                                                                >  3:268778/1‑149 (MQ=255)

                        gTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACgactgact                                                                                                                              <  4:1232297/149‑1 (MQ=255)

                             cAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCg                                                                                                                                               >  3:496305/1‑127 (MQ=255)

                                 gCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCa                                                                                                                     <  2:429727/149‑1 (MQ=255)

                                 gCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCa                                                                                                                     <  4:701537/149‑1 (MQ=255)

                                          ggAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTAtt                                                                                                            >  2:1648050/1‑149 (MQ=255)

                                             aaGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTc                                                                                                         <  3:288012/149‑1 (MQ=255)

                                                  acggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACtt                                                                                                    <  4:1393435/149‑1 (MQ=255)

                                                  acggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACtt                                                                                                    >  4:722907/1‑149 (MQ=255)

                                                    ggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTg                                                                                                  >  4:1520234/1‑149 (MQ=255)

                                                      acggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGAAATTCAAATTGTATTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGcc                                                                                                >  4:1606983/1‑149 (MQ=255)

                                                      acggCTTCACAAAAAGAAGCCAGCGCATCGTCCGTAGTAATCATTTGGTAATTCAAATTGTTTTCGCTTTAATGGCCGAC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACAACTGACGTAACGCTCAGGCTTTATTGTCCACTTTGcc                                                                                                <  4:3233/149‑1 (MQ=255)

                                                        ggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgc                                                                                              >  3:1613301/1‑149 (MQ=255)

                                                        ggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgc                                                                                              >  3:854210/1‑149 (MQ=255)

                                                        ggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgc                                                                                              >  2:1713182/1‑149 (MQ=255)

                                                          cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑c                                                                                                                                                 >  1:971068/1‑96 (MQ=255)

                                                           ttCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                                                           >  1:1002753/1‑149 (MQ=255)

                                                             cacacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCt                                                                                         <  3:1520234/149‑1 (MQ=255)

                                                             cacacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGAGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCt                                                                                                                 >  1:306089/1‑125 (MQ=255)

                                                             cacacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGAGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCt                                                                                                                 <  2:306089/125‑1 (MQ=255)

                                                               cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTc                                                                                       <  1:1635642/149‑1 (MQ=255)

                                                               cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTc                                                                                       >  4:340480/1‑149 (MQ=255)

                                                                acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                      >  2:1687341/1‑149 (MQ=255)

                                                                acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                      >  3:119831/1‑149 (MQ=255)

                                                                 caAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGt                                                                                     >  3:1241353/1‑149 (MQ=255)

                                                                 caAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGt                                                                                     <  4:1139655/149‑1 (MQ=255)

                                                                  aaaGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTAtt                                                                                                            <  3:889574/125‑1 (MQ=255)

                                                                  aaaGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTAtt                                                                                                            >  4:889574/1‑125 (MQ=255)

                                                                  aaaGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTc                                                                                    >  4:38537/1‑149 (MQ=255)

                                                                    agaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                  >  4:418549/1‑149 (MQ=255)

                                                                     gaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACg                                                                                 >  2:1375171/1‑149 (MQ=255)

                                                                     gaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACg                                                                                 >  3:1440072/1‑149 (MQ=255)

                                                                      aagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGt                                                                                >  3:1259330/1‑149 (MQ=255)

                                                                         ccAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAAt                                                                             >  3:990431/1‑149 (MQ=255)

                                                                          cAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAAtt                                                                            <  1:1281003/149‑1 (MQ=255)

                                                                             cgcATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATtctc                                                                         <  1:792437/149‑1 (MQ=255)

                                                                               cATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgt                                                                       <  4:511530/149‑1 (MQ=255)

                                                                                aTCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGt                                                                                                                                           >  1:73661/1‑80 (MQ=255)

                                                                                aTCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGt                                                                                                                                           <  2:73661/80‑1 (MQ=255)

                                                                                aTCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTACTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgtc                                                                      <  1:89243/149‑1 (MQ=255)

                                                                                       gTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAAT‑CCGGCGtt                                                                                                                                          >  1:1573054/1‑73 (MQ=39)

                                                                                       gTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAAT‑CCGGCGtt                                                                                                                                          <  2:1573054/74‑2 (MQ=39)

                                                                                           tAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                           <  4:119831/149‑1 (MQ=255)

                                                                                            aaTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTc                                                          <  4:1641687/149‑1 (MQ=255)

                                                                                             aTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTcc                                                         <  4:854210/149‑1 (MQ=255)

                                                                                              tCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCg                                                        >  1:4230/1‑149 (MQ=255)

                                                                                                   tGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGtt                                                   <  4:636835/149‑1 (MQ=255)

                                                                                                   tGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGtt                                                   <  4:1620000/149‑1 (MQ=255)

                                                                                                        aTTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAtt                                              <  4:1311611/149‑1 (MQ=255)

                                                                                                          tCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAACAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTc                                            <  4:526291/149‑1 (MQ=255)

                                                                                                             aaTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                         >  3:686005/1‑149 (MQ=255)

                                                                                                             aaTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                         >  3:1665934/1‑149 (MQ=255)

                                                                                                             aaTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                         <  2:4230/149‑1 (MQ=255)

                                                                                                              aTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGc                                                                                                                                             <  2:671788/48‑1 (MQ=255)

                                                                                                              aTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGc                                                                                                                                             >  1:671788/1‑48 (MQ=255)

                                                                                                              aTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTaa                                                                              <  3:586177/111‑1 (MQ=255)

                                                                                                              aTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTa                                        >  1:1016004/1‑149 (MQ=255)

                                                                                                               ttGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTaa                                       <  3:1593668/149‑1 (MQ=255)

                                                                                                                  ttttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtc                                    >  4:172608/1‑149 (MQ=255)

                                                                                                                   tttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtca                                   <  4:1025733/149‑1 (MQ=255)

                                                                                                                    ttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcat                                  <  4:268778/149‑1 (MQ=255)

                                                                                                                       tctTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc                               >  4:96072/1‑149 (MQ=255)

                                                                                                                         tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                  >  4:1292004/1‑96 (MQ=255)

                                                                                                                         tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                  <  3:1292004/96‑1 (MQ=255)

                                                                                                                         tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCg                                           >  1:1146169/1‑135 (MQ=255)

                                                                                                                           tAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaaa                           <  2:994824/149‑1 (MQ=255)

                                                                                                                           tAGTGAGCGTCAAAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaa                            <  2:718734/149‑1 (MQ=255)

                                                                                                                                     tcgaaaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCATAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAg                  >  4:555243/4‑149 (MQ=255)

                                                                                                                                        aaaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccac                     <  2:925692/143‑1 (MQ=255)

                                                                                                                                         aaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                     <  2:1134923/126‑1 (MQ=255)

                                                                                                                                         aaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                     >  1:1134923/1‑126 (MQ=255)

                                                                                                                                         aaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc              >  3:930106/1‑149 (MQ=255)

                                                                                                                                         aaaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCACAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc              >  1:1129374/1‑149 (MQ=255)

                                                                                                                                          aaaaaCGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                  >  1:1467181/1‑80 (MQ=255)

                                                                                                                                             gaCGCCGGTTAAT‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa          >  4:1068022/2‑149 (MQ=255)

                                                                                                                                             gaCGCCGG‑TTAATCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa          >  3:821627/2‑149 (MQ=255)

                                                                                                                                              aCGCCGGTTAAT‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTAc         >  1:619138/1‑149 (MQ=255)

                                                                                                                                                gCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACtt       <  1:1462467/149‑1 (MQ=255)

                                                                                                                                                  cGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                           >  4:494478/1‑95 (MQ=255)

                                                                                                                                                     aTTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTAt  <  2:1016004/149‑1 (MQ=255)

TAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAA‑CCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTAT  >  NC_000913/1887882‑1888178

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