BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,411,047	C→T	22.3%	S165S (AGC→AGT)	yjfJ →	PspA/IM30 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,411,047	0	C	T	22.3%	365.0 / 89.0	179	S165S (AGC→AGT)	yjfJ	PspA/IM30 family protein
Reads supporting (aligned to +/- strand): ref base C (66/73); new base T (21/19); total (87/92)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

CGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGAC  >  NC_000913/4410948‑4411146

tgctcaaactGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGc                                                                                                     <  2:605877/91‑1 (MQ=255)

cGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGc                                                                                                     <  1:2971488/100‑1 (MQ=255)

 gCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGc                                                                                                     >  2:1252785/1‑99 (MQ=255)

 gCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCg                                                                                                    <  1:1483366/100‑1 (MQ=255)

 gCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCg                                                                                                    <  1:54721/100‑1 (MQ=255)

  cTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGt                                                                                                   >  1:102480/1‑100 (MQ=255)

   tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGtt                                                                                                  <  2:3071431/100‑1 (MQ=255)

   tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGtt                                                                                                  >  1:1769312/1‑100 (MQ=255)

    cAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGttt                                                                                                 <  2:549424/100‑1 (MQ=255)

    cAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGttt                                                                                                 <  1:2735663/100‑1 (MQ=255)

     aGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGttt                                                                                                 >  2:3215496/1‑99 (MQ=255)

     aGTTTGAGCAGCAAATGGAAGTCGTCAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTc                                                                                                >  2:3680448/1‑100 (MQ=255)

      gTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACATGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCg                                                                                               <  1:1391771/100‑1 (MQ=255)

      gTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCg                                                                                               <  2:921430/100‑1 (MQ=255)

      gTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCg                                                                                               <  1:447726/100‑1 (MQ=255)

      gTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCg                                                                                               <  1:1942763/100‑1 (MQ=255)

        ttGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGAc                                                                                             <  1:2087308/100‑1 (MQ=255)

        ttGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGAc                                                                                             >  2:2166153/1‑100 (MQ=255)

        ttGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGAc                                                                                             <  2:2618662/100‑1 (MQ=255)

          gAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAg                                                                                           <  2:1665282/100‑1 (MQ=255)

             cagcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCgg                                                                                        <  1:249660/100‑1 (MQ=255)

             cagcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCgg                                                                                        <  1:1327520/100‑1 (MQ=255)

              agcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCgg                                                                                        <  2:1572779/99‑1 (MQ=255)

              agcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGc                                                                                       <  1:710958/100‑1 (MQ=255)

              agcaAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTATACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGc                                                                                       >  1:3245052/1‑100 (MQ=255)

                caAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAg                                                                                     >  1:3297221/1‑100 (MQ=255)

                        agacGTTAACGCCACTGAAGGCCTGCAGCGTGAACAACAGGCTGTACGAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTa                                                                             >  2:2066911/4‑100 (MQ=255)

                        aGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTa                                                                             <  2:1178138/100‑1 (MQ=255)

                        aGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTa                                                                             <  1:1265312/100‑1 (MQ=255)

                         gTCGTTCAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTaa                                                                            >  2:2498642/1‑100 (MQ=255)

                         gTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTa                                                                             >  2:365266/1‑99 (MQ=255)

                         gTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTaa                                                                            >  1:1891463/1‑100 (MQ=255)

                         gTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTaa                                                                            >  1:1797108/1‑100 (MQ=255)

                         gTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTaa                                                                            >  1:3799087/1‑100 (MQ=255)

                           cGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTaaaa                                                                          >  2:2006172/1‑100 (MQ=255)

                            gTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAAc                                                                         <  2:1529723/100‑1 (MQ=255)

                                aaGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCt                                                                     >  2:474551/1‑100 (MQ=255)

                                 aGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTg                                                                    <  2:1412195/100‑1 (MQ=255)

                                 aGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTg                                                                    <  1:726757/100‑1 (MQ=255)

                                  gCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGc                                                                   >  2:293595/1‑100 (MQ=255)

                                    cACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCaa                                                                 >  1:262173/1‑100 (MQ=255)

                                    cACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCaa                                                                 >  2:3164566/1‑100 (MQ=255)

                                     aCTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCaaa                                                                <  1:3059775/100‑1 (MQ=255)

                                     aCTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCaaa                                                                >  1:1836340/1‑100 (MQ=255)

                                      cTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACACCGGACGAATCCTTAAAACGCCTGCAAAc                                                               >  2:3308953/1‑100 (MQ=255)

                                      cTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCCCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAAc                                                               <  2:171066/100‑1 (MQ=255)

                                      cTGAAGCCATACAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAAc                                                               <  1:2000098/100‑1 (MQ=255)

                                       tGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAAcg                                                              >  2:749253/1‑100 (MQ=255)

                                         aaGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAAcgcg                                                            <  2:3745070/100‑1 (MQ=255)

                                         aaGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAAcgcg                                                            >  2:1374753/1‑100 (MQ=255)

                                          aGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGt                                                           >  2:303168/1‑100 (MQ=255)

                                           gCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTc                                                          >  1:3105559/1‑100 (MQ=255)

                                             cATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAg                                                        <  1:1554559/100‑1 (MQ=255)

                                              aTGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAgg                                                       >  2:516343/1‑100 (MQ=255)

                                               tGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGc                                                      <  2:2183194/100‑1 (MQ=255)

                                                gCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGcn                                                     <  2:3535356/100‑2 (MQ=255)

                                                gCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGcc                                                     <  1:1026704/100‑1 (MQ=255)

                                                gCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGcc                                                     >  1:2323347/1‑100 (MQ=255)

                                                 cAGCGTGCACAACAGGCTGTAACAACCTCTAGAGTTTGCGCATCTGCCATCGTTTCGACAGCGGCAGCATCCTCAAAACGCCTGCAAACGCGGCAGGc                                                      <  2:2697890/98‑1 (MQ=255)

                                                 cAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCg                                                    <  2:1891452/100‑1 (MQ=255)

                                                  aGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGa                                                   >  2:1093454/1‑100 (MQ=255)

                                                   gCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGaa                                                  <  1:1232023/100‑1 (MQ=255)

                                                    cGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAAc                                                 <  2:833259/100‑1 (MQ=255)

                                                        cacaACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCa                                             <  1:932603/100‑1 (MQ=255)

                                                        cacaACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCa                                             >  2:1328712/1‑100 (MQ=255)

                                                         acaacaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAg                                            <  1:2955837/100‑1 (MQ=255)

                                                          caacaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAgg                                           >  1:3197690/1‑100 (MQ=255)

                                                            acaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCt                                         >  2:3064843/1‑100 (MQ=255)

                                                              aGGCTGTGACAACCTCTACCTGTGGCCAATTTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAATCGCCTGCAAACGCGTCAGGCCGCACGTCAGGc                                          <  2:3585675/97‑1 (MQ=255)

                                                                gCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGcc                                     <  2:2867372/100‑1 (MQ=255)

                                                                  tGTAACAACCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTg                                   <  1:2166167/100‑1 (MQ=255)

                                                                   gTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTgg                                  >  2:1969362/1‑100 (MQ=255)

                                                                     aacaacCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCGAGCCTgacg                                <  1:2811915/100‑4 (MQ=255)

                                                                     aacaacCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGAt                                >  2:1571856/1‑100 (MQ=255)

                                                                       caacCTCTACCGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGAtgc                              <  1:2080225/100‑1 (MQ=255)

                                                                       caacCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGAtgc                              <  1:1359095/100‑1 (MQ=255)

                                                                           ctctACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGcc                          <  2:687050/100‑1 (MQ=255)

                                                                           ctctACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGcc                          <  2:859032/100‑1 (MQ=255)

                                                                           ctctACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGcc                          >  1:3384877/1‑100 (MQ=255)

                                                                             ctACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGc                        <  1:3678676/100‑1 (MQ=255)

                                                                               aCCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGcac                      >  2:3783377/1‑100 (MQ=255)

                                                                                ccGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAACCGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGca                       <  2:3318829/98‑1 (MQ=255)

                                                                                ccGTTGGCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGcaca                     >  2:3295766/1‑100 (MQ=255)

                                                                                 cGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAg                    <  1:111839/100‑1 (MQ=255)

                                                                                 cGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAg                    <  1:529373/100‑1 (MQ=255)

                                                                                    tGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTg                 <  1:2233474/100‑1 (MQ=255)

                                                                                     ggCGCATCTTCCAGTGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTgg                >  1:3580521/1‑100 (MQ=255)

                                                                                     ggCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTgg                >  1:2665593/1‑100 (MQ=255)

                                                                                     ggCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTgg                >  1:797325/1‑100 (MQ=255)

                                                                                       cgcATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGgag              >  2:2546540/1‑100 (MQ=255)

                                                                                         cATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGaa            <  1:1333814/100‑1 (MQ=255)

                                                                                            cTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGt         >  2:1297640/1‑100 (MQ=255)

                                                                                            cTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGt         <  2:2421511/100‑1 (MQ=255)

                                                                                               ccAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGc      <  1:2035197/100‑1 (MQ=255)

                                                                                               ccAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGc      <  1:2921466/100‑1 (MQ=255)

                                                                                                cAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCa     >  2:2906158/1‑100 (MQ=255)

                                                                                                cAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCa     <  1:880398/100‑1 (MQ=255)

                                                                                                cAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCa     >  2:3416141/1‑100 (MQ=255)

                                                                                                 aGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGTATGCTGCCGCACAGTTGGAGAAAGTCGCAg    >  2:3591787/1‑100 (MQ=255)

                                                                                                   cGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGAc  >  2:300933/1‑100 (MQ=255)

CGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGAC  >  NC_000913/4410948‑4411146

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