BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,559,239	T→A	8.7%	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	8.7%	693.8 / 4.5	242	S235C (AGT→TGT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base T (108/104); new base A (20/6); total (128/114)
Fisher's exact test for biased strand distribution p-value = 1.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

ACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAGCGCGTGT  >  CP009273/3559090‑3559388

aCGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACt                                                                                                                                                       <  1:3477195/150‑1 (MQ=255)

 cGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGACGTAGTCACCCTGGTTGAATTTACCGAGGTTACTGTCGCTACTGGCGTGCGCCCTc                                                                                                                                                      >  2:2632698/1‑150 (MQ=255)

 cGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTTTTTTCCACATCCCCGCAGTAGTGACCCAGGTTTAATTTAAAGAGGGTAATTTAGCTACAGGCGTGCGCACTc                                                                                                                                                      >  1:65907/1‑150 (MQ=255)

   cAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACAGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTAGCTACTGGCTTGCGCACTCca                                                                                                                                                    >  1:1554695/1‑150 (MQ=255)

        gcAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTATTCCACAGGCCCGAACTAGTAACCCAGGTTTAATTTTAACAGGTTAATTTAGCAACTGGCATGAGCACACCACaaaa                                                                                                                                               >  1:1105690/1‑149 (MQ=255)

        gcAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCCCCCTGGTTGAATTTACCGCGGGTAATTTCGCTACTGGCTTGCGCACTCcacccac                                                                                                                                               >  2:4512621/1‑146 (MQ=255)

          aGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGc                                                                                                                                             >  2:2679578/1‑150 (MQ=255)

          aGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGc                                                                                                                                             >  1:4606114/1‑150 (MQ=255)

          aGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACAGCCGCGAAGTAGTACACCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGc                                                                                                                                             >  2:1712317/1‑150 (MQ=255)

          aGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGGTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGc                                                                                                                                             >  2:2536878/1‑150 (MQ=255)

           gCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCa                                                                                                                                            <  1:2098193/150‑1 (MQ=255)

            cTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTCCACACCCCCGCCCTGGTCACCCGGGTTGAAATTAACCAGGTTAATTTAGCTACAGGCTTGGGCACACCACAAAAGCac                                                                                                                                           >  1:2244/1‑149 (MQ=255)

              cACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCg                                                                                                                                         >  2:354559/1‑150 (MQ=255)

                cGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCgtg                                                                                                                                       <  1:1069127/150‑1 (MQ=255)

                 gCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTAACGAGGGTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCgtgt                                                                                                                                      >  1:1947299/1‑150 (MQ=255)

                 gCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACAGCCGCGAAGTAGTCACCCTGGTTGAATTTAACGAGGTTAATTTCGCTAATGGCTTGCGCACACCCCAAACGCAGCgtgt                                                                                                                                      >  2:1629159/1‑150 (MQ=255)

                  cccGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCCCGAAGTAGTCACCCAGGGTGAATTTACCGATGTTAATTTCGCTACTGGCTTGCGCACACCCCAAACGCAGCGGGtt                                                                                                                                     >  1:683602/1‑150 (MQ=255)

                   ccGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTg                                                                                                                                    <  2:3707300/150‑1 (MQ=255)

                    cGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTTTCTTACCCAGCCGCGAAGTAGTCACCCGGGTTGTATTTAAACAGGTTAAATTAGCTACTGGCTTGCGCACTCCACAAACCCACCGTGTTgg                                                                                                                                   >  1:4133172/1‑150 (MQ=255)

                       aGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTcg                                                                                                                                <  1:3828878/150‑1 (MQ=255)

                        gTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAAAGCAGCGTGTTGGTcgc                                                                                                                               >  2:4092761/1‑150 (MQ=255)

                                aGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATAATCCACAGCCGCGAACTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTAGCTACTGGCTTGCGCACACCACAAACCCACCGTGTTGGTCGCGGCGgtgt                                                                                                                       >  1:2993836/1‑150 (MQ=255)

                                 gTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTc                                                                                                                      <  2:2356717/150‑1 (MQ=255)

                                 gTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCCCCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAAAGCAGCGGGTTGGTCGCGTCGGTGTc                                                                                                                      >  2:617131/1‑150 (MQ=255)

                                 gTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACACCACAAACGCAGCGTGTTGGTCGCGTCGGTGTc                                                                                                                      >  2:945300/1‑150 (MQ=255)

                                  tCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCg                                                                                                                     <  1:3360694/150‑1 (MQ=255)

                                   cATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACAGCCGCGAAGTAGTCACCCGGGTTGAATTTAACGAGGTTAATTTCGCTAAAGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGt                                                                                                                    >  2:2897531/1‑150 (MQ=255)

                                    aTCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTa                                                                                                                   <  2:531655/150‑1 (MQ=255)

                                        ggATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCa                                                                                                               <  1:3202958/150‑1 (MQ=255)

                                          aTACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTACCCAGCCGCGAAGTAGTCACCCTGGTTGTATTTACCGAGGTGAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGGGTTGGTCGCGGCGGTGTCGTTACCAgg                                                                                                             >  2:2089352/1‑150 (MQ=255)

                                           tACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAggg                                                                                                            <  2:295159/150‑1 (MQ=255)

                                               gTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTa                                                                                                        <  2:1386646/150‑1 (MQ=255)

                                               gTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTa                                                                                                        <  2:4575755/150‑1 (MQ=255)

                                                tACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTAt                                                                                                       <  1:4342065/150‑1 (MQ=255)

                                                      agaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATc                                                                                                 <  2:1414577/150‑1 (MQ=255)

                                                       gaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCg                                                                                                <  1:3327243/150‑1 (MQ=255)

                                                         tACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTa                                                                                              <  1:3914500/150‑1 (MQ=255)

                                                         tACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTa                                                                                              <  1:2235424/150‑1 (MQ=255)

                                                           cGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAg                                                                                            <  1:987134/150‑1 (MQ=255)

                                                             ttCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAAATTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAAAGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGAACGTAAGCg                                                                                          >  2:1127759/1‑150 (MQ=255)

                                                              tctcGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGa                                                                                         <  2:1105740/150‑1 (MQ=255)

                                                                tcGGAGTGGTTTTTATAATCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACACCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGAca                                                                                       >  2:4747450/1‑149 (MQ=255)

                                                                tcGGAGTGGTTTTTATAATCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTAACGAGGTTAATTTAGCTAATGGCTTGCGCACACCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATAGTAAGCGACt                                                                                       >  2:4746097/1‑150 (MQ=255)

                                                                     gtgtTTTTTATCTTCCACTGCCGCGATGTAGTCACCCTGGTTTAATTTACCGTGGTTATTTTCGCTAGTGGGTTGCGCACTCAACAAACGCAGCTTGTTGGTCGCGTGGGTTTGGTAACCAGGGATTATCTTATCTTAAGCGACTCCCAg                                                                                  <  2:1334921/150‑1 (MQ=255)

                                                                     gTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTTATTTTACCGTGTTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTGGGTGTGGTAACCAGGGATTATCTTATCGTAAGCGACTCCCAg                                                                                  <  2:2597348/150‑1 (MQ=255)

                                                                          ttttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtct                                                                             >  2:3575672/1‑150 (MQ=255)

                                                                          ttttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtct                                                                             >  2:1017847/1‑150 (MQ=255)

                                                                              aTCTTCCACGGCCGCGATGTAGACACCCTGTTTTATTTTCCCGTGGTTATTTTCGCTAGTGGCTTGCGCACTCCAAAAACGCTGCGTGTTGGTCGCGTCGGTGTCGTAACCCGGTATTATCTTATCGTAAGCGACTCCCAGAATCTCTTc                                                                         <  2:4105303/150‑1 (MQ=255)

                                                                                   ccACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGTATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGttt                                                                    <  2:3424944/150‑1 (MQ=255)

                                                                                       tGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAAt                                                                >  1:3769153/1‑150 (MQ=255)

                                                                                         ccGCGATTTAGTCACCCGGGTTGATTTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATcc                                                              <  2:3464434/150‑1 (MQ=255)

                                                                                         ccGCGAAGTAGTCACCCTGGTTGAATTTACCGAGTTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATcc                                                              <  2:2328965/150‑1 (MQ=255)

                                                                                          cgctatgtagtCACCCTGTTTTAATTTACCGTGTTTATTTTCGCTACGGGCTTGCGCACTCCACAAACGCCGGGTGTTGGGCGCGTGGGTGTCGTAACCTGGTATTATCTGATCTTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCa                                                             <  1:4615527/144‑1 (MQ=255)

                                                                                          cgcgAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCa                                                             >  1:1352913/1‑150 (MQ=255)

                                                                                             gAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAgcg                                                          >  2:4636484/1‑150 (MQ=255)

                                                                                                 ttgtCACCGTGGTTTATTTTACCGTGGTTATTTTCGCTACGGGCTTGCGCACCCCAAAACGGCTGGGTGTTGGTGGCGTCGGTGTCTAAACCTGGTATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtt                                                      <  1:386453/148‑1 (MQ=255)

                                                                                                 tagtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtt                                                      >  1:1219156/1‑150 (MQ=255)

                                                                                                 tagtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtt                                                      >  1:905019/1‑150 (MQ=255)

                                                                                                   gtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtttt                                                    >  2:3966771/1‑150 (MQ=255)

                                                                                                   gtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtttt                                                    >  2:3966887/1‑150 (MQ=255)

                                                                                                      aCCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGttttttt                                                 >  2:2707153/1‑150 (MQ=255)

                                                                                                      aCCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGttttttt                                                 >  1:1130914/1‑150 (MQ=255)

                                                                                                      aCCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGttttttt                                                 >  1:803167/1‑150 (MQ=255)

                                                                                                         gtgttttaaTTTACCGTGTTTATTTGCGCTAGGGGTTTGCGCACCCCAAAAACGCTGGTTTTGGGTCGGGTGGTTGTCTTACCCATTTTTTATATTATCTTACGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTAcc                                              <  1:2500025/143‑1 (MQ=255)

                                                                                                         cTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTAcc                                              >  1:4164210/1‑150 (MQ=255)

                                                                                                           ggTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCtt                                            >  1:2342617/1‑150 (MQ=255)

                                                                                                            gTTTAATTTACCTTGTTTAATTTCGCTAGTGTCTTGCGCACTCAAAAAACGCAGCTTGTTGTTCGCGTCGGTGTCGTAACCAGGTTTTATATGTTCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTc                                           <  2:3956345/150‑1 (MQ=255)

                                                                                                              tGATTTTCCCGTGTTTATTTTCGCTACTGGCTTGCGCACTCCAAAAACGCTGCGTTTGGGTCGCTTCGGTGTCTTAACCTGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCct                                         <  2:2845206/150‑1 (MQ=255)

                                                                                                              tGAATTTACCTAGGTTATTTTCGCTACTGGCTTGCGCACTCCACAAACGCTGCGTGTTGGTCGCGTGGGTGTCTTAACCAGGTTTTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCct                                         <  2:4606178/150‑1 (MQ=255)

                                                                                                              tGAATTTACCGTGTTTAATTTCGCTCCTGGCTTGCGCACTCCACAAACGCTGGGTGTTGGTCGCGTCGGTGTCTTAACCAGGGTTTATCTGATCGAAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCct                                         <  1:2353893/150‑1 (MQ=255)

                                                                                                              tGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCct                                         <  2:1708043/150‑1 (MQ=255)

                                                                                                               gAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctg                                        >  2:4410888/1‑150 (MQ=255)

                                                                                                               gAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctg                                        >  2:3498131/1‑150 (MQ=255)

                                                                                                                 aTTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctgct                                      >  2:1062670/1‑150 (MQ=255)

                                                                                                                   ttACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGa                                    >  1:1682462/1‑150 (MQ=255)

                                                                                                                     acgttgtttATTTTCGCTAGTGGCTTGCGCACTCCACAACCGCCGCGTGTTGGTCGCGTGGGTGTCTTAACCTGGTTTTATCTTATCTAACGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAAt                                  <  1:1991432/143‑1 (MQ=255)

                                                                                                                     accttgtttATTTTCGGTAGTGGCTTGCGCACTCCAAAAACGCAGTTTGTTGGGCGCGTCGTGGTGTTAACCCGGGTTTATGTGTTCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAAt                                  <  2:1233611/143‑1 (MQ=255)

                                                                                                                     aCCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAAt                                  >  2:3675857/1‑150 (MQ=255)

                                                                                                                      ccGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATg                                 >  1:2877787/1‑150 (MQ=255)

                                                                                                                        gAGTTTTATTTCGCGAGTGGCTTGCGAACTCAAAAAACGCTGGGGGTTGGTCGCGTCGTTGTCTTAACCAGGTATTATCTTATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCg                               <  1:2608803/150‑1 (MQ=255)

                                                                                                                        gAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCg                               >  1:237576/1‑150 (MQ=255)

                                                                                                                           gTTATTTTCGCTACGGGCTGGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCCGGGTTTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGAcc                            <  2:4189425/150‑1 (MQ=255)

                                                                                                                              aaTTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGcc                         >  1:3575538/1‑150 (MQ=255)

                                                                                                                              aaTTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTAACTTCCTGCTGAATGCGACCGcc                         >  1:2395150/1‑150 (MQ=255)

                                                                                                                              aaTTTCGATACGGGCTTGCGCACTCAAAAAACGCTGCGTGTTGGTGGCGTCGGTGTCGTAACCCGGTTTTATCTTATCTAAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGcc                         <  2:841650/150‑1 (MQ=255)

                                                                                                                               aTTTCGCTACTGGCTTGCGCCCTCCACAAACGCAGGGTTTTGGTCGCGTCGGGGGCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCa                        <  2:3409180/150‑1 (MQ=255)

                                                                                                                               aTTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCa                        >  2:2971749/1‑150 (MQ=255)

                                                                                                                                 ttCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAACGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCaaa                      >  1:4505855/1‑150 (MQ=255)

                                                                                                                                 ttCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAACGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCaaa                      >  1:4505942/1‑150 (MQ=255)

                                                                                                                                  tCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCaaaa                     >  1:2665652/1‑150 (MQ=255)

                                                                                                                                      tccTGGCTTGCGCACTCCACAAACGCTGGGTGTTGGTCGGGTGGTTGTCTAAACCAGTTATTATCTGATCTTACGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGa                 <  1:4054780/148‑1 (MQ=255)

                                                                                                                                         tGGTTTGCTCACTCCACAAACGCCGCGTTTTGGTCGCTTCGGTGTCTTAACCAGTTATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACtt              <  2:1892836/150‑1 (MQ=255)

                                                                                                                                          ggCTTGCGCACTCAAAAAACGCTGGTTGTTGGTCGGGTGGGTGTCGTAACCCGGGATTATCTTTTGGTACGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACttt             <  1:728816/150‑1 (MQ=255)

                                                                                                                                           gCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTta            >  1:3308813/1‑150 (MQ=255)

                                                                                                                                            cTTGCGCACTCCAAAACCGCAGGGTGTTGGTCGCGTCGGTGTCGTAACCAGGTATTATCTTATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTtat           <  2:1247048/150‑1 (MQ=255)

                                                                                                                                              tGCGCCCACCACACAGGCTGGGTGTGGGGCGCGTGGGTGTGTAAACCTGGGTTTATCTTAGCGTACGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAg         <  1:853347/150‑1 (MQ=255)

                                                                                                                                                  cACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATATcgcg     >  1:1782521/1‑150 (MQ=255)

                                                                                                                                                  cACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTGTATAGcgcg     >  1:428837/1‑150 (MQ=255)

                                                                                                                                                    cTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAGCGCgtg   >  1:930348/1‑150 (MQ=255)

                                                                                                                                                     tCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAGCGCgtgt  <  2:1554695/150‑1 (MQ=255)

ACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGACCGCCAAAACGGACTTTATAGCGCGTGT  >  CP009273/3559090‑3559388

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