BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,609,856	A→T	9.2%	H146L (CAC→CTC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,609,856	0	A	T	9.2%	515.5 / 33.7	171	H146L (CAC→CTC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand): ref base A (68/87); new base T (7/9); total (75/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

ATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGA  >  CP009273/4609709‑4610003

aTACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa                                                                                                                                                     >  2:2284435/1‑148 (MQ=255)

aTACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa                                                                                                                                                     >  2:2284580/1‑148 (MQ=255)

  aCCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCctct                                                                                                                                                   >  2:3175008/1‑148 (MQ=255)

   ccccGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCActc                                                                                                                                                  <  1:1506685/148‑1 (MQ=255)

   ccccGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCActc                                                                                                                                                  <  1:1758549/148‑1 (MQ=255)

    cccGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCActct                                                                                                                                                 <  1:2390867/148‑1 (MQ=255)

     ccGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCtt                                                                                                                                                <  1:2112854/148‑1 (MQ=255)

         cTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCgg                                                                                                                                            <  1:2734837/148‑1 (MQ=255)

              cACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACgg                                                                                                                                       >  2:954837/1‑148 (MQ=255)

                cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTc                                                                                                                                     >  2:241619/1‑148 (MQ=255)

                 tGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCa                                                                                                                                    <  1:2533660/148‑1 (MQ=255)

                  gggAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCAt                                                                                                                                   <  1:3268174/148‑1 (MQ=255)

                  gggAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCAt                                                                                                                                   <  2:2204899/148‑1 (MQ=255)

                   ggAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCAtt                                                                                                                                  <  1:2138331/148‑1 (MQ=255)

                            cGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCaa                                                                                                                         <  1:1326016/148‑1 (MQ=255)

                             ggTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAAc                                                                                                                        <  1:1639018/148‑1 (MQ=255)

                             ggTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAAc                                                                                                                        <  1:1639145/148‑1 (MQ=255)

                              gtgtCCCGGTTCTGTTTGAGTGGGTATATTTCTACGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACTGTCATTCTGGATCAACt                                                                                                                       <  2:3255561/148‑1 (MQ=255)

                              gtgtCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACt                                                                                                                       >  2:2760085/1‑148 (MQ=255)

                              gtgtCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACt                                                                                                                       <  2:3113852/148‑1 (MQ=255)

                              gtgtCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACt                                                                                                                       <  2:3113887/148‑1 (MQ=255)

                                 tCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTggg                                                                                                                    <  1:2587237/148‑1 (MQ=255)

                                      gTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAag                                                                                                               <  2:2387532/148‑1 (MQ=255)

                                       ttCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAaga                                                                                                              <  1:2592525/148‑1 (MQ=255)

                                           gTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGcac                                                                                                          >  2:2025200/1‑148 (MQ=255)

                                            tttGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGcaca                                                                                                         >  2:1647721/1‑148 (MQ=255)

                                                aGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGaa                                                                                                     >  1:475559/1‑148 (MQ=255)

                                                 gTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGaaa                                                                                                    <  1:1396154/148‑1 (MQ=255)

                                                   ggggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAAc                                                                                                  >  1:200892/1‑148 (MQ=255)

                                                   ggggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAAc                                                                                                  <  1:2696046/148‑1 (MQ=255)

                                                   ggggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAAc                                                                                                  <  1:2696184/148‑1 (MQ=255)

                                                     ggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCg                                                                                                <  1:1508802/148‑1 (MQ=255)

                                                          tttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAg                                                                                           >  2:3109630/1‑148 (MQ=255)

                                                           ttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCTCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGc                                                                                          <  1:3175008/148‑1 (MQ=255)

                                                           ttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGc                                                                                          <  2:2606411/148‑1 (MQ=255)

                                                           ttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGc                                                                                          <  2:1473640/148‑1 (MQ=255)

                                                           ttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGc                                                                                          <  1:2870588/148‑1 (MQ=255)

                                                                 gATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCTCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAttt                                                                                    <  1:622304/148‑1 (MQ=255)

                                                                  aTCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAtttt                                                                                   <  2:2705143/148‑1 (MQ=255)

                                                                  aTCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAtttt                                                                                   <  1:101714/148‑1 (MQ=255)

                                                                   tCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                                                                                  >  1:2861596/1‑148 (MQ=255)

                                                                   tCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                                                                                  >  1:2790969/1‑148 (MQ=255)

                                                                   tCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                                                                                  <  2:692816/148‑1 (MQ=255)

                                                                     aCGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCta                                                                                >  1:3271627/1‑148 (MQ=255)

                                                                      cGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCtat                                                                               <  2:576278/148‑1 (MQ=255)

                                                                         aaaCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATAcc                                                                            <  2:1000096/148‑1 (MQ=255)

                                                                         aaaCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATAcc                                                                            >  2:752457/1‑148 (MQ=255)

                                                                          aaCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCt                                                                           >  1:2145214/1‑148 (MQ=255)

                                                                           aCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCTCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTc                                                                          >  1:551911/1‑148 (MQ=255)

                                                                           aCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCTCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTc                                                                          >  1:2212173/1‑148 (MQ=255)

                                                                           aCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTc                                                                          >  1:1632320/1‑148 (MQ=255)

                                                                           aCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTc                                                                          >  1:1627073/1‑148 (MQ=255)

                                                                           aCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTc                                                                          <  2:208989/148‑1 (MQ=255)

                                                                             aGCTTCCCGCAAGAGCTGCTGTATAAACTGGTCGAACGCGATAATCTGCCGGGTTACCTCGTTAACTGCCTCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCg                                                                        <  1:515627/148‑1 (MQ=255)

                                                                             aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCg                                                                        <  1:26233/148‑1 (MQ=255)

                                                                               cTTCCCTCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCt                                                                      <  2:1557857/148‑1 (MQ=255)

                                                                               cTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCt                                                                      >  1:2051763/1‑148 (MQ=255)

                                                                               cTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACAAGAAAACCGGCAAGCCGATTTTCTATACCTCCGCt                                                                      >  1:918495/1‑148 (MQ=255)

                                                                                 tCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGa                                                                    >  1:539727/1‑148 (MQ=255)

                                                                                  cccGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCTCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGAc                                                                   <  1:2818660/148‑1 (MQ=255)

                                                                                     gCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTcc                                                                >  1:2364519/1‑148 (MQ=255)

                                                                                      cAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCg                                                               <  2:403536/148‑1 (MQ=255)

                                                                                        agagCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCgtg                                                             >  2:2480054/1‑148 (MQ=255)

                                                                                        agagCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCgtg                                                             >  2:2480011/1‑148 (MQ=255)

                                                                                        agagCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCgtg                                                             <  1:553218/148‑1 (MQ=255)

                                                                                              gctgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAg                                                       >  1:1456282/1‑148 (MQ=255)

                                                                                               ctgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCTCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGa                                                      <  2:3051820/148‑1 (MQ=255)

                                                                                               ctgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGa                                                      <  2:876837/148‑1 (MQ=255)

                                                                                                  gATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtg                                                   <  1:800614/148‑1 (MQ=255)

                                                                                                  gATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtg                                                   >  1:1137016/1‑148 (MQ=255)

                                                                                                   aTAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgc                                                  >  1:2247009/1‑148 (MQ=255)

                                                                                                      aaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcctg                                               >  2:973155/1‑148 (MQ=255)

                                                                                                      aaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcctg                                               <  1:1291408/148‑1 (MQ=255)

                                                                                                      aaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcctg                                               <  1:669061/148‑1 (MQ=255)

                                                                                                      aaCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcctg                                               <  1:668998/148‑1 (MQ=255)

                                                                                                        cttgTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcctgcc                                             >  1:1748719/4‑148 (MQ=255)

                                                                                                            tCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATga                                         >  1:1763516/1‑148 (MQ=255)

                                                                                                             cGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATgaa                                        <  1:1031481/148‑1 (MQ=255)

                                                                                                              gAACGCGCTAATCGGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATgaag                                       <  2:1600383/148‑1 (MQ=255)

                                                                                                               aaCGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATgaaga                                      <  2:374029/148‑1 (MQ=255)

                                                                                                                aCGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATgaagaa                                     >  2:1939271/1‑148 (MQ=255)

                                                                                                                 cgcgCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGaaa                                    >  2:14522/1‑148 (MQ=255)

                                                                                                                 cgcgCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCAGCCATGACGaaa                                    >  2:14573/1‑148 (MQ=255)

                                                                                                                   cgcTAATCTGCCGGGTTACCTCGGTAACTGCCTCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACt                                  >  1:1781887/1‑148 (MQ=255)

                                                                                                                   cgcTAATCTGCCGGGTTACCTCGGTAACTGCCTCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACt                                  >  1:2557290/1‑148 (MQ=255)

                                                                                                                   cgcTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACt                                  <  1:2873047/148‑1 (MQ=255)

                                                                                                                       aaTCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCg                              <  2:2144042/148‑1 (MQ=255)

                                                                                                                          cTGCCGTGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCTTCTGGGCGAAGATCAGATGATAAGCGGCAATCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGTTTGCCTGCCATGTAGAAACTTTAGGTc                           <  2:1639782/148‑1 (MQ=255)

                                                                                                                          cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTc                           <  1:3109630/148‑1 (MQ=255)

                                                                                                                          cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTc                           <  2:140971/148‑1 (MQ=255)

                                                                                                                           tGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCt                          <  1:2929403/148‑1 (MQ=255)

                                                                                                                                 gTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATaa                    <  1:752457/148‑1 (MQ=255)

                                                                                                                                  ttACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATaaa                   <  1:2284580/148‑1 (MQ=255)

                                                                                                                                     ccTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAActc                <  2:266060/148‑1 (MQ=255)

                                                                                                                                      caccgTAACTGCCTCTCTTCCGGTACGGTCATTCTGGATCAACTCCGCGAATAGAACATGAAAACAGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATCAAGAAACTTTCGGTCTGGATAAActct               >  2:16366/5‑148 (MQ=255)

                                                                                                                                          gTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGa           >  1:1369671/1‑148 (MQ=255)

                                                                                                                                           tAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGaa          >  1:9492/1‑148 (MQ=255)

                                                                                                                                             aCTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACt        >  2:2038480/1‑148 (MQ=255)

                                                                                                                                              cTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACtg       >  1:2784664/1‑148 (MQ=255)

                                                                                                                                                   aCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGa  >  2:3302137/1‑148 (MQ=255)

ATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGA  >  CP009273/4609709‑4610003

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