BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	8,415	A→C	31.5%	K60Q (AAA→CAA) ‡	talB →	transaldolase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	8,415	0	A	C	31.5%	144.8 / 104.4	95	K60Q (AAA→CAA) ‡	talB	transaldolase B
Reads supporting (aligned to +/- strand): ref base A (32/33); new base C (16/14); total (48/47)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01

CACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCA  >  CP009273/8270‑8561

caccacCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGaaa                                                                                                                                                  >  1:672247/1‑148 (MQ=255)

caccacCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAAc                                                                                                                                                 <  1:633300/149‑1 (MQ=255)

 accaccGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAAcag                                                                                                                                               <  2:1052834/150‑1 (MQ=255)

  ccaccGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAAcagc                                                                                                                                              <  2:785236/150‑1 (MQ=255)

     ccGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCaga                                                                                                                                           <  2:817462/150‑1 (MQ=255)

        tAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCa                                                                                                                                        <  1:1026921/150‑1 (MQ=255)

        tAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCa                                                                                                                                        <  1:771332/150‑1 (MQ=255)

         aGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCaa                                                                                                                                       >  2:918938/1‑150 (MQ=255)

          gTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAAc                                                                                                                                      >  1:528165/1‑150 (MQ=255)

          gTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAAc                                                                                                                                      <  2:1075961/150‑1 (MQ=255)

               cGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATcg                                                                                                                                 <  2:834331/150‑1 (MQ=255)

                gACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATcgc                                                                                                                                <  2:720596/150‑1 (MQ=255)

                   acTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATcgcgcg                                                                                                                             >  1:1142939/1‑150 (MQ=255)

                       gggACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCgcagc                                                                                                                         <  1:725406/150‑1 (MQ=255)

                           tattgcgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCgcag                                                                                                                          >  1:253692/5‑145 (MQ=255)

                           tattgcgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCgcag                                                                                                                          >  1:1132482/5‑145 (MQ=255)

                           tattgcgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCgcag                                                                                                                          >  2:420077/5‑145 (MQ=255)

                           cATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGAt                                                                                                                     <  1:779517/150‑1 (MQ=255)

                            attgcgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATc                                                                                                                    >  2:955598/4‑150 (MQ=255)

                            aTCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATc                                                                                                                    <  1:918938/150‑1 (MQ=255)

                                cgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctcg                                                                                                                <  2:315741/150‑5 (MQ=255)

                                cgGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctc                                                                                                                 >  1:510553/1‑146 (MQ=255)

                                 ggCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGa                                                                                                               >  1:325155/1‑150 (MQ=255)

                                  gCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGAc                                                                                                              <  1:787570/150‑1 (MQ=255)

                                   cAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctcgatg                                                                                                             >  1:839056/1‑143 (MQ=255)

                                   cAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATcctcgatg                                                                                                             >  2:402266/1‑143 (MQ=255)

                                   cAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGAcg                                                                                                             >  1:157184/1‑150 (MQ=255)

                                      tGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCgca                                                                                                                           >  1:133896/1‑133 (MQ=255)

                                       gAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCgcag                                                                                                                          <  2:133896/133‑1 (MQ=255)

                                       gAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGcgac                                                                                                         >  2:17371/1‑150 (MQ=255)

                                        aaGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGcgacc                                                                                                        <  2:5709/150‑1 (MQ=255)

                                        aaGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGcgacc                                                                                                        <  1:1122061/150‑1 (MQ=255)

                                              tATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACaaa                                                                                                  >  2:425519/1‑150 (MQ=255)

                                                  aaCCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTg                                                                                               <  1:1034559/149‑1 (MQ=255)

                                                   gcctcAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGc                                                                                             >  2:904138/5‑150 (MQ=255)

                                                    cctcAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGcc                                                                                            <  1:277255/147‑2 (MQ=255)

                                                    ccGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTgg                                                                                              >  2:363146/1‑148 (MQ=255)

                                                       cAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTa                                                                                         <  2:706123/150‑1 (MQ=255)

                                                       cAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTa                                                                                         <  2:313671/150‑1 (MQ=255)

                                                         ggATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                                                                       >  2:594738/1‑150 (MQ=255)

                                                         ggATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                                                                       <  2:962763/150‑1 (MQ=255)

                                                          gATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                                                                       >  2:589407/1‑149 (MQ=255)

                                                          gATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTaaa                                                                                       >  1:967438/1‑149 (MQ=255)

                                                           aTGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTgaaca                                                                                     >  2:738123/1‑148 (MQ=255)

                                                           aTGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTgaaca                                                                                     >  2:1109305/1‑148 (MQ=255)

                                                           aTGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAata                                                                                     <  1:979128/150‑1 (MQ=255)

                                                            tGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAacat                                                                                    <  1:610249/150‑1 (MQ=255)

                                                              ccACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTg                                                                                  >  1:965700/1‑150 (MQ=255)

                                                               caccacTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCgg                                                                                 >  2:442801/1‑150 (MQ=255)

                                                               cacaACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTgg                                                                                 >  2:582833/1‑150 (MQ=255)

                                                                accactaaccCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGt                                                                                <  2:182023/144‑1 (MQ=255)

                                                                acaACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGt                                                                                >  2:248138/1‑150 (MQ=255)

                                                                   actaaccCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCtc                                                                             >  2:950846/4‑149 (MQ=255)

                                                                   actaaccCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCtc                                                                             >  1:443385/4‑149 (MQ=255)

                                                                        cccTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAAt                                                                        <  2:1137026/150‑1 (MQ=255)

                                                                        cccTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAAt                                                                        <  1:796092/150‑1 (MQ=255)

                                                                         ccTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATc                                                                       >  1:479670/1‑150 (MQ=255)

                                                                             ctctcATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTg                                                                    <  2:678858/149‑1 (MQ=255)

                                                                                  aTTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAGCAGCAGAGCAACGATCGCGCGCAGCAGGTCGTGGGCGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAAc                                                                <  2:442265/148‑1 (MQ=255)

                                                                                    tCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGt                                                            >  2:123565/1‑150 (MQ=255)

                                                                                     ctcggcgcggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTgat                                                            <  2:398679/144‑3 (MQ=255)

                                                                                     cTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGt                                                            <  2:1191317/149‑1 (MQ=255)

                                                                                        ggcgcggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCg                                                        <  1:353151/148‑1 (MQ=255)

                                                                                        ggcgcggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCg                                                        <  2:41865/148‑1 (MQ=255)

                                                                                        ggcgcggctcagATCCCTGAGTAACGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCg                                                        <  2:1057072/148‑1 (MQ=255)

                                                                                        aaCGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCg                                                        >  2:303597/1‑150 (MQ=255)

                                                                                         aCGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCgg                                                       <  1:1192639/150‑1 (MQ=255)

                                                                                            cggctcagATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGcc                                                    <  1:825723/145‑1 (MQ=255)

                                                                                               cgcagATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTa                                                 >  2:749744/1‑150 (MQ=255)

                                                                                                gcagATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTAt                                                >  1:376412/1‑150 (MQ=255)

                                                                                                  agATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTAtct                                              <  1:155379/150‑1 (MQ=255)

                                                                                                   gATCCCTGAGTACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTAtctc                                             >  1:803817/1‑150 (MQ=255)

                                                                                                     ttCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCaa                                           >  1:347520/1‑150 (MQ=255)

                                                                                                        ctgagtACCGTAAGCTGATTGACGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCAccg                                        <  2:258949/145‑3 (MQ=255)

                                                                                                         ggAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTg                                        >  2:67474/1‑149 (MQ=255)

                                                                                                         ggAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGa                                       <  1:1126570/150‑1 (MQ=255)

                                                                                                           aaTACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAg                                     <  2:564387/150‑1 (MQ=255)

                                                                                                               ccGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGa                                 <  1:1007092/150‑1 (MQ=255)

                                                                                                               ccGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGa                                 >  1:388733/1‑150 (MQ=255)

                                                                                                                   aaGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATgcg                             >  1:583578/1‑150 (MQ=255)

                                                                                                                   aaGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATgcg                             >  2:379514/1‑150 (MQ=255)

                                                                                                                    agctgattgacGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACgcac                            <  2:724435/147‑3 (MQ=255)

                                                                                                                    agctgattgacGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACgcac                            <  1:976071/147‑3 (MQ=255)

                                                                                                                       tgattgacGACGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTc                         >  2:675492/1‑150 (MQ=255)

                                                                                                                         attgatGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCtt                       >  2:488927/1‑150 (MQ=255)

                                                                                                                           tgatgaTGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTc                     >  1:341799/1‑150 (MQ=255)

                                                                                                                           tgacgaCGCTGTCGCCTGGGCTCGTAGCCAGAGCAGCGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGtctctc                     >  1:201575/1‑150 (MQ=255)

                                                                                                                            gatgatGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTcc                    >  1:422122/1‑150 (MQ=255)

                                                                                                                               gatgCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTAt                 >  1:1095627/1‑150 (MQ=255)

                                                                                                                                   cTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGaca             >  2:1190830/1‑150 (MQ=255)

                                                                                                                                      tCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCg          <  2:672247/150‑1 (MQ=255)

                                                                                                                                      tCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTAGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACAcc           <  2:1142939/149‑1 (MQ=255)

                                                                                                                                          cTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGc      <  2:967438/150‑1 (MQ=255)

                                                                                                                                              gCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTc   >  1:237286/1‑149 (MQ=255)

                                                                                                                                              gCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCa  <  1:409302/150‑1 (MQ=255)

                                                                                                                                                gAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGa               >  2:1026255/1‑135 (MQ=255)

                                                                                                                                                 aaaCAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGac              <  1:1026255/135‑1 (MQ=255)

CACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCGCCTGGGCGAAACAGCAGAGCAACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCA  >  CP009273/8270‑8561

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`