BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,649,460	T→A	49.8%	K65* (AAA→TAA)	hdeB ←	acid‑resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,649,460	0	T	A	49.8%	‑3.3 / inf	251	K65* (AAA→TAA)	hdeB	acid‑resistance protein
Reads supporting (aligned to +/- strand): ref base T (62/64); new base A (69/56); total (131/120)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 260 total aligned reads displayed.

ACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGCCGCTA  >  CP009273/3649314‑3649606

aCTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTta                                                                                                                                                   >  2:9265/1‑148 (MQ=255)

          cATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGACGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGttt                                                                                                                                         >  1:619358/1‑148 (MQ=255)

          cATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGttt                                                                                                                                         <  1:603632/148‑1 (MQ=255)

          cATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGttt                                                                                                                                         >  1:1297807/1‑148 (MQ=255)

           aTTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTttc                                                                                                                                        <  2:1729585/148‑1 (MQ=255)

            ttAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTttct                                                                                                                                       >  1:1674073/1‑148 (MQ=255)

              aaTTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTttcttc                                                                                                                                     <  2:3038039/148‑1 (MQ=255)

               aTTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCa                                                                                                                                    <  1:1419153/148‑1 (MQ=255)

               aTTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCa                                                                                                                                    >  2:1338289/1‑148 (MQ=255)

                        gTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATc                                                                                                                           <  2:769110/148‑1 (MQ=255)

                           aTTATATGATTGATTTTTGAAGGTATACATATTTTTCTGCTGTTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTATAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATccac                                                                                                                        <  2:3141519/148‑1 (MQ=255)

                            ttAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATccacc                                                                                                                       <  1:96533/148‑1 (MQ=255)

                            ttAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATccacc                                                                                                                       >  1:1766322/1‑148 (MQ=255)

                               gATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATg                                                                                                                    >  1:1386390/1‑148 (MQ=255)

                               gATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATg                                                                                                                    >  2:1918964/1‑148 (MQ=255)

                               gATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATg                                                                                                                    <  2:2883070/148‑1 (MQ=255)

                                  gCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCaa                                                                                                                 >  1:1072662/1‑148 (MQ=255)

                                      gATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATccacc                                                                                                                       <  2:1304982/138‑1 (MQ=255)

                                      gATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCgg                                                                                                             >  1:2256549/1‑148 (MQ=255)

                                         ttttGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGt                                                                                                          <  2:1753091/148‑1 (MQ=255)

                                         ttttGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGt                                                                                                          <  1:2763140/148‑1 (MQ=255)

                                               aGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGc                                                                                                    <  1:190785/148‑1 (MQ=255)

                                               aGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGc                                                                                                    <  1:3153863/148‑1 (MQ=255)

                                                ggTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCt                                                                                                   <  2:3240662/148‑1 (MQ=255)

                                                 gTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCtt                                                                                                  <  1:3121108/148‑1 (MQ=255)

                                                 gTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCtt                                                                                                  <  1:2822964/148‑1 (MQ=255)

                                                 gTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCtt                                                                                                  >  2:375559/1‑148 (MQ=255)

                                                  tataCAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCttt                                                                                                 >  1:1187483/1‑148 (MQ=255)

                                                  tataCAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCttt                                                                                                 >  1:2858167/1‑148 (MQ=255)

                                                         aTTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTc                                                                                          >  2:1875034/1‑148 (MQ=255)

                                                          tttgtCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCa                                                                                         <  1:1137052/148‑1 (MQ=255)

                                                           ttttCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAg                                                                                        >  1:2687470/1‑148 (MQ=255)

                                                             ttCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGAt                                                                                      >  2:2648354/1‑148 (MQ=255)

                                                                tGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCaa                                                                                   <  2:1539573/148‑1 (MQ=255)

                                                                 gCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAAt                                                                                  <  2:679683/148‑1 (MQ=255)

                                                                       ttttCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTc                                                                            <  2:899750/148‑1 (MQ=255)

                                                                       ttttCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTc                                                                            <  1:1270571/148‑1 (MQ=255)

                                                                       ttttCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTc                                                                            >  1:1746368/1‑148 (MQ=255)

                                                                       ttttCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTc                                                                            >  1:1746427/1‑148 (MQ=255)

                                                                       ttttCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTc                                                                            <  1:199936/148‑1 (MQ=255)

                                                                           cTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTgg                                                                        <  1:2971890/148‑1 (MQ=255)

                                                                             tACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCa                                                                      <  1:2164074/148‑1 (MQ=255)

                                                                             tACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCa                                                                      >  2:491808/1‑148 (MQ=255)

                                                                             tACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCa                                                                      >  1:473474/1‑148 (MQ=255)

                                                                              aCAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAg                                                                     >  2:3545693/1‑148 (MQ=255)

                                                                                 gTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTc                                                                  >  1:2927077/1‑148 (MQ=255)

                                                                                 gTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTc                                                                  <  2:1249271/148‑1 (MQ=255)

                                                                                  tATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCa                                                                 >  1:1041505/1‑148 (MQ=255)

                                                                                   aTTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCat                                                                >  2:2859409/1‑148 (MQ=255)

                                                                                   aTTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCat                                                                >  2:2649900/1‑148 (MQ=255)

                                                                                    ttCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCata                                                               >  2:583561/1‑148 (MQ=255)

                                                                                     tCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCACCCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCatat                                                              >  1:1831651/1‑148 (MQ=255)

                                                                                     tCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCatat                                                              <  2:3238588/148‑1 (MQ=255)

                                                                                     tCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCatat                                                              <  1:2372540/148‑1 (MQ=255)

                                                                                        aTCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCtt                                                           >  1:1934307/1‑148 (MQ=255)

                                                                                        aTCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCtt                                                           >  1:3347535/1‑148 (MQ=255)

                                                                                         tCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCttt                                                          >  2:3523070/1‑148 (MQ=255)

                                                                                          cACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTa                                                         <  2:545179/148‑1 (MQ=255)

                                                                                           aCTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAg                                                        <  1:409633/148‑1 (MQ=255)

                                                                                              ttAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCgg                                                     <  1:2982038/148‑1 (MQ=255)

                                                                                              ttAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCgg                                                     <  2:655524/148‑1 (MQ=255)

                                                                                               tAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGa                                                    <  2:246447/148‑1 (MQ=255)

                                                                                                 ggAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCATTAAATTCCTGGCAGGTGATCTCTTTAGCGGAtt                                                  >  2:1821905/1‑148 (MQ=255)

                                                                                                 ggAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGAtt                                                  <  2:1837033/148‑1 (MQ=255)

                                                                                                  gAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTc                                                 <  1:853066/148‑1 (MQ=255)

                                                                                                  gAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTc                                                 <  1:346187/148‑1 (MQ=255)

                                                                                                  gAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTc                                                 >  2:2655027/1‑148 (MQ=255)

                                                                                                  gAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTc                                                 >  1:2467658/1‑148 (MQ=255)

                                                                                                  gAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGGATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTc                                                 <  1:1707683/148‑1 (MQ=255)

                                                                                                     tttGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCAtt                                              >  1:2308125/1‑148 (MQ=255)

                                                                                                       tgagtgagATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTgg                                            <  1:1548297/148‑1 (MQ=255)

                                                                                                       tgagtgagATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTgg                                            <  2:2568397/148‑1 (MQ=255)

                                                                                                            gagaTCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGcc                                       >  2:2748681/1‑148 (MQ=255)

                                                                                                             agaTCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCa                                      >  2:504888/1‑148 (MQ=255)

                                                                                                             agaTCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCa                                      >  1:2137176/1‑148 (MQ=255)

                                                                                                              gaTCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCaa                                     >  2:2315500/1‑148 (MQ=255)

                                                                                                                tCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACg                                   <  1:2353232/148‑1 (MQ=255)

                                                                                                                    tttCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGa                               <  2:1011705/148‑1 (MQ=255)

                                                                                                                    tttCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGa                               <  2:1155702/148‑1 (MQ=255)

                                                                                                                       cATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATtg                            >  2:845346/1‑148 (MQ=255)

                                                                                                                        aTTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATtgt                           >  2:3436474/1‑148 (MQ=255)

                                                                                                                          ttAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGc                         >  1:1671234/1‑148 (MQ=255)

                                                                                                                          ttAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGc                         >  1:789839/1‑148 (MQ=255)

                                                                                                                            aaaTTAACGTTATCGCCAACTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGt                       <  2:2146060/148‑1 (MQ=255)

                                                                                                                             aattAAAGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTTGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCATATTGTGCGtt                      <  2:2248249/145‑1 (MQ=255)

                                                                                                                               gTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCa                    >  1:3562370/1‑148 (MQ=255)

                                                                                                                                  aCGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCa                 <  1:3376459/148‑1 (MQ=255)

                                                                                                                                   cGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAg                <  2:2728490/148‑1 (MQ=255)

                                                                                                                                     gTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTg              <  1:9265/148‑1 (MQ=255)

                                                                                                                                       aTCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGAc            >  1:1637171/1‑148 (MQ=255)

                                                                                                                                       aTCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGAc            >  1:1637280/1‑148 (MQ=255)

                                                                                                                                        tCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACa           <  2:2210833/148‑1 (MQ=255)

                                                                                                                                         cGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACaa          >  2:2528911/1‑148 (MQ=255)

                                                                                                                                           ccACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAg        <  1:529258/148‑1 (MQ=255)

                                                                                                                                            cACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGc       >  2:2029797/1‑148 (MQ=255)

                                                                                                                                              ccTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGTCg     <  1:889553/148‑1 (MQ=255)

                                                                                                                                              ccTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGCCg     >  2:2614228/1‑148 (MQ=255)

                                                                                                                                               cTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGCCGc    <  1:2402282/148‑1 (MQ=255)

                                                                                                                                               cTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGCCGc    <  1:1221870/148‑1 (MQ=255)

                                                                                                                                                tttATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGCCGCt   >  1:3318014/1‑148 (MQ=255)

ACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGCCGCTA  >  CP009273/3649314‑3649606

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