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	99.8%	1346.8 / NA	422	G75G (GGT→GGC)	yieF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (205/215); total (206/216)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 453 total aligned reads displayed.

CATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAC  >  NC_000913/3892415‑3892593

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

cATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGgcgt                                                                                        >  2:7702324/1‑93 (MQ=255)

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

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

   tCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtggt                                                                                     <  2:184373/93‑1 (MQ=255)

    cccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtggtg                                                                                    <  1:8852278/93‑1 (MQ=255)

     ccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGa                                                                                   >  1:2155311/1‑93 (MQ=255)

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

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

        tGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCg                                                                                >  1:3380337/1‑93 (MQ=255)

          tATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTc                                                                              <  1:6188724/93‑1 (MQ=255)

           aTGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCa                                                                             <  2:7291532/93‑1 (MQ=255)

            tGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCAc                                                                            <  2:9736673/93‑1 (MQ=255)

            tGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCAc                                                                            <  2:5029083/93‑1 (MQ=255)

            tGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCAc                                                                            <  2:3219284/93‑1 (MQ=255)

             gACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACg                                                                           >  1:9947073/1‑93 (MQ=255)

              aCGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGc                                                                          >  1:8148445/1‑93 (MQ=255)

                gCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCg                                                                        >  1:3280610/1‑93 (MQ=255)

                 cTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCgg                                                                       >  2:2149077/1‑93 (MQ=255)

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

                   gACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGaa                                                                     <  1:4782056/93‑1 (MQ=255)

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

                     cGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAata                                                                   >  2:8586902/1‑93 (MQ=255)

                      gTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAGCAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAatat                                                                  <  1:2496731/93‑1 (MQ=255)

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

                       tACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAAtata                                                                 >  2:9394/1‑93 (MQ=255)

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

                          agcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAact                                                              >  1:5568547/1‑93 (MQ=255)

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

                           gcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAacta                                                             >  1:1083537/1‑93 (MQ=255)

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

                             agGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAacttct                                                           <  2:6099591/93‑1 (MQ=255)

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

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

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

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

                               gaagaaGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCg                                                         <  1:6970289/93‑1 (MQ=255)

                               gaagaaGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCg                                                         >  2:10000776/1‑93 (MQ=255)

                               gaagaaGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCg                                                         <  1:4942116/93‑1 (MQ=255)

                               gaagaaGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCg                                                         <  1:4681424/93‑1 (MQ=255)

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

                                 agaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCgcg                                                       >  1:3924128/1‑91 (MQ=255)

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

                                  gaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATAGTCACGCCGGAATATAACTACTCGGTa                                                      <  2:338630/93‑1 (MQ=255)

                                   aagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAc                                                     >  2:4190361/1‑93 (MQ=255)

                                   aagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAc                                                     >  1:4447867/1‑93 (MQ=255)

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

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

                                      gTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCgg                                                  <  1:3570376/93‑1 (MQ=255)

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

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

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

                                           ccAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTggg                                             >  2:6014912/1‑93 (MQ=255)

                                           ccAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTggg                                             >  2:441398/1‑93 (MQ=255)

                                            cAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGc                                            <  2:3988986/93‑1 (MQ=255)

                                             aGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCt                                           <  2:6498421/93‑1 (MQ=255)

                                             aGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCt                                           <  1:6850356/93‑1 (MQ=255)

                                              gCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTg                                          <  1:7571044/93‑1 (MQ=255)

                                              gCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTg                                          <  1:1243054/93‑1 (MQ=255)

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

                                                  cGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaa                                      <  2:2618004/93‑1 (MQ=255)

                                                   ggTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaaa                                     >  1:4853385/1‑93 (MQ=255)

                                                    gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAt                                    >  1:9171680/1‑93 (MQ=255)

                                                     ttGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTCCCGGGTGGGCTGAAAAATg                                   >  1:3193270/1‑93 (MQ=255)

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

                                                      tGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAtcc                                  <  2:2522296/93‑3 (MQ=255)

                                                      tGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGc                                  <  2:814651/93‑1 (MQ=255)

                                                       gAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGcc                                 >  1:7172168/1‑93 (MQ=255)

                                                       gAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGcc                                 >  1:4318268/1‑93 (MQ=255)

                                                       gAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGcc                                 >  2:4829730/1‑93 (MQ=255)

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

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

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

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

                                                               ggCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTg                         >  2:6564036/1‑93 (MQ=255)

                                                                gCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGAACCGGGTGGGCTGAAAAATGCCATCGACTgg                        >  2:4684864/1‑93 (MQ=255)

                                                                 cGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGc                       <  2:5295450/93‑1 (MQ=255)

                                                                 cGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGc                       <  1:9056893/93‑1 (MQ=255)

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

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

                                                                   gatcAGATCCGCCAGGCGGACGGCGTGGTGATCGTCACGCCGGAAGATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACAGGCtt                     >  2:8967119/4‑93 (MQ=255)

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

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

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

                                                                    aaCAGATCCGTCAGGCTGACGGCGTGGCGATAGCCAGGCCGGAATATAACTACTCGGTACCGCAGGGGCTGCAAAATGACATCGACTAGCttt                    >  2:877563/1‑93 (MQ=255)

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

                                                                     aCAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTc                   <  2:1356370/93‑1 (MQ=255)

                                                                      cAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTcc                  >  2:5168760/1‑93 (MQ=255)

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

                                                                        aaTCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCg                <  1:7779084/92‑1 (MQ=255)

                                                                         aTCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGc               >  1:9295979/1‑93 (MQ=255)

                                                                         aTCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGc               <  2:801339/93‑1 (MQ=255)

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

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

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

                                                                                   cTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATc     >  1:2546974/1‑93 (MQ=255)

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

                                                                                     gACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCaa   >  2:5614628/1‑93 (MQ=255)

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

                                                                                      aCGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc  <  2:3715116/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` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`