BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,743,163	T→C	5.4%	E27G (GAG→GGG)	ribC ←	riboflavin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,743,163	0	T	C	5.4%	460.8 / 10.9	147	E27G (GAG→GGG)	ribC	riboflavin synthase
Reads supporting (aligned to +/- strand): ref base T (65/74); new base C (3/5); total (68/79)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 148 total aligned reads displayed.

AAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAG  >  NC_000913/1743016‑1743312

  aTCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTccac                                                                                                                                                   <  1:1505323/150‑1 (MQ=255)

    cGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTccacca                                                                                                                                                 <  2:650144/150‑1 (MQ=255)

      ccAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACcaca                                                                                                                                               <  1:233341/150‑1 (MQ=255)

         aGATTGGTAATGCGTAACGTTTCTTGCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGc                                                                                                                                            <  2:1300920/150‑1 (MQ=255)

           aTTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGctg                                                                                                                                          <  1:1670185/150‑3 (MQ=255)

               gTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGctggacg                                                                                                                                      <  1:708720/150‑7 (MQ=255)

                   tGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAAt                                                                                                                                  >  1:2204694/1‑150 (MQ=255)

                    gCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAAtt                                                                                                                                 >  2:1574285/1‑150 (MQ=255)

                     cGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAAttt                                                                                                                                >  2:929288/1‑150 (MQ=255)

                      gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCATGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTg                                                                                                                               <  1:1821238/150‑1 (MQ=255)

                      gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGTCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGTTCGGGTAACTCCACCACATGCGTACTAAAATTTg                                                                                                                               <  2:509783/150‑1 (MQ=255)

                      gTAACGGTTCCTTCATCAGCTCAAAACTGACATGGTTGCCGTTAAATTCCGGCACGGGCAGGCAGCAACCGTTAAGCGCCACGGAACCACCGGTTTCCAGGCCGTCCAGCATGTGGGCGGGTAACTCCACCACATGCGTACGAAAATTTg                                                                                                                               <  1:817361/150‑1 (MQ=255)

                       tAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTgg                                                                                                                              >  1:823889/1‑150 (MQ=255)

                       tAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTgg                                                                                                                              >  2:416703/1‑150 (MQ=255)

                         aCGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGtt                                                                                                                            <  2:240481/150‑1 (MQ=255)

                              tctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCg                                                                                                                       <  1:1300706/150‑1 (MQ=255)

                               ctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGt                                                                                                                      >  1:1946471/1‑150 (MQ=255)

                                   catcaGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAAt                                                                                                                  >  1:2335417/1‑150 (MQ=255)

                                   catcaGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCATCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAAt                                                                                                                  >  2:2457235/1‑150 (MQ=255)

                                    atcaGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATc                                                                                                                 <  1:1750825/150‑1 (MQ=255)

                                        ggTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGaca                                                                                                             <  2:285982/150‑1 (MQ=255)

                                         gTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGacac                                                                                                            >  1:2447286/1‑150 (MQ=255)

                                           cAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCa                                                                                                          >  2:598595/1‑150 (MQ=255)

                                             aaaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGt                                                                                                        >  1:1317462/1‑150 (MQ=255)

                                             aaaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGt                                                                                                        >  2:2164380/1‑150 (MQ=255)

                                                  gACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGc                                                                                                   <  1:1441480/150‑1 (MQ=255)

                                                    cATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCgg                                                                                                 <  2:1593396/150‑1 (MQ=255)

                                                    cATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCgg                                                                                                 <  2:1493954/150‑1 (MQ=255)

                                                      tGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAAGCGTTATGCGCCTCGGAAGCACCGGTGTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACGAGTTTTGCGGTg                                                                                               >  2:1946945/1‑150 (MQ=255)

                                                      tGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTg                                                                                               >  2:1310308/1‑150 (MQ=255)

                                                      tGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTg                                                                                               <  1:2164380/150‑1 (MQ=255)

                                                        gTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGcc                                                                                             >  1:2266369/1‑150 (MQ=255)

                                                          tGCCGTTAATTTCCGTCACGGTCAGGGAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCt                                                                                           <  2:1657628/150‑1 (MQ=255)

                                                           gCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCgg                                                                                                 <  2:41650/143‑1 (MQ=255)

                                                           gCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCgg                                                                                                 >  1:41650/1‑143 (MQ=255)

                                                           gCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTg                                                                                          >  1:605123/1‑150 (MQ=255)

                                                            ccGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGt                                                                                         >  2:1201162/1‑150 (MQ=255)

                                                             cGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTa                                                                                        <  1:733342/150‑1 (MQ=255)

                                                              gTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTAc                                                                                       <  1:105025/150‑1 (MQ=255)

                                                              gTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTAc                                                                                       <  1:1307218/150‑1 (MQ=255)

                                                              gTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTAc                                                                                       <  1:406534/150‑1 (MQ=255)

                                                              gTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTAc                                                                                       <  2:1788731/150‑1 (MQ=255)

                                                              gTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGGCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTAc                                                                                       <  1:2072006/150‑1 (MQ=255)

                                                                  aTTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATc                                                                                   >  1:2268733/1‑150 (MQ=255)

                                                                   tttCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACAAGTTTTGCGGTGCCCTGTACAATcc                                                                                  >  1:2036300/1‑150 (MQ=255)

                                                                    ttCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATccc                                                                                 >  1:1596222/1‑150 (MQ=255)

                                                                     tCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATcccc                                                                                <  1:2127253/150‑1 (MQ=255)

                                                                      ccGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCg                                                                               <  1:991908/150‑1 (MQ=255)

                                                                      ccGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCg                                                                               <  2:265173/150‑1 (MQ=255)

                                                                          cACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTaaa                                                                           >  2:379483/1‑150 (MQ=255)

                                                                           aCGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAAc                                                                          <  1:899278/150‑1 (MQ=255)

                                                                            cGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACa                                                                         <  2:922509/150‑1 (MQ=255)

                                                                            cGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACa                                                                         <  1:1238813/150‑1 (MQ=255)

                                                                            cGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACa                                                                         >  1:525653/1‑150 (MQ=255)

                                                                             ggTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACAt                                                                        <  1:1419109/150‑1 (MQ=255)

                                                                               tCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCCCCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGGCCTGTCCAATCCCCGTAACCATAc                                                                      >  2:2409563/1‑150 (MQ=255)

                                                                               tCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATAc                                                                      <  1:909163/150‑1 (MQ=255)

                                                                                cAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACt                                                                     <  1:1849970/150‑1 (MQ=255)

                                                                                cAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGGTTCCAGGCCGGCCAGAATGTCGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACt                                                                     <  2:2073460/150‑1 (MQ=255)

                                                                                  ggCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTa                                                                   <  2:339146/150‑1 (MQ=255)

                                                                                  ggCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTa                                                                   >  2:1217214/1‑150 (MQ=255)

                                                                                  ggCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTa                                                                   >  2:1274204/1‑150 (MQ=255)

                                                                                  ggCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTa                                                                   >  1:2054763/1‑150 (MQ=255)

                                                                                        aaCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTc                                                             >  1:1084675/1‑150 (MQ=255)

                                                                                           cGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTg                                                          <  1:1718340/150‑1 (MQ=255)

                                                                                           cGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTg                                                          <  1:1975418/150‑1 (MQ=255)

                                                                                           cGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCCTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTg                                                          <  1:2457235/150‑1 (MQ=255)

                                                                                            attATGCGCCACCGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGa                                                         >  2:2027283/2‑150 (MQ=255)

                                                                                             ttATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGaa                                                        <  1:379483/150‑1 (MQ=255)

                                                                                                  cgcCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAg                                                   >  2:2299962/1‑150 (MQ=255)

                                                                                                     cACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTa                                                <  1:316552/150‑1 (MQ=255)

                                                                                                     cACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACAACTGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTTCAACTCCTGTAATCATTTa                                                <  2:2113726/150‑1 (MQ=255)

                                                                                                         gAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGAc                                            >  2:1251287/1‑150 (MQ=255)

                                                                                                            gCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACAtt                                         >  2:1610861/1‑150 (MQ=255)

                                                                                                              aCCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCt                                       <  1:332059/150‑1 (MQ=255)

                                                                                                               ccGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTg                                      >  1:531537/1‑150 (MQ=255)

                                                                                                               ccGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTg                                      >  2:803344/1‑150 (MQ=255)

                                                                                                                cGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGt                                     <  1:1549767/150‑1 (MQ=255)

                                                                                                                 ggTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACGCCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGtt                                    >  1:483881/1‑150 (MQ=255)

                                                                                                                  gTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTc                                   <  2:2266369/150‑1 (MQ=255)

                                                                                                                    ttCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAg                                 <  1:105601/150‑1 (MQ=255)

                                                                                                                    ttCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAg                                 <  1:598595/150‑1 (MQ=255)

                                                                                                                    ttCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAg                                 <  1:1961352/150‑1 (MQ=255)

                                                                                                                    ttCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGACCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAg                                 <  2:823889/150‑1 (MQ=255)

                                                                                                                        aGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAGATCAGTTAAGACATTCTGTTCAGcaca                             >  1:1928698/1‑150 (MQ=255)

                                                                                                                          gCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAAt                           <  2:330114/150‑1 (MQ=255)

                                                                                                                           ttcgtacgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATa                          >  1:510994/8‑150 (MQ=39)

                                                                                                                           ccGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATa                          >  1:1827425/1‑150 (MQ=255)

                                                                                                                              gtacgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCa                       >  2:1014365/5‑150 (MQ=255)

                                                                                                                              gtacgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCa                       >  2:2185897/5‑150 (MQ=255)

                                                                                                                               ccAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAg                      >  1:1610858/1‑150 (MQ=255)

                                                                                                                                 aGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGt                    >  1:1331204/1‑150 (MQ=255)

                                                                                                                                  gCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTg                   >  1:1960475/1‑150 (MQ=255)

                                                                                                                                   cATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTgg                  <  1:1217214/150‑1 (MQ=255)

                                                                                                                                      gtgGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGaaa               >  1:271363/1‑150 (MQ=255)

                                                                                                                                       tgGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGaacc              >  1:407667/1‑148 (MQ=255)

                                                                                                                                           cGGGTAACTCCACCACATGCGTACGAAAACTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGcc          <  1:1889619/150‑1 (MQ=255)

                                                                                                                                            gggTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGccc         >  2:526977/1‑150 (MQ=255)

                                                                                                                                             ggTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCt        >  1:1741492/1‑150 (MQ=255)

                                                                                                                                                   tCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAg  <  2:1471240/150‑1 (MQ=255)

AAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAG  >  NC_000913/1743016‑1743312

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`