BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,856,333	T→G	5.0%	E199A (GAG→GCG) ‡	ydjH ←	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,856,333	0	T	G	5.0%	608.4 / 23.3	179	E199A (GAG→GCG) ‡	ydjH	putative kinase
Reads supporting (aligned to +/- strand): ref base T (101/69); new base G (5/4); total (106/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-03

Only 100 of 180 total aligned reads displayed.

GCGACCGCCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTACC  >  NC_000913/1856186‑1856481

    ccgccgGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCa                                                                                                                                                 >  3:1088471/1‑149 (MQ=255)

       ccgGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCaaaa                                                                                                                                              <  2:1564765/149‑1 (MQ=255)

        cgGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAAt                                                                                                                                             >  4:1357602/1‑149 (MQ=255)

         ggCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAAtt                                                                                                                                            <  2:810680/149‑1 (MQ=255)

          gCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTATCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTa                                                                                                                                           <  4:602000/149‑1 (MQ=255)

               ttCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGaaa                                                                                                                                      >  4:363988/1‑149 (MQ=255)

               ttCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGaaa                                                                                                                                      >  1:807834/1‑149 (MQ=255)

               ttCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGaaa                                                                                                                                      >  2:726854/1‑149 (MQ=255)

                  aTCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAg                                                                                                                                   >  1:1098460/1‑149 (MQ=255)

                      tCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATaa                                                                                                                               >  4:1489/1‑149 (MQ=255)

                              acaaggttGATAAAGCAGCCGTCTTTACCCGTTTTAATAACCACCGTTTTTACGCCGCACGCTAGAAAGCAGTCAGCAATTGCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAAACGACATa                                                                                                                       >  4:1096760/7‑149 (MQ=255)

                               cacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATaa                                                                                                                      >  1:69291/1‑149 (MQ=255)

                                acGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAAc                                                                                                                     >  3:569912/1‑149 (MQ=255)

                                  gCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTc                                                                                                                   >  2:1203374/1‑149 (MQ=255)

                                     tGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAAt                                                                                                                <  1:1680278/149‑1 (MQ=255)

                                     tGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAAt                                                                                                                <  1:1240598/149‑1 (MQ=255)

                                      gATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATg                                                                                                               >  1:1228335/1‑149 (MQ=255)

                                       aTAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGc                                                                                                              >  4:731407/1‑149 (MQ=255)

                                         aaaGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCtt                                                                                                            >  4:1045453/1‑149 (MQ=255)

                                          aaGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                                                                           >  1:1032839/1‑149 (MQ=255)

                                          aaGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                                                                           >  1:749840/1‑149 (MQ=255)

                                          aaGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                                                                           <  3:1371593/149‑1 (MQ=255)

                                          aaGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                                                                           <  1:347746/149‑1 (MQ=255)

                                            gcagcCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAAt                                                                                                                              >  1:644680/1‑128 (MQ=255)

                                            gcagcCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAAt                                                                                                                              <  2:644680/128‑1 (MQ=255)

                                              agcCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCaa                                                                                                       >  2:1158044/1‑149 (MQ=255)

                                               gcCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCaaa                                                                                                      >  2:471722/1‑149 (MQ=255)

                                                 cGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAata                                                                                                    <  3:1707038/149‑1 (MQ=255)

                                                 cGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAata                                                                                                    >  3:1197955/1‑149 (MQ=255)

                                                  gTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAatat                                                                                                   >  3:116055/1‑149 (MQ=255)

                                                  gTCTTTAACCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAatat                                                                                                   >  2:428761/1‑149 (MQ=255)

                                                   tCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatc                                                                                                  <  4:40193/149‑1 (MQ=255)

                                                     tttACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcat                                                                                                >  2:666413/1‑149 (MQ=255)

                                                      ttACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcatc                                                                                               <  2:1465045/149‑1 (MQ=255)

                                                      ttACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCAAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcatc                                                                                               <  4:1197955/149‑1 (MQ=255)

                                                        aCCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCa                                                                                             <  2:1017089/149‑1 (MQ=255)

                                                        aCCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCa                                                                                             >  4:802312/1‑149 (MQ=255)

                                                             ttttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAAt                                                                                                                              <  1:5737/111‑1 (MQ=255)

                                                             ttttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGtt                                                                                        >  4:694122/1‑149 (MQ=255)

                                                              tttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGttt                                                                                       <  2:1681274/149‑1 (MQ=255)

                                                                   tCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCaa                                                                                                                 >  2:1633866/1‑118 (MQ=255)

                                                                   tCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcattc                                                                                  >  4:255988/1‑149 (MQ=255)

                                                                   tCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcattc                                                                                  <  1:324543/149‑1 (MQ=255)

                                                                     accaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCaa                                                                                >  2:1571371/1‑149 (MQ=255)

                                                                     accaccGGTTGTACGCCTCACGCCAGAAAGCAGTCAGCAAATTCATCCAGTGTCTCTCTCCCGGTGGGTAATTTTGCCTCAGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGCTTCATTCaa                                                                                >  4:1229200/1‑149 (MQ=255)

                                                                      ccaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAAc                                                                               <  3:515970/149‑1 (MQ=255)

                                                                       caccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                                                                           <  4:325209/104‑1 (MQ=255)

                                                                       caccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                                                                           >  3:325209/1‑104 (MQ=255)

                                                                       caccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAAt                                                                                                                >  3:992177/1‑115 (MQ=255)

                                                                       caccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAAcc                                                                              >  4:778709/1‑149 (MQ=255)

                                                                       caccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAAcc                                                                              >  3:1134822/1‑149 (MQ=255)

                                                                            tttttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGt                                                                         >  3:1410189/1‑149 (MQ=255)

                                                                             ttttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGtt                                                                        <  4:98967/149‑1 (MQ=255)

                                                                                 aCGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAt                                                                    <  3:555755/149‑1 (MQ=255)

                                                                                 aCGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAt                                                                    <  1:1696726/149‑1 (MQ=255)

                                                                                       cacgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATc                                                              <  2:1736280/149‑1 (MQ=255)

                                                                                          gcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGc                                                           >  2:310776/1‑149 (MQ=255)

                                                                                            aaGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTg                                                                      <  3:1470502/136‑1 (MQ=255)

                                                                                            aaGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTg                                                                      >  4:1470502/1‑136 (MQ=255)

                                                                                                aaGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAAcc                                                                              <  4:470518/124‑1 (MQ=255)

                                                                                                  gCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATaa                                                   <  4:1384716/149‑1 (MQ=255)

                                                                                                     gTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAatca                                                <  2:317383/149‑1 (MQ=255)

                                                                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAt                                                                    <  3:699391/126‑1 (MQ=255)

                                                                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCt                                             >  1:1034090/1‑149 (MQ=255)

                                                                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCt                                             >  3:403782/1‑149 (MQ=255)

                                                                                                         gCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTg                                            >  1:1496197/1‑149 (MQ=255)

                                                                                                         gCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTg                                            >  2:1440421/1‑149 (MQ=255)

                                                                                                          cAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGc                                           >  3:1602604/1‑149 (MQ=255)

                                                                                                          cAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGc                                           >  2:986431/1‑149 (MQ=255)

                                                                                                             tttCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCgg                                        <  1:262499/149‑1 (MQ=255)

                                                                                                             tttCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCgg                                        >  1:709588/1‑149 (MQ=255)

                                                                                                             tttCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCgg                                        >  2:334855/1‑149 (MQ=255)

                                                                                                               acATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCCGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGc                                      <  4:759129/148‑1 (MQ=255)

                                                                                                                cATCCAGTGGCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCCTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCt                                     <  2:269516/149‑1 (MQ=255)

                                                                                                                  tCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCttt                                   >  3:1269553/1‑149 (MQ=255)

                                                                                                                  tCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCttt                                   >  4:1452300/1‑149 (MQ=255)

                                                                                                                   ccAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTa                                  >  1:1704301/1‑149 (MQ=255)

                                                                                                                      gtgtCTCTTTCCCGGTGAGTAATTTTGCCTCGGCACAATTAGGAAACAGATACTCGACATAACTCAATGCTTCGCGAATATCATCCAGCGTTTCATTCAACCGCGGCTTGATCATATCGACACAGATAATCATCTGCCGGGCTTTAGCt                               >  4:1131443/1‑149 (MQ=255)

                                                                                                                      gtgtCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                                                                           <  1:1143518/73‑1 (MQ=255)

                                                                                                                            ctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTa                         >  1:419555/1‑149 (MQ=255)

                                                                                                                             tttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAatcat                                               >  3:780766/1‑126 (MQ=255)

                                                                                                                             tttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTaa                        >  1:1381013/1‑149 (MQ=255)

                                                                                                                                   ggTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATGACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAAttt                  <  3:1458359/149‑1 (MQ=255)

                                                                                                                                   ggTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAAttt                  <  3:1096760/149‑1 (MQ=255)

                                                                                                                                     tGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCt                <  2:15911/149‑1 (MQ=255)

                                                                                                                                     tGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCt                >  3:982863/1‑149 (MQ=255)

                                                                                                                                     tGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCt                <  1:666413/149‑1 (MQ=255)

                                                                                                                                       cctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt              >  2:628168/3‑149 (MQ=255)

                                                                                                                                       aGTAATTTTGCCTCGGAAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt              <  4:1204563/149‑1 (MQ=255)

                                                                                                                                        gTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGtt             <  2:1032839/149‑1 (MQ=255)

                                                                                                                                         tAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTa            <  3:1045453/149‑1 (MQ=255)

                                                                                                                                         tAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTa            >  3:907231/1‑149 (MQ=255)

                                                                                                                                         tAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt              >  4:1580607/1‑147 (MQ=255)

                                                                                                                                         tAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt              <  3:1580607/147‑1 (MQ=255)

                                                                                                                                         tAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTa            >  2:1352804/1‑149 (MQ=255)

                                                                                                                                           aTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAgc          <  3:275455/149‑1 (MQ=255)

                                                                                                                                             tttGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAgcgc        <  2:237894/149‑1 (MQ=255)

                                                                                                                                                gCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGtt             >  1:1586356/1‑141 (MQ=255)

                                                                                                                                                 ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTa    <  1:986431/149‑1 (MQ=255)

                                                                                                                                                   tCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTGCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTAcc  <  1:527697/149‑1 (MQ=255)

GCGACCGCCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTACC  >  NC_000913/1856186‑1856481

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