BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,655,360	(C)_5→6	coding (314/327 nt)	ynfA ←	conserved inner membrane protein YnfA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,655,355	1	.	C	100.0%	1101.7 / NA	263	R107G (CGC→GGC)	ynfA	conserved inner membrane protein YnfA
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (174/89); total (174/89)

Only 100 of 298 total aligned reads displayed.

TTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCG‑CCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCA  >  NC_000913/1655255‑1655450

ttGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTgcg                                                                                                  >  1:3181661/1‑101 (MQ=255)

ttGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTgcg                                                                                                  <  2:3181155/101‑1 (MQ=255)

 tGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTgcg‑c                                                                                                <  1:3181235/101‑1 (MQ=255)

 tGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTgcg‑c                                                                                                >  2:3181664/1‑101 (MQ=255)

  gcgcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCG‑cc                                                                                               <  1:3181178/101‑1 (MQ=255)

    gcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCG‑cccc                                                                                             >  2:3181671/1‑101 (MQ=255)

    gcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCG‑cccc                                                                                             >  2:3181667/1‑101 (MQ=255)

     cTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCG‑ccccc                                                                                            <  1:3181015/101‑1 (MQ=255)

     cTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCG‑ccccc                                                                                            >  1:3181672/1‑101 (MQ=255)

       gtgAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCG‑cccccca                                                                                          <  2:3181023/101‑3 (MQ=255)

       gtgAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCG‑cccccca                                                                                          >  2:3181677/1‑99 (MQ=255)

         gAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCGTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCG‑ccccccagc                                                                                        >  1:3181679/1‑97 (MQ=255)

           aTGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGccc                                                                                       >  1:3181682/1‑101 (MQ=255)

            tGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCg                                                                                      <  2:3181284/101‑1 (MQ=255)

            tGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCg                                                                                      <  2:3180587/101‑1 (MQ=255)

            tGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCg                                                                                      >  2:3181684/1‑101 (MQ=255)

              gTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGca                                                                                    <  1:3181144/101‑1 (MQ=255)

                tGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGcaac                                                                                  >  1:3181689/1‑101 (MQ=255)

                 gACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGcaaca                                                                                 <  2:3181180/101‑1 (MQ=255)

                  aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGcaacaa                                                                                >  2:3181690/1‑101 (MQ=255)

                    aaaaGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATg                                                                              >  1:3181695/1‑101 (MQ=255)

                     aaaGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGa                                                                             >  1:3181698/1‑101 (MQ=255)

                      aaGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGAt                                                                            <  2:3181133/101‑1 (MQ=255)

                      aaGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGAt                                                                            >  2:3181701/1‑101 (MQ=255)

                       aGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATc                                                                           >  1:3181702/1‑101 (MQ=255)

                         aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCaa                                                                         >  1:3181707/1‑101 (MQ=255)

                         aTGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCaa                                                                         >  2:3181709/1‑101 (MQ=255)

                          tGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAAc                                                                        >  1:3181712/1‑101 (MQ=255)

                           gAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACa                                                                       >  2:3181714/1‑101 (MQ=255)

                           gAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACa                                                                       >  2:3181717/1‑101 (MQ=255)

                              aTAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGc                                                                    >  2:3181722/1‑101 (MQ=255)

                               tAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGcc                                                                   >  1:3181726/1‑101 (MQ=255)

                               tAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGcc                                                                   >  2:3181725/1‑101 (MQ=255)

                                 gAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGc                                                                 >  1:3181731/1‑101 (MQ=255)

                                 gAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGc                                                                 <  2:3181292/101‑1 (MQ=255)

                                   aTACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCaa                                                               >  1:3181735/1‑101 (MQ=255)

                                   aTACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCaa                                                               <  1:3181521/101‑1 (MQ=255)

                                    tACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCaaa                                                              >  2:3181736/1‑101 (MQ=255)

                                     aCCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAg                                                             <  1:3181309/101‑1 (MQ=255)

                                      ccTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAgc                                                            >  2:3181740/1‑101 (MQ=255)

                                      ccTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAgc                                                            <  1:3181482/101‑1 (MQ=255)

                                        ttttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAgcgc                                                          <  2:3181195/101‑1 (MQ=255)

                                        ttttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAgcgc                                                          >  2:3181744/1‑101 (MQ=255)

                                        ttttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAgcgc                                                          >  1:3181747/1‑101 (MQ=255)

                                         tttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCa                                                         >  2:3181749/1‑101 (MQ=255)

                                           tGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAAt                                                       >  1:3181752/1‑101 (MQ=255)

                                           tGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAAt                                                       >  2:3181754/1‑101 (MQ=255)

                                            gTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAATCGCAATc                                                      <  2:3181080/101‑1 (MQ=255)

                                             tCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCa                                                     <  2:3181301/101‑1 (MQ=255)

                                              cAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCaa                                                    >  2:3181758/1‑101 (MQ=255)

                                               aGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAAc                                                   <  2:3181440/101‑1 (MQ=255)

                                                gCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACg                                                  >  1:3181762/1‑101 (MQ=255)

                                                 cTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGc                                                 >  2:3181763/1‑101 (MQ=255)

                                                  tGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCa                                                <  1:3181272/101‑1 (MQ=255)

                                                  tGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCa                                                >  2:3181766/1‑101 (MQ=255)

                                                  tGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCa                                                <  2:3181361/101‑1 (MQ=255)

                                                  tGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCa                                                <  2:3181271/101‑1 (MQ=255)

                                                      acTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCg                                            >  2:3181771/1‑101 (MQ=255)

                                                       cTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGt                                           >  1:3181774/1‑101 (MQ=255)

                                                        ttCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTc                                          >  2:3181775/1‑101 (MQ=255)

                                                         tCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTcc                                         <  2:3181206/101‑1 (MQ=255)

                                                          ccTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                        >  2:3181778/1‑101 (MQ=255)

                                                          ccTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                        <  2:3181391/101‑1 (MQ=255)

                                                          ccTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                        <  1:3181363/101‑1 (MQ=255)

                                                            tcttatcttaTTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGt                                      <  1:3181405/101‑1 (MQ=255)

                                                            tcttatcttaTTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGt                                      <  1:3181437/101‑1 (MQ=255)

                                                            tcttatcttaTTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGt                                      >  2:3181781/1‑101 (MQ=255)

                                                             cttatcttatTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTc                                     >  2:3181784/1‑101 (MQ=255)

                                                               tatcttatTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCAt                                   <  2:3181247/101‑1 (MQ=255)

                                                               tatcttatTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCAt                                   >  2:3181789/1‑101 (MQ=255)

                                                                atcttatTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATa                                  >  1:3181790/1‑101 (MQ=255)

                                                                  cttatTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATaaa                                >  1:3181794/1‑101 (MQ=255)

                                                                    tatTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAaga                              >  1:3181797/1‑101 (MQ=255)

                                                                     atTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAagag                             >  2:3181799/1‑101 (MQ=255)

                                                                     atTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAagag                             <  2:3181229/101‑1 (MQ=255)

                                                                      ttGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTcc                                         >  2:3181804/1‑88 (MQ=255)

                                                                       tGATAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAagag                             >  1:3181805/1‑99 (MQ=255)

                                                                         aTAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAAc                                                                        >  1:3181811/1‑54 (MQ=255)

                                                                         aTAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAagag                             >  2:3181810/1‑97 (MQ=255)

                                                                         aTAAAATGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAg                         <  1:3181290/101‑1 (MQ=255)

                                                                           aaaaTGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGtt                       >  1:3181816/1‑101 (MQ=255)

                                                                            aaaTGGATTTATGTTCCTACGGGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATa                                  >  2:3181818/1‑89 (MQ=255)

                                                                             aaTGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTc                     >  1:3181820/1‑101 (MQ=255)

                                                                              aTGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCa                    <  1:3181640/101‑1 (MQ=255)

                                                                              aTGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCa                    <  2:3181177/101‑1 (MQ=255)

                                                                               tGGATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCAc                   <  2:3181542/101‑1 (MQ=255)

                                                                                ggATTTATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACg                  >  1:3181825/1‑101 (MQ=255)

                                                                                   tttATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCa               >  2:3181828/1‑101 (MQ=255)

                                                                                    ttATGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCAt              <  1:3181346/101‑1 (MQ=255)

                                                                                      aTGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATcc            >  1:3181831/1‑101 (MQ=255)

                                                                                       tGTTCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCca           >  2:3181834/1‑101 (MQ=255)

                                                                                         ttCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCcaca         >  1:3181836/1‑101 (MQ=255)

                                                                                          tCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACaa        >  2:3181839/1‑101 (MQ=255)

                                                                                          tCCTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACaa        <  2:3181326/101‑1 (MQ=255)

                                                                                           ccTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATc                                                      <  1:3181060/54‑1 (MQ=255)

                                                                                            cTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAAc                                                   >  2:3181844/1‑56 (MQ=255)

                                                                                            cTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAAcg      >  1:3181842/1‑101 (MQ=255)

                                                                                            cTACGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAAcg      <  1:3181377/101‑1 (MQ=255)

                                                                                              aCGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAAcgcg    >  2:3181849/1‑101 (MQ=255)

                                                                                               cGTGCGCCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACgcgc   >  1:3181852/1‑101 (MQ=255)

TTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCG‑CCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCA  >  NC_000913/1655255‑1655450

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