BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,745,930	G→A	E406K (GAA→AAA)	paoC →	aromatic aldehyde dehydrogenase molybdopterin‑binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,745,930	0	G	A	100.0%	506.0 / NA	173	E406K (GAA→AAA)	paoC	aromatic aldehyde dehydrogenase molybdopterin‑binding subunit
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (94/79); total (94/79)

Only 100 of 175 total aligned reads displayed.

TTGCCTGTGCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCGAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGCCAGCGCCA  >  NC_002947/3745808‑3746050

ttGCCTGTGCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCaaa                                                                                                                        <  1:3181350/125‑1 (MQ=255)

ttGCCTGTGCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCaaa                                                                                                                        <  2:3181351/125‑1 (MQ=255)

 tGCCTGTGCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAAc                                                                                                                       <  1:3181352/125‑1 (MQ=255)

 tGCCTGTGCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAAc                                                                                                                       >  2:3181470/1‑125 (MQ=255)

   ccTGTGCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGt                                                                                                                     <  2:3181353/125‑1 (MQ=255)

    cTGTGCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGtt                                                                                                                    <  1:3181355/125‑1 (MQ=255)

     tgtgCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGttt                                                                                                                   <  1:3181357/125‑1 (MQ=255)

     tgtgCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGttt                                                                                                                   >  1:3181590/1‑125 (MQ=255)

       tgCAGCGCAGATTGACCCGCTTGACCTGCGGCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTtg                                                                                                                 <  2:3181359/125‑1 (MQ=255)

        gCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTtg                                                                                                                 <  2:3181361/124‑1 (MQ=255)

         cAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTtggc                                                                                                               <  1:3181363/125‑1 (MQ=255)

           gcgcAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTtggctg                                                                                                             >  2:3181681/1‑125 (MQ=255)

           gcgcAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTtggctg                                                                                                             >  2:3181615/1‑125 (MQ=255)

           gcgcAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTtggctg                                                                                                             <  2:3181365/125‑1 (MQ=255)

             gcAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTtggctggc                                                                                                           <  2:3181366/125‑1 (MQ=255)

                 aTTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCaaaa                                                                                                       <  1:3181367/125‑1 (MQ=255)

                 aTTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCaaaa                                                                                                       >  2:3181423/1‑125 (MQ=255)

                  ttGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAAc                                                                                                      >  1:3181665/1‑125 (MQ=255)

                     aCCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGt                                                                                                   >  2:3181674/1‑125 (MQ=255)

                     aCCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGt                                                                                                   <  1:3181369/125‑1 (MQ=255)

                      cccGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTa                                                                                                  <  2:3181371/125‑1 (MQ=255)

                       ccGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTaa                                                                                                 <  2:3181373/125‑1 (MQ=255)

                       ccGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTaa                                                                                                 >  1:3181580/1‑125 (MQ=255)

                         gCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATc                                                                                               <  1:3181375/125‑1 (MQ=255)

                         gCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATc                                                                                               <  1:3181377/125‑1 (MQ=255)

                          cTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATcc                                                                                              >  1:3181645/1‑125 (MQ=255)

                           ttGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCg                                                                                             >  2:3181568/1‑125 (MQ=255)

                             gACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCa                                                                                           >  2:3181517/1‑125 (MQ=255)

                             gACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCa                                                                                           <  1:3181378/125‑1 (MQ=255)

                              aCCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCaa                                                                                          >  2:3181475/1‑125 (MQ=255)

                              aCCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCaa                                                                                          >  2:3181551/1‑125 (MQ=255)

                              aCCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGACAAACGTTTTGGCTGGCAAAACCGTAATCCGCaa                                                                                          >  2:3181663/1‑125 (MQ=255)

                               ccTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAAc                                                                                         <  1:3181381/125‑1 (MQ=255)

                                cTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAAcc                                                                                        >  2:3181569/1‑125 (MQ=255)

                                cTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAAcc                                                                                        <  2:3181384/125‑1 (MQ=255)

                                cTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGACAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAAcc                                                                                        <  2:3181383/125‑1 (MQ=255)

                                  gcgcCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACcgc                                                                                      >  2:3181540/1‑125 (MQ=255)

                                   cgcCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACcggg                                                                                     <  2:3181386/125‑3 (MQ=255)

                                   cgcCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACcgcg                                                                                     <  1:3181387/125‑1 (MQ=255)

                                   cgcCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACcgcg                                                                                     >  2:3181500/1‑125 (MQ=255)

                                        gcgcAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCgcagca                                                                                >  2:3181508/1‑125 (MQ=255)

                                        gcgcAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCgcagca                                                                                >  2:3181394/1‑125 (MQ=255)

                                        gcgcAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCgcagca                                                                                <  2:3181391/125‑1 (MQ=255)

                                        gcgcAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCgcagca                                                                                <  1:3181390/125‑1 (MQ=255)

                                          gcAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATg                                                                              <  1:3181392/125‑1 (MQ=255)

                                           cAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATgc                                                                             >  1:3181574/1‑125 (MQ=255)

                                           cAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATgc                                                                             >  1:3181443/1‑125 (MQ=255)

                                            aaTTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATgcg                                                                            <  2:3181395/125‑1 (MQ=255)

                                             aTTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATgcgc                                                                           >  1:3181516/1‑125 (MQ=255)

                                               tttGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCa                                                                         >  1:3181743/1‑125 (MQ=255)

                                               tttGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCa                                                                         >  2:3181637/1‑125 (MQ=255)

                                               tttGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCa                                                                         <  1:3181397/125‑1 (MQ=255)

                                                ttGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCAt                                                                        <  1:3181399/125‑1 (MQ=255)

                                                ttGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCAt                                                                        >  1:3181466/1‑125 (MQ=255)

                                                 tGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATg                                                                       >  2:3181557/1‑125 (MQ=255)

                                                   ccgccAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGc                                                                     >  2:3181541/1‑125 (MQ=255)

                                                   ccgccAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGc                                                                     <  2:3181401/125‑1 (MQ=255)

                                                   ccgccAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGc                                                                     >  2:3181504/1‑125 (MQ=255)

                                                     gccAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCaa                                                                   >  2:3181455/1‑125 (MQ=255)

                                                      ccAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAAc                                                                  <  2:3181404/125‑1 (MQ=255)

                                                        aGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCa                                                                >  1:3181483/1‑125 (MQ=255)

                                                        aGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCa                                                                >  1:3181705/1‑125 (MQ=255)

                                                        aGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCa                                                                <  1:3181405/125‑1 (MQ=255)

                                                           aaTGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCt                                                             >  1:3181468/1‑125 (MQ=255)

                                                            aTGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTg                                                            >  1:3181788/1‑125 (MQ=255)

                                                              gCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGt                                                          <  1:3181406/125‑1 (MQ=255)

                                                                cAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCg                                                        >  2:3181589/1‑125 (MQ=255)

                                                                cAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCg                                                        <  2:3181409/125‑1 (MQ=255)

                                                                   cGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCt                                                     <  2:3181410/125‑1 (MQ=255)

                                                                    gAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCCGCATGCGCCATGGCAACCAGCTGGTCgg                                                       >  2:3181630/1‑122 (MQ=255)

                                                                     tgggCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCACGCAAAACCGCAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGc                                                              >  2:3181662/2‑114 (MQ=255)

                                                                      gggCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTggg                                                  <  2:3181412/125‑1 (MQ=255)

                                                                       ggCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTgggg                                                 >  1:3181640/1‑125 (MQ=255)

                                                                         cAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGa                                               <  1:3181413/125‑1 (MQ=255)

                                                                          acgCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGc                                                      >  2:3181586/3‑117 (MQ=255)

                                                                            gCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATgg                                            <  1:3181416/125‑1 (MQ=255)

                                                                                tACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCcgg                                        >  1:3181458/1‑125 (MQ=255)

                                                                                  cTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCcggcg                                      >  1:3181729/1‑125 (MQ=255)

                                                                                       gCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTc                                 <  2:3181418/125‑1 (MQ=255)

                                                                                          aGGAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTgg                              >  1:3181604/1‑125 (MQ=255)

                                                                                            gAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGa                            >  1:3181571/1‑125 (MQ=255)

                                                                                            gAACTGCTGACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGa                            <  1:3181420/125‑1 (MQ=255)

                                                                                               ctgctgACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGc                         >  2:3181722/1‑125 (MQ=255)

                                                                                                 gctgACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCa                       <  1:3181421/125‑1 (MQ=255)

                                                                                                  ctgACGTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCAt                      >  1:3181837/1‑125 (MQ=255)

                                                                                                       gTGCTACCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGa                 <  1:3181423/125‑1 (MQ=255)

                                                                                                          cTACCGGCAAGGTGCCAAACGTATTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTggg                                                  >  2:3181732/1‑89 (MQ=255)

                                                                                                            aCCGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAg            <  2:3181425/125‑1 (MQ=255)

                                                                                                             ccGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAgg           <  1:3181426/125‑1 (MQ=255)

                                                                                                             ccGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAgg           <  1:3181428/125‑1 (MQ=255)

                                                                                                              cGGCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGc          >  2:3181488/1‑125 (MQ=255)

                                                                                                               ggCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGcc         >  2:3181610/1‑125 (MQ=255)

                                                                                                               ggCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGcc         >  2:3181717/1‑125 (MQ=255)

                                                                                                                gCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGc                                                                                            >  2:3181739/1‑41 (MQ=255)

                                                                                                                gCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGCCa        >  1:3181657/1‑125 (MQ=255)

                                                                                                                gCAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGCCa        >  1:3181611/1‑125 (MQ=255)

                                                                                                                 cAAGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGCCAg       >  2:3181537/1‑125 (MQ=255)

                                                                                                                  aaGGTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGCCAgc      <  2:3181430/125‑1 (MQ=255)

                                                                                                                     gTGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGCCAGCGcc   >  1:3181688/1‑125 (MQ=255)

                                                                                                                      tGCCAAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGCCAGCGCCa  <  1:3181431/125‑1 (MQ=255)

TTGCCTGTGCAGCGCAGATTGACCCGCTTGACCTGCGCCGGCGCAATTTTGCCGCCAGCAATGCCAACGAGGGCAAGCCCTACTCCAGCAAGGAACTGCTGACGTGCTACCGGCAAGGTGCCGAACGTTTTGGCTGGCAAAACCGTAATCCGCAACCGCGCAGCATGCGCCATGGCAACCAGCTGGTCGGCTGGGGGATGGCCGGCGGCGTCTGGGAAGCCATGCAGATGAAGGCCAGCGCCA  >  NC_002947/3745808‑3746050

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 28 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`