BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,626,389	G→T	23.0%	I347I (ATC→ATA)	guaB ←	IMP dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,626,389	0	G	T	23.0%	1015.5 / 213.1	411	I347I (ATC→ATA)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (152/164); new base T (43/52); total (195/216)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 417 total aligned reads displayed.

CTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGCC  >  CP009273/2626241‑2626538

    ggACATCGCGCCCTGTTAACCAATACCACGGTAATTTTTGTTAGAACGGACCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTc                                                                                                                                                  <  1:4638546/150‑1 (MQ=255)

    ggACATCGCGCCCAGGGAACCCATACCACGTTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTc                                                                                                                                                  <  1:4508468/150‑1 (MQ=255)

    ggACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTc                                                                                                                                                  <  1:2774198/150‑1 (MQ=255)

    ggACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTc                                                                                                                                                  <  2:2101026/150‑1 (MQ=255)

       caccgcgCCCAGGGAACCCATCCCACGGAATGATTTGTAACAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGcc                                                                                                                                               <  1:2262844/147‑1 (MQ=255)

       cATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGcc                                                                                                                                               >  2:3863228/1‑150 (MQ=255)

       cATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGcc                                                                                                                                               <  2:1291094/150‑1 (MQ=255)

       cATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGcc                                                                                                                                               >  1:1896954/1‑150 (MQ=255)

       cATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGCCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGcc                                                                                                                                               >  2:1896641/1‑150 (MQ=255)

         tCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCgg                                                                                                                                             <  2:3326015/150‑1 (MQ=255)

              cccAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAg                                                                                                                                        <  1:704787/150‑1 (MQ=255)

               ccAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGc                                                                                                                                       <  1:998324/150‑1 (MQ=255)

               ccAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGc                                                                                                                                       >  2:683531/1‑150 (MQ=255)

                ctggtaaCCCATACCACGGTAAGTTTTGTACGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCg                                                                                                                                      <  1:2391395/145‑1 (MQ=255)

                cAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCg                                                                                                                                      >  1:5108085/1‑150 (MQ=255)

                cAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCg                                                                                                                                      <  2:832559/150‑1 (MQ=255)

                cAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCg                                                                                                                                      >  2:285251/1‑150 (MQ=255)

                 aGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGACGCGa                                                                                                                                     >  2:4796610/1‑150 (MQ=255)

                 aGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGa                                                                                                                                     >  1:953500/1‑150 (MQ=255)

                  gggAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGaa                                                                                                                                    >  2:1825971/1‑150 (MQ=255)

                   ggAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAAt                                                                                                                                   <  2:4418443/150‑1 (MQ=255)

                     aaCCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATAc                                                                                                                                 <  1:2017395/150‑1 (MQ=255)

                     aaCCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATAc                                                                                                                                 >  1:409437/1‑150 (MQ=255)

                      aCCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATAcc                                                                                                                                >  1:1670385/1‑150 (MQ=255)

                      aCCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATAcc                                                                                                                                >  1:1660735/1‑150 (MQ=255)

                       cccATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCg                                                                                                                               >  2:4177427/1‑150 (MQ=255)

                       cccATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCg                                                                                                                               >  1:2614788/1‑150 (MQ=255)

                              cACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAg                                                                                                                        >  1:1094563/1‑150 (MQ=255)

                               aCGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGc                                                                                                                       <  2:3259559/150‑1 (MQ=255)

                               aCGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGc                                                                                                                       >  1:2346332/1‑150 (MQ=255)

                                  gTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGAt                                                                                                                    <  2:677414/150‑1 (MQ=255)

                                   tAAGAATTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGAATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATATCTATGGCTATGTCGCCGGAGAAGCGAAAACCGCCATCAGCGATa                                                                                                                   >  1:3875650/1‑150 (MQ=255)

                                     aGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAAc                                                                                                                 <  2:3537244/150‑1 (MQ=255)

                                        tttGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCgg                                                                                                              <  2:1776197/150‑1 (MQ=255)

                                        tttGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCgg                                                                                                              >  2:3945851/1‑150 (MQ=255)

                                          tGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGaa                                                                                                            >  2:2370675/1‑150 (MQ=255)

                                           gTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAAt                                                                                                           <  2:37302/150‑1 (MQ=255)

                                           gTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAAt                                                                                                           >  1:566542/1‑150 (MQ=255)

                                              aGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATAcc                                                                                                        <  1:4740069/150‑1 (MQ=255)

                                              aGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATAcc                                                                                                        <  2:3458503/150‑1 (MQ=255)

                                              aGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATAcc                                                                                                        >  1:1130580/1‑150 (MQ=255)

                                               gAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCg                                                                                                       >  2:2779710/1‑150 (MQ=255)

                                                atcGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCgg                                                                                                      <  1:5015116/148‑1 (MQ=255)

                                                aaCGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCgg                                                                                                      <  1:2489240/150‑1 (MQ=255)

                                                  cGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTg                                                                                                    <  2:2048153/150‑1 (MQ=255)

                                                  cGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTg                                                                                                    <  1:1765536/150‑1 (MQ=255)

                                                  cGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTg                                                                                                    <  1:4435397/150‑1 (MQ=255)

                                                   ggCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGc                                                                                                   <  2:774325/150‑1 (MQ=255)

                                                    gCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGcc                                                                                                  >  1:3134663/1‑150 (MQ=255)

                                                     cccTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCt                                                                                                 <  2:3664231/150‑1 (MQ=255)

                                                          gTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAg                                                                                            >  2:4516958/1‑150 (MQ=255)

                                                           tAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAgg                                                                                           >  2:3368886/1‑150 (MQ=255)

                                                           tAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAgg                                                                                           <  1:1007394/150‑1 (MQ=255)

                                                               gTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCtt                                                                                       >  1:1470126/1‑150 (MQ=255)

                                                               gTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCtt                                                                                       <  1:3996296/150‑1 (MQ=255)

                                                                 tCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCt                                                                                     <  2:226058/150‑1 (MQ=255)

                                                                   gATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCAACAGCGATAACCGGAACACCGGTGCCTTCCAGGGCTTCTAc                                                                                   >  1:2707607/1‑150 (MQ=255)

                                                                    aTTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACt                                                                                  >  2:302013/1‑150 (MQ=255)

                                                                             ggAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCa                                                                         <  1:794647/150‑1 (MQ=255)

                                                                               agaTTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAAc                                                                       <  1:2045443/150‑1 (MQ=255)

                                                                                gaTTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACa                                                                      >  1:1056251/1‑150 (MQ=255)

                                                                                 aTTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAg                                                                     <  2:4020408/150‑1 (MQ=255)

                                                                                      tcAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTa                                                                <  1:1848322/150‑1 (MQ=255)

                                                                                       cAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTaa                                                               >  1:715698/1‑150 (MQ=255)

                                                                                              cGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCATTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGGGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCgg                                                        >  2:2607304/1‑150 (MQ=255)

                                                                                               gCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGa                                                       <  2:4987752/150‑1 (MQ=255)

                                                                                                 cAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAac                                                     >  2:4310459/1‑150 (MQ=255)

                                                                                                  aGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAaca                                                    <  2:1605884/150‑1 (MQ=255)

                                                                                                   gCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAacac                                                   <  1:1538415/150‑1 (MQ=255)

                                                                                                    cATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACAcc                                                  <  1:1282177/150‑1 (MQ=255)

                                                                                                    cATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACAcc                                                  <  1:3697305/150‑1 (MQ=255)

                                                                                                    cATGGAACCTACCATCACCGCGCGTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACAcc                                                  <  2:55144/150‑1 (MQ=255)

                                                                                                       ggAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGAc                                               >  2:115270/1‑150 (MQ=255)

                                                                                                            cTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAg                                          <  1:1606989/150‑1 (MQ=255)

                                                                                                               ccATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCa                                       <  2:2674245/150‑1 (MQ=255)

                                                                                                                 aTCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACg                                     <  2:5200236/150‑1 (MQ=255)

                                                                                                                 aTCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACg                                     <  2:1065929/150‑1 (MQ=255)

                                                                                                                   cACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGAt                                   >  2:1645015/1‑150 (MQ=255)

                                                                                                                    aCCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATa                                  >  2:1992458/1‑150 (MQ=255)

                                                                                                                          cTTGCGCCAGCGGCGATATCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGtt                            <  2:3222194/150‑1 (MQ=255)

                                                                                                                           ttGCGCCAGCGGCGATAGCTTTGGCTTTGTCGCCGGAGAAGCGAATACCGCCATCAGCGATATCCGGAATACCGGTGCCTTCCAGGGCTTCTTCTGCGTCAGCAACAGCGGTAATCTGCGGAACTCCGACGCCAGTCACGATACGAGTTg                           <  1:5319838/150‑1 (MQ=255)

                                                                                                                           ttGCGCCAGCGGCGATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTg                           <  1:5135669/150‑1 (MQ=255)

                                                                                                                                  aGCGGCGATAGCTTTGGCGATGTCTCCGGAGAAGCTAATACCGCCATCAGCGATAACCGGAATACCTGTGCCTTCCAGTTCTTCTACTGAGTCAGCAACAGCGTTAATCTGCGGAAAACCGAATCTAGTCACGATACGAGTTGTACAGAt                    <  1:3441408/150‑1 (MQ=255)

                                                                                                                                    cggcgATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATag                  <  1:656509/150‑1 (MQ=255)

                                                                                                                                    cggcgATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATag                  <  2:4726899/150‑1 (MQ=255)

                                                                                                                                    cggcgATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATag                  <  2:1907069/150‑1 (MQ=255)

                                                                                                                                     ggcgATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATaga                 <  2:1417584/150‑1 (MQ=255)

                                                                                                                                     ggcgATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATaga                 >  2:4612048/1‑150 (MQ=255)

                                                                                                                                       cgATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGc               >  1:2163070/1‑150 (MQ=255)

                                                                                                                                        gATAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGcc              >  2:4218246/1‑150 (MQ=255)

                                                                                                                                         aTAGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCa             <  1:2860746/150‑1 (MQ=255)

                                                                                                                                         aTAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCa             >  2:2912578/1‑150 (MQ=255)

                                                                                                                                           aGCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATCCGAGTTGTACAGATAGAGCCAgg           >  2:730551/1‑150 (MQ=255)

                                                                                                                                            gCTTTGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAggg          <  2:1684862/150‑1 (MQ=255)

                                                                                                                                            gCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAggg          <  2:4902623/150‑1 (MQ=255)

                                                                                                                                              tttGGCTATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGcc        >  2:2522227/1‑150 (MQ=255)

                                                                                                                                                tGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCaa      <  1:2373647/150‑1 (MQ=255)

                                                                                                                                                 ggCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAAt     <  1:3059087/150‑1 (MQ=255)

                                                                                                                                                  gCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATg    <  1:4177427/150‑1 (MQ=255)

                                                                                                                                                  gCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATg    <  2:1488567/150‑1 (MQ=255)

CTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGCC  >  CP009273/2626241‑2626538

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