BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,892,505	T→C	G75G (GGT→GGC)	yieF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,892,505	0	T	C	98.0%	461.9 / ‑2.0	151	G75G (GGT→GGC)	yieF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (71/77); total (73/78)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 163 total aligned reads displayed.

CATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAC  >  NC_000913/3892415‑3892593

cATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGgcgt                                                                                        <  2:974471/93‑1 (MQ=255)

cATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAAGCTGACGgcgt                                                                                        >  2:2549807/1‑93 (MQ=255)

  ttCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtgg                                                                                      >  1:2651788/1‑93 (MQ=255)

  ttCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtgg                                                                                      >  1:282071/1‑93 (MQ=255)

  ttCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtgg                                                                                      >  2:3315085/1‑93 (MQ=255)

  ttCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtgg                                                                                      <  1:488645/93‑1 (MQ=255)

      cTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGAt                                                                                  <  1:4494194/93‑1 (MQ=255)

      cTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGAt                                                                                  <  1:528431/93‑1 (MQ=255)

              aCGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGc                                                                          <  1:2180492/93‑1 (MQ=255)

                gCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTTGTGATCGTCccgccg                                                                        >  1:3781215/1‑93 (MQ=255)

                gCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCg                                                                        <  1:4339339/93‑1 (MQ=255)

                 cTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCgg                                                                       <  1:626950/93‑1 (MQ=255)

                 cTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCgg                                                                       <  1:2761534/93‑1 (MQ=255)

                   gACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGaa                                                                     >  2:292480/1‑93 (MQ=255)

                    aCGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAat                                                                    <  1:162511/93‑1 (MQ=255)

                     cGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAata                                                                   <  2:2618282/93‑1 (MQ=255)

                      gTACAGCAGGAAGAAGGGTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGCCGGCGTGGTGATCGTCACGCCGGAatat                                                                  >  2:4307068/1‑93 (MQ=255)

                       tACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAAtata                                                                 <  1:4110493/93‑1 (MQ=255)

                       tACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGGGATCGTCACGCCGGAAtata                                                                 >  1:1246641/1‑93 (MQ=255)

                        aCAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATaa                                                                <  1:4076999/93‑1 (MQ=255)

                         cagcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAac                                                               <  2:4631586/93‑1 (MQ=255)

                         cagcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAac                                                               >  2:672527/1‑93 (MQ=255)

                           gcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAacta                                                             <  2:2277121/93‑1 (MQ=255)

                            cagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAactac                                                            <  2:3744175/93‑1 (MQ=255)

                             agGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAactact                                                           <  2:4188269/93‑1 (MQ=255)

                             agGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAactact                                                           >  2:2998273/1‑93 (MQ=255)

                             agGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAactact                                                           <  1:1520924/93‑1 (MQ=255)

                              ggAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTc                                                          <  1:204130/93‑1 (MQ=255)

                              ggAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTc                                                          <  2:1606860/93‑1 (MQ=255)

                              ggAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAAATACTc                                                          >  2:3122664/1‑93 (MQ=255)

                                aagaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCgg                                                        <  1:2267395/93‑1 (MQ=255)

                                aagaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCgg                                                        >  2:1959845/1‑93 (MQ=255)

                                aagaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCgg                                                        >  2:2823322/1‑93 (MQ=255)

                                 agaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGt                                                       >  2:2871039/1‑93 (MQ=255)

                                 agaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCAGTCAGGCTGACGGCGTGGTGAGCGGCACGACGGAATATCACTACTCGGt                                                       >  1:2821713/1‑93 (MQ=255)

                                   aagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAc                                                     <  2:4507507/93‑1 (MQ=255)

                                   aagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGGAc                                                     >  1:873247/1‑93 (MQ=255)

                                    agGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAcc                                                    >  2:2337389/1‑93 (MQ=255)

                                    agGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAcc                                                    <  1:3218652/93‑1 (MQ=255)

                                    agGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAcc                                                    <  2:4296920/93‑1 (MQ=255)

                                    agGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAcc                                                    >  1:1854309/1‑93 (MQ=255)

                                       ttttCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCggg                                                 >  2:2706505/1‑93 (MQ=255)

                                        tttCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGt                                                <  1:2770364/93‑1 (MQ=255)

                                        tttCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGt                                                <  1:247704/93‑1 (MQ=255)

                                          tCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTgg                                              >  1:1061366/1‑93 (MQ=255)

                                            cAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGc                                            >  2:3529616/1‑93 (MQ=255)

                                             aGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCt                                           >  1:3888036/1‑93 (MQ=255)

                                                aaCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaa                                        >  1:3082357/1‑93 (MQ=255)

                                                  cGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaa                                      >  2:1106510/1‑93 (MQ=255)

                                                   ggTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaaa                                     <  2:2098467/93‑1 (MQ=255)

                                                    gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAt                                    <  1:4282343/93‑1 (MQ=255)

                                                    gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAt                                    <  1:1225698/93‑1 (MQ=255)

                                                     ttGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATg                                   >  2:1848660/1‑93 (MQ=255)

                                                     ttGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATg                                   <  2:1638312/93‑1 (MQ=255)

                                                     ttGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATg                                   >  2:2914660/1‑93 (MQ=255)

                                                       gAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGGGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGcc                                 >  1:1286271/1‑93 (MQ=255)

                                                        aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCa                                <  1:236043/93‑1 (MQ=255)

                                                        aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCa                                <  2:505194/93‑1 (MQ=255)

                                                        aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCa                                >  1:3306128/1‑93 (MQ=255)

                                                        aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCa                                <  1:3577869/93‑1 (MQ=255)

                                                        aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTCATCGTCACGCCGGAACATCACTACTCGGTACCGGGGGGGCTGAGAAATGcct                                >  1:819242/1‑92 (MQ=255)

                                                         aGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCAt                               <  2:3906112/93‑1 (MQ=255)

                                                          gCTGTGGCGGAGCAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATc                              <  1:2059155/93‑1 (MQ=255)

                                                          gCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATc                              >  2:1348986/1‑93 (MQ=255)

                                                          gCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATc                              <  2:69440/93‑1 (MQ=255)

                                                          gCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATc                              >  1:2384671/1‑93 (MQ=255)

                                                            tctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGa                            >  1:3014746/1‑93 (MQ=255)

                                                            tctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGa                            >  1:2231468/1‑93 (MQ=255)

                                                             ctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTTAAAAATGCCATCGAc                           >  1:4716452/1‑93 (MQ=255)

                                                             ctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGAc                           <  1:2077157/93‑1 (MQ=255)

                                                              tGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACt                          <  2:4066258/93‑1 (MQ=255)

                                                               ggCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTg                         <  2:3354373/93‑1 (MQ=255)

                                                               ggCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTg                         <  1:402585/93‑1 (MQ=255)

                                                                  ggAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCt                      >  1:3778472/1‑93 (MQ=255)

                                                                  ggAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCt                      >  2:3639581/1‑93 (MQ=255)

                                                                   gAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCtt                     >  1:2682881/1‑93 (MQ=255)

                                                                     aCAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTc                   >  2:3170917/1‑93 (MQ=255)

                                                                      cAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTcc                  <  2:3104236/93‑1 (MQ=255)

                                                                      cAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTcc                  <  1:2985287/93‑1 (MQ=255)

                                                                       aGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGCGTGGTCTGAAAAATGCCATCGACTGGCTTTccc                 <  2:4652316/93‑1 (MQ=255)

                                                                       aGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTccc                 >  2:1292500/1‑93 (MQ=255)

                                                                        gATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCg                >  1:1659301/1‑93 (MQ=255)

                                                                        gATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCg                >  2:1434011/1‑93 (MQ=255)

                                                                          tCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGcc              >  1:1095318/1‑93 (MQ=255)

                                                                          tCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGcc              <  1:1554858/93‑1 (MQ=255)

                                                                          tCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGcc              >  1:3266606/1‑93 (MQ=255)

                                                                           ccGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCt             >  2:376373/1‑93 (MQ=255)

                                                                           ccGTCACGCAGACGGCGTGGTGATCGTCACGTCATAATATAACTATTGGGTACCGGGTGGGCTGAAAAATTCCATCGACTGGCTTTCCCGcca             <  2:2876957/93‑2 (MQ=255)

                                                                             gTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGc           >  1:917991/1‑93 (MQ=255)

                                                                              tCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGcc          <  2:3783385/93‑1 (MQ=255)

                                                                                 ggCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGa       >  1:629970/1‑93 (MQ=255)

                                                                                  gCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGAt      <  1:3336542/93‑1 (MQ=255)

                                                                                  gCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGAt      >  2:2242805/1‑93 (MQ=255)

                                                                                   cTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATc     <  1:1737592/93‑1 (MQ=255)

                                                                                    tGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCa    >  2:1050077/1‑93 (MQ=255)

                                                                                    tGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCa    >  1:1348081/1‑93 (MQ=255)

                                                                                      gcGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc  <  2:2024211/92‑1 (MQ=255)

                                                                                      aCGGCGTGGTGATCGTCACGCCGGAATATAACTGCTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc  <  2:3077885/93‑1 (MQ=255)

                                                                                      aCGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc  >  2:2958875/1‑93 (MQ=255)

                                                                                      aCGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc  <  1:3830097/93‑1 (MQ=255)

CATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAC  >  NC_000913/3892415‑3892593

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