BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,559,239	T→A	7.2%	S235C (AGT→TGT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,559,239	0	T	A	7.2%	445.1 / 2.5	179	S235C (AGT→TGT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base T (71/89); new base A (6/9); total (79/100)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

GCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAGCGCGT  >  CP009273/3559092‑3559386

gCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGTAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCACt                                                                                                                                                     >  2:1757158/1‑148 (MQ=255)

         gCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACg                                                                                                                                            >  2:1318349/1‑148 (MQ=255)

            cACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAg                                                                                                                                         <  2:3200968/148‑1 (MQ=255)

             aCGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTCTTCCACAGCCCCGAAGTAGTCACCCGGGTTGAATTTAACGAGGTTAATTTAGCTACTGGCGTGCGCACTCCCCAAACGCAGc                                                                                                                                        >  1:213389/1‑148 (MQ=255)

             aCGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACAGCCGCCAAGTAGTCACCCTGGTTGAATTTAACGAGGGTAATTTAGCTAAAGGCTTTCGCACTCCACAAACCCAGc                                                                                                                                        >  1:1412968/1‑148 (MQ=255)

               gCCCGGAGTAGGTGGAGTCATCCGGATACCGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCgt                                                                                                                                      <  1:2204062/148‑1 (MQ=255)

                cccGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCgtg                                                                                                                                     <  2:1714054/148‑1 (MQ=255)

                    gAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTgg                                                                                                                                 <  1:79434/148‑1 (MQ=255)

                       tAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTcg                                                                                                                              <  1:763753/148‑1 (MQ=255)

                          gTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACTGCCGCGAAGTAGTCACCCGGGGTTAATTTAACGAGGTTTATTTCGCTACTGGCTTGCGCCCTCCACCAACGCCGCGGGTTGGGCGCGt                                                                                                                           >  1:3058088/1‑148 (MQ=255)

                            ggAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGGTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCg                                                                                                                         >  2:1632735/1‑148 (MQ=255)

                             gAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACACTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCgg                                                                                                                        <  1:2821986/148‑1 (MQ=255)

                               gTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgtg                                                                                                                      <  2:568691/148‑1 (MQ=255)

                               gTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgtg                                                                                                                      <  1:113765/148‑1 (MQ=255)

                                  aTCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTGGTCACCCGGGTTGTATTTACCCAGGTTTATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGGGTTGGGCGCGTCGgtgtgg                                                                                                                   >  1:2040727/1‑146 (MQ=255)

                                  aTCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACAGCCGCGAAGTATTCCCCCGGGTTGAATTTACCGAGGGTAAGTTGGCTACAGGGTTGCGCCCGCCACAAAAGCAGCGGGGTGGTCGCGTGGGGGTCg                                                                                                                   >  1:2000068/1‑148 (MQ=255)

                                    ccGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTa                                                                                                                 <  1:1042250/148‑1 (MQ=255)

                                    ccGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGGTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGGGTTGGTCGCGTCGGTGTCGTa                                                                                                                 >  2:2317743/1‑148 (MQ=255)

                                      ggATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAAc                                                                                                               <  1:263596/148‑1 (MQ=255)

                                       gATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAAcc                                                                                                              <  2:1452779/148‑1 (MQ=255)

                                       gATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAAcc                                                                                                              <  2:1451596/148‑1 (MQ=255)

                                        aTACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCa                                                                                                             <  2:3348116/148‑1 (MQ=255)

                                        aTACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGCGGTTTATTTCGCTACAGGCGTGCGCACTCCACCAACGCCGCGGGTTGGTCGCGTCGGTGTCGTAAccc                                                                                                             >  2:1505551/1‑147 (MQ=255)

                                        aTACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTATTTTCGCTACTGGCTTTCGCCCTCCACAAACGCAGCGTGTTGGTCGCGTCGGGGTCGTAAACa                                                                                                             >  2:990000/1‑148 (MQ=255)

                                           cAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAggg                                                                                                          >  2:293557/1‑148 (MQ=255)

                                           cAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAAGTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAggg                                                                                                          <  2:3432062/148‑1 (MQ=255)

                                           cAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGGAGTCACCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAggg                                                                                                          >  2:3083189/1‑148 (MQ=255)

                                           cAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAAACAggg                                                                                                          >  2:293512/1‑148 (MQ=255)

                                           cAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACTGCCGCGACGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGGCGCGTCGGTGTCGTAACCAggg                                                                                                          >  2:2278781/1‑148 (MQ=255)

                                            aGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGa                                                                                                         <  1:1248465/148‑1 (MQ=255)

                                              tACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGAtt                                                                                                       <  1:1482355/148‑1 (MQ=255)

                                                cgcgAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTAt                                                                                                     <  2:1768326/148‑1 (MQ=255)

                                                  cgAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCCGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCt                                                                                                   >  2:1076239/1‑148 (MQ=255)

                                                   gagaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTg                                                                                                  <  1:433163/148‑1 (MQ=255)

                                                   gagaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTg                                                                                                  <  1:2454876/148‑1 (MQ=255)

                                                    agaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGa                                                                                                 <  2:3286111/148‑1 (MQ=255)

                                                    agaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGa                                                                                                 <  1:2111161/148‑1 (MQ=255)

                                                    agaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGa                                                                                                 <  1:1573203/148‑1 (MQ=255)

                                                     gaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGAt                                                                                                <  2:810888/148‑1 (MQ=255)

                                                      aTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATc                                                                                               <  2:2687432/148‑1 (MQ=255)

                                                        aCGTTCTCGGAGTGGTTTTTATCTTACACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTGGTAACCACGGATTATATGATCGt                                                                                             >  2:3274088/1‑148 (MQ=255)

                                                        aCGTTCTCGGAGTGGTTTTTATCTTACACAGCCCCGAAGTAGTCACCCGGGTTGAATTTAACGAGGTTAATTTCGCTACAGGCTTGCGCACTCCACAAAAGCAGCGGGTTGGTCGCGTCGGGGTCGTAAACAGGGATTATATGATCGt                                                                                             >  2:632062/1‑148 (MQ=255)

                                                         cGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTa                                                                                            >  2:3211077/1‑148 (MQ=255)

                                                         cGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTa                                                                                            <  1:642611/148‑1 (MQ=255)

                                                           ttCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAg                                                                                          <  2:3168271/148‑1 (MQ=255)

                                                             ctcGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCg                                                                                        <  2:3414453/148‑1 (MQ=255)

                                                               gggtgTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTTAATTTACCGAGGTTATTTTCGCTAGTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGTATTATCTGATCGTAAGCGAc                                                                                      <  2:175787/144‑1 (MQ=255)

                                                                  aGTGGTTTTTATATTCCACATCCCCCAAGTAGTCACCCAGGTTGAATTTAACGAGGTTAATTTAGCTAATGGCTTTCGCACACCACAAAACCAGCGTGTTGGTCGCGTCGGGGTAGTAAACAGGGATTATATTATCATAACAGACtac                                                                                   >  2:3265448/1‑146 (MQ=255)

                                                                   gtgtTTTTTATCTTCCACTGCCGCTAATTAGTAACCGTTTTTTATTTTACCGTGTTTATTTTCGCTACTGGCTTGCGCACTCAACAAACGCAGGGTTTTGGTCGCTTGGTTTTCGTAACCAGGTATTATCTGATCGTAAGCGACTccc                                                                                  <  2:2786602/148‑1 (MQ=255)

                                                                   gTGGTTTTTAGCTTCCACTGCCGCTAATTAGTCACCCTGGTTGATTTTACCGAGGTTAATTTCGCTACTGGGTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCTTAACCAGGGTTTATCTGATCGTAAGCGACTccc                                                                                  <  2:526945/148‑1 (MQ=255)

                                                                       tttttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGaa                                                                              >  1:3008799/1‑148 (MQ=255)

                                                                           tATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtctc                                                                          >  2:481205/1‑148 (MQ=255)

                                                                           tATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtctc                                                                          <  2:2621434/148‑1 (MQ=255)

                                                                             tCTTCCACTGCCGCGATTTTTTAACCATGTTTGAATTTACCGAGGTTAATTTCGCTACTGTGTTGCGCAATCCAAAAACGATGGTTGTTGGTCGCGTCTGTTTCGTAACCTGTTTTTATCTTATCGTAAGCGACTCCCAGAATCTCtt                                                                        <  2:2034642/148‑1 (MQ=255)

                                                                               ttcccgtGCCGCGATGTTGTCACCGGGGTTGAATTTACCGTGGTTATTTGCGCTACGGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCg                                                                      <  2:3015287/142‑1 (MQ=255)

                                                                                   aCTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGTTTTTCTGATCGTAAGCGACTCCCAGAATCTCTTCGGttt                                                                  <  2:1434168/148‑1 (MQ=255)

                                                                                    cTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTc                                                                 >  1:2888997/1‑148 (MQ=255)

                                                                                    cTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTc                                                                 >  2:457795/1‑148 (MQ=255)

                                                                                          cgAAGTAGTCACCCTGGTTTATTTTACCGAGGTTATTTTCGCTACTGGCTTGCGCACTCCACAAACGCCGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGTATTATCTTATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCa                                                           <  1:2921952/148‑1 (MQ=255)

                                                                                               tagtCACCCTGTTTGAATTTACCGAGGTTAATTTCGCTAGGGGCTTGCGCACACCACAAACGCAGCGTTTTGGTCGCGTCGGTGTCGTAACCAGGTATTATCTGATCTTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGcgcg                                                      <  1:840950/148‑1 (MQ=255)

                                                                                               tagtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGcgcg                                                      >  1:335909/1‑148 (MQ=255)

                                                                                                agtCACCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCGGGCGCACACCACAAACGCAGCGTGTTGGGCGGGTCGGTGTCGTAACCAGGGATTATCTGATCGTACGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGt                                                     <  2:347964/148‑1 (MQ=255)

                                                                                                 gtCACCCTGGTTTAATTTACCGGGGTTAATTTCGCTACTGGCTGGCGCACTCCACAAACGCCGCGTGTTGGTCGCGGCGGTGTCTTAACCAGGTATTATCTGATCGTACGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtt                                                    <  2:492266/148‑1 (MQ=255)

                                                                                                   aaCCCTGTTTGATTTTACCGTGGTTAATTTCGATACTGGCTTGGGCACTCCACAAACGCAGCGTGTTGGACGCGTGGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtttt                                                  <  2:181885/147‑1 (MQ=255)

                                                                                                    ccccTGTTTGAATTTACCGTGTTTATTTTCGCTACTGGGTTGCCCACACAAAAAACGCAGGGTGTTGGTCGCGTGGGTGTCTTAACCAGGTTTTATCTGATCGAAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGttttt                                                 <  1:2461571/147‑1 (MQ=255)

                                                                                                    aCCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGttttt                                                 >  2:2452770/1‑148 (MQ=255)

                                                                                                         ggTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGAAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTAcc                                            <  2:2883139/148‑1 (MQ=255)

                                                                                                              attttACCGTGGTTAATTTCGCTCCTGGCTTGCGCACTCCACAAACGCAGGGTGTGGGTCGCGTCGGTGTCGTAACCAGGGATTATCTTATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCct                                       <  1:2204926/146‑1 (MQ=255)

                                                                                                               ttttACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAACCGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctg                                      <  2:1601132/147‑1 (MQ=255)

                                                                                                               ttttACCGAGGTTAATTTCGCTACGGGCTGGCGCACTCCACAAACGCAGCGTGTTGGACGCGGCGGTGTCGTAACCAGGTATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctg                                      <  2:639399/147‑1 (MQ=255)

                                                                                                               aTTTACCGAGTTTTATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctg                                      <  2:2815590/148‑1 (MQ=255)

                                                                                                               aTTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctg                                      >  1:133297/1‑148 (MQ=255)

                                                                                                                tttACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctgc                                     >  2:261428/1‑148 (MQ=255)

                                                                                                                tttACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctgc                                     >  1:3377781/1‑148 (MQ=255)

                                                                                                                   aCCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGa                                  >  1:2592143/1‑148 (MQ=255)

                                                                                                                   aCCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGa                                  >  1:2764975/1‑148 (MQ=255)

                                                                                                                   aCCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGa                                  >  2:242533/1‑148 (MQ=255)

                                                                                                                     cGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAAt                                >  1:313784/1‑148 (MQ=255)

                                                                                                                       aGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGc                              >  1:3363665/1‑148 (MQ=255)

                                                                                                                        ggTTAATTTCGCTACTGGCTTGCGCACTCCACACACGCAGCGGGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCg                             >  2:1660742/1‑148 (MQ=255)

                                                                                                                        ggTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCg                             >  1:1697114/1‑148 (MQ=255)

                                                                                                                         gTTAATTTCGCTACTGGCTGGCGCACCCCACAACCGCAGCGTGTTGGGCGCGTCGGTGTCGTAACCAGGTATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGa                            <  2:3363665/148‑1 (MQ=255)

                                                                                                                                     agtGGCTTGCGCACTCCACAAACGCAGGGTGTTGGTGGCGTGGGTGTCTTAACCTGGTTTTATCTGATCTTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACgg                <  1:92227/146‑1 (MQ=255)

                                                                                                                                     aCTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACgg                >  1:29077/1‑148 (MQ=255)

                                                                                                                                        gggttGCCCAATCCAAAAACGCTGGGTGTTGGTCGCGTGGGTGTCTTACCCTTGTTTTATATGATCGAACGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACt             <  1:902879/145‑1 (MQ=255)

                                                                                                                                        ggCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACt             >  1:3259744/1‑148 (MQ=255)

                                                                                                                                        ggCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACt             <  2:1411705/148‑1 (MQ=255)

                                                                                                                                        ggCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACt             <  2:213277/148‑1 (MQ=255)

                                                                                                                                        ggCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACt             <  2:213389/148‑1 (MQ=255)

                                                                                                                                         gCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACtt            >  2:2678553/1‑148 (MQ=255)

                                                                                                                                          cTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACttt           >  2:754462/1‑148 (MQ=255)

                                                                                                                                           ttGCGCACTCCACAACCGCAGGGTGTTGGTCGCGTCGGTGTCGTAACCTGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTta          <  1:2521700/148‑1 (MQ=255)

                                                                                                                                           ttGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTta          <  2:3377781/148‑1 (MQ=255)

                                                                                                                                           ttGCGCACTCAACAACCGCTGGGTGTTGGTCGGGTGGTTGTGTTAACCTGGTATTATATGATCTTAAGCGACTCCCAGAATCTCGTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTta          <  1:994232/148‑1 (MQ=255)

                                                                                                                                            tGCGCACACCAAAACCGCAGGTTGTGGGTGGGGTGGTTTTCTTAACCTGGTTTTTTGTGATCGAAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTtat         <  1:68348/148‑1 (MQ=255)

                                                                                                                                              cgcACTCCACAAACGCAGGTTGTTGGGCGCGTCGGTGTGGTAACCTGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAg       <  2:174441/148‑1 (MQ=255)

                                                                                                                                              cgcACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAg       <  2:3460664/148‑1 (MQ=255)

                                                                                                                                                  cTCCAAAAACGCCGGTTGTTGGTCGGGTCGGTGTCGTAACCTGGGTTTTTCTGATCTTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAGcgcg   <  1:1757305/148‑1 (MQ=255)

                                                                                                                                                   tCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAGCGCgt  >  1:1355284/1‑148 (MQ=255)

                                                                                                                                                   tCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAGCGCgt  >  1:440276/1‑148 (MQ=255)

GCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAGCGCGT  >  CP009273/3559092‑3559386

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