BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,834,958	G→T	57.2%	Q48K (CAG→AAG)	nepI ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,834,958	0	G	T	57.2%	21.3 / inf	299	Q48K (CAG→AAG)	nepI	putative transporter
Reads supporting (aligned to +/- strand): ref base G (75/53); new base T (90/81); total (165/134)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 300 total aligned reads displayed.

AGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTA  >  CP009273/3834810‑3835106

agCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTgg                                                                                                                                                     >  2:2966983/1‑150 (MQ=255)

   gTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGcc                                                                                                                                                  <  2:1245193/150‑1 (MQ=255)

     cAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCAt                                                                                                                                                <  2:1380247/150‑1 (MQ=255)

      aGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCAtt                                                                                                                                               <  2:820033/150‑1 (MQ=255)

      aGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCAtt                                                                                                                                               >  1:3638676/1‑150 (MQ=255)

         aaaaCGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGc                                                                                                                                            >  2:2219351/1‑150 (MQ=255)

         aaaaCGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGc                                                                                                                                            >  2:4007286/1‑150 (MQ=255)

         aaaaCGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGc                                                                                                                                            <  1:1774520/150‑1 (MQ=255)

         aaaaCGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGc                                                                                                                                            <  2:105934/150‑1 (MQ=255)

           aaCGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGt                                                                                                                                          <  2:2563375/150‑1 (MQ=255)

                caaacaaaATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTcaaca                                                                                                                                     <  1:2045660/150‑1 (MQ=255)

                caaacaaaATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTcaaca                                                                                                                                     <  1:2222250/150‑1 (MQ=255)

                 aaacaaaATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTcaacaa                                                                                                                                    >  2:446453/1‑150 (MQ=255)

                 aaacaaaATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTcaacaa                                                                                                                                    >  2:2740767/1‑150 (MQ=255)

                  aacaaaATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACaaa                                                                                                                                   >  2:172742/1‑150 (MQ=255)

                    caaaATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAACTGCCTAAATCCTTGGCCATTGGCGTCAACAAACt                                                                                                                                 <  2:4563129/150‑1 (MQ=255)

                    caaaATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACt                                                                                                                                 <  1:4613415/150‑1 (MQ=255)

                     aaaaTAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTg                                                                                                                                >  1:1432225/1‑150 (MQ=255)

                      aaaTAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGa                                                                                                                               <  1:1382904/150‑1 (MQ=255)

                       aaTAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGAc                                                                                                                              <  1:3208759/150‑1 (MQ=255)

                            caacGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAACCTGGCAAACATTGCCACAACGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCa                                                                                                                         <  2:173292/150‑1 (MQ=255)

                            caacGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCa                                                                                                                         >  2:1295705/1‑150 (MQ=255)

                             aacGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAACCTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCaa                                                                                                                        <  2:3078289/150‑1 (MQ=255)

                             aacGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCaa                                                                                                                        <  2:3676644/150‑1 (MQ=255)

                              acGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCaaa                                                                                                                       <  1:2390630/150‑1 (MQ=255)

                                gTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCaaaaa                                                                                                                     <  2:2353167/150‑1 (MQ=255)

                                 tAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAAc                                                                                                                    <  2:283512/150‑1 (MQ=255)

                                  aGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACt                                                                                                                   >  2:537854/1‑150 (MQ=255)

                                    cggcggTCAGTAGCCTGAATTGTCTGGGAAATAAACAAACTGGCAAACATTGCCACAAAGTCTGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCAGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACtaa                                                                                                                 >  2:939065/1‑148 (MQ=255)

                                     ggcggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCaa                                                                                                                >  1:2455207/1‑150 (MQ=255)

                                     ggcggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCaa                                                                                                                >  1:2750386/1‑150 (MQ=255)

                                     ggcggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCaa                                                                                                                >  1:329491/1‑150 (MQ=255)

                                         gTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGAt                                                                                                            >  1:4042940/1‑150 (MQ=255)

                                         gTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGAt                                                                                                            >  1:1106354/1‑150 (MQ=255)

                                          tCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAAACCTGGGCCATTTGCGTCAACAAACTTACGGGCAAAAAATAAACGATa                                                                                                           >  1:4349499/1‑150 (MQ=255)

                                            aGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAAt                                                                                                         >  2:2484117/1‑150 (MQ=255)

                                              tAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCa                                                                                                       >  1:1049786/1‑150 (MQ=255)

                                               aGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAg                                                                                                      <  2:907445/150‑1 (MQ=255)

                                               aGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAg                                                                                                      <  2:4228048/150‑1 (MQ=255)

                                               aGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAg                                                                                                      <  1:3307545/150‑1 (MQ=255)

                                                gCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGa                                                                                                     <  2:2480360/150‑1 (MQ=255)

                                                   tGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAg                                                                                                  >  1:2661493/1‑150 (MQ=255)

                                                   tGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAg                                                                                                  >  1:343470/1‑150 (MQ=255)

                                                        tGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGTCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGac                                                                                             >  1:975273/1‑150 (MQ=255)

                                                        tGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGac                                                                                             >  2:3363500/1‑150 (MQ=255)

                                                         gTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGaca                                                                                            <  1:3788000/150‑1 (MQ=255)

                                                          tCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGacac                                                                                           >  2:1982167/1‑150 (MQ=255)

                                                           cTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGAcaca                                                                                          >  1:3122260/1‑150 (MQ=255)

                                                            tGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACaa                                                                                         >  2:470549/1‑150 (MQ=255)

                                                            tGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACaa                                                                                         >  2:2789219/1‑150 (MQ=255)

                                                            tGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACaa                                                                                         >  1:762329/1‑150 (MQ=255)

                                                            tGGGTAATAAACAAACTGGCAAACAGTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACaa                                                                                         <  2:2884339/150‑1 (MQ=255)

                                                              ggTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACaaaa                                                                                       >  2:2292767/1‑150 (MQ=255)

                                                               gTAATAAACAAACTTTCAAACATTGCCACTAAGGAGGTCACGGTCTCCGATTGCCCGTCAAACCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTTACGGTCAATAACTCAACGATAATCAGACATTCGacaaattac                                                                                      <  2:3216981/150‑7 (MQ=255)

                                                               gTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAAc                                                                                      <  1:201625/150‑1 (MQ=255)

                                                                       caaacTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCga                                                                              >  2:3696827/1‑150 (MQ=255)

                                                                       caaacTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCga                                                                              >  2:1371245/1‑150 (MQ=255)

                                                                        aaacTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCgag                                                                             >  2:2060544/1‑150 (MQ=255)

                                                                         aacTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCgaga                                                                            >  1:1959094/1‑150 (MQ=255)

                                                                          acTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGaa                                                                           >  2:3833644/1‑150 (MQ=255)

                                                                               cAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACgg                                                                      <  1:1775986/150‑1 (MQ=255)

                                                                                aaaCATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGc                                                                     >  1:2833637/1‑150 (MQ=255)

                                                                                    aTTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGAc                                                                 >  2:3637922/1‑150 (MQ=255)

                                                                                      tGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCa                                                               >  1:4500578/1‑150 (MQ=255)

                                                                                         cacaAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGtt                                                            <  2:418722/150‑1 (MQ=255)

                                                                                          actaaGGCGGTCAGGGTCACCGTTTGCCCGGAAACCCCTTCCGAAATGCCTAAATCCTTGGCCTTTGGCGTCAACACACTGACGGGCAAAAAATCACCGATAATCAGACAGGCTACACAAAACGCCACCGAGAACACGGCTGACCAGTtc                                                           <  1:2489812/147‑1 (MQ=255)

                                                                                           caAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTtcg                                                          >  2:1503474/1‑150 (MQ=255)

                                                                                           caAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTtcg                                                          <  2:597565/150‑1 (MQ=255)

                                                                                             aaGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTtcggt                                                        >  1:846704/1‑150 (MQ=255)

                                                                                             aaGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTtcggt                                                        >  1:1783125/1‑150 (MQ=255)

                                                                                                  ggtcacggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGt                                                   >  1:3873436/1‑150 (MQ=255)

                                                                                                  ggtcacggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGt                                                   >  2:2643508/1‑150 (MQ=255)

                                                                                                   gtcacggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTg                                                  >  1:432049/1‑150 (MQ=255)

                                                                                                   gtcacggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTg                                                  <  1:1559894/150‑1 (MQ=255)

                                                                                                   gtcacggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTg                                                  >  1:649042/1‑150 (MQ=255)

                                                                                                      acggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATg                                               <  1:4432148/150‑1 (MQ=255)

                                                                                                      acggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATg                                               >  2:3750461/1‑150 (MQ=255)

                                                                                                       cggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATgg                                              >  1:2047821/1‑150 (MQ=255)

                                                                                                        ggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGc                                             >  1:4239256/1‑150 (MQ=255)

                                                                                                         gtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCa                                            >  1:692177/1‑150 (MQ=255)

                                                                                                            acCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATcc                                         <  2:572631/150‑1 (MQ=255)

                                                                                                               gATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccgcg                                      <  2:2449544/150‑1 (MQ=255)

                                                                                                                 ttGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccgcgcc                                    <  1:926142/150‑1 (MQ=255)

                                                                                                                  tGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccgcgccg                                   <  1:2966983/150‑1 (MQ=255)

                                                                                                                     ccGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCgg                                <  1:2155512/150‑1 (MQ=255)

                                                                                                                      cGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGt                               >  1:1017011/1‑150 (MQ=255)

                                                                                                                             cccTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGc                        <  1:409397/150‑1 (MQ=255)

                                                                                                                              ccTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCa                       <  1:3241486/150‑1 (MQ=255)

                                                                                                                                ttCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAAt                     >  2:1660323/1‑150 (MQ=255)

                                                                                                                                   cGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATaaa                  >  1:2434607/1‑150 (MQ=255)

                                                                                                                                     aaaTGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAAtt                <  1:4007286/150‑1 (MQ=255)

                                                                                                                                      aaTGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTc               <  1:2134801/150‑1 (MQ=255)

                                                                                                                                       aTGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCa              >  1:3544404/1‑150 (MQ=255)

                                                                                                                                        tGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCAc             <  1:3284668/150‑1 (MQ=255)

                                                                                                                                        tGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTGTCGGCAATAAATTCAc             >  2:337110/1‑150 (MQ=255)

                                                                                                                                            tAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCAt         <  1:1667657/150‑1 (MQ=255)

                                                                                                                                              aaTCCTTGGCCATTGGCGTCAACAAACTGAAGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATgg       >  1:4289307/1‑150 (MQ=255)

                                                                                                                                              aaTCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATgg       >  2:156597/1‑150 (MQ=255)

                                                                                                                                               aTCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGgt      <  1:2315267/150‑1 (MQ=255)

                                                                                                                                                   tGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTa  >  2:2085873/1‑150 (MQ=255)

AGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTA  >  CP009273/3834810‑3835106

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