BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	52,760	A→G	G219G (GGT→GGC)	pdxA ←	PdxA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	52,760	0	A	G	99.3%	473.1 / ‑4.8	149	G219G (GGT→GGC)	pdxA	PdxA
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (94/54); total (95/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 157 total aligned reads displayed.

AGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGT  >  NC_000913/52671‑52848

aGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGccc                                                                                       >  1:2657566/1‑93 (MQ=255)

aGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGccc                                                                                       >  1:4445285/1‑93 (MQ=255)

 gggTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCa                                                                                      <  1:356274/93‑1 (MQ=255)

 gggTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCa                                                                                      >  2:1684445/1‑93 (MQ=255)

 gggTATCGGCAGGCAGCGGCCCGTTGAGGTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCCTGCCCa                                                                                      >  2:2621151/1‑93 (MQ=255)

   gTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCata                                                                                    >  2:1175187/1‑93 (MQ=255)

    tATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGGATAATGGTGTCTATCTCTTTCGTGCCCatat                                                                                   >  1:1060018/1‑93 (MQ=255)

     aTCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATg                                                                                  >  1:4453696/1‑93 (MQ=255)

       cGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGAc                                                                                <  2:3796466/93‑1 (MQ=255)

         gcaggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCg                                                                              >  1:3910904/1‑93 (MQ=255)

         gcaggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCg                                                                              >  1:3910902/1‑93 (MQ=255)

           aggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGcc                                                                            >  1:1442913/1‑93 (MQ=255)

            ggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCt                                                                           <  2:1932042/93‑1 (MQ=255)

            ggcagAGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCt                                                                           >  2:2002322/1‑93 (MQ=255)

             gcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCtt                                                                          <  1:4201492/93‑1 (MQ=255)

             gcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCtt                                                                          <  1:1184910/93‑1 (MQ=255)

              cagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTc                                                                         <  1:4072519/93‑1 (MQ=255)

              cagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTc                                                                         >  1:874711/1‑93 (MQ=255)

                gCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGc                                                                       >  1:1396120/1‑93 (MQ=255)

                gCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGc                                                                       >  1:3568377/1‑93 (MQ=255)

                 cGGCCCGTTGAGTTTCATCACCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGcc                                                                      <  2:2497830/93‑1 (MQ=255)

                  ggCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGccc                                                                     <  2:370892/93‑1 (MQ=255)

                     ccGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGACTAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCcgcg                                                                  >  1:105969/1‑93 (MQ=255)

                     ccGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCcgcg                                                                  >  2:3542712/1‑93 (MQ=255)

                      cGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCcgcgg                                                                 >  2:1516103/1‑92 (MQ=255)

                      cGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGt                                                                 >  2:3063585/1‑93 (MQ=255)

                      cGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGt                                                                 >  1:4322347/1‑93 (MQ=255)

                      cGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGt                                                                 >  2:3997916/1‑93 (MQ=255)

                      cGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGt                                                                 <  2:3276047/93‑1 (MQ=255)

                       gTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTg                                                                >  2:147986/1‑93 (MQ=255)

                          gAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCgg                                                             <  2:2695624/93‑1 (MQ=255)

                            gTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGAt                                                           >  1:3206355/1‑93 (MQ=255)

                            gTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGAt                                                           >  1:532936/1‑93 (MQ=255)

                            gTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGAt                                                           <  2:651173/93‑1 (MQ=255)

                              ttCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTc                                                         <  1:2309285/93‑1 (MQ=255)

                               tCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGCTTCTTTCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCa                                                        >  2:1067194/1‑93 (MQ=255)

                                cATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAg                                                       <  1:3398328/93‑1 (MQ=255)

                                 aTCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGc                                                      >  2:177219/1‑93 (MQ=255)

                                  tCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGcc                                                     <  1:557746/93‑1 (MQ=255)

                                  tCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGcc                                                     <  2:3959303/93‑1 (MQ=255)

                                   ccccTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGccc                                                    <  1:2898415/93‑1 (MQ=255)

                                      cTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCa                                                 <  1:2591218/93‑1 (MQ=255)

                                      cTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCa                                                 <  1:2706792/93‑1 (MQ=255)

                                      cTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCa                                                 >  2:2340457/1‑93 (MQ=255)

                                      cTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGCGCGGATTCAGCCCGCa                                                 >  2:659963/1‑93 (MQ=255)

                                          ggccGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGAcc                                             <  2:1585642/91‑1 (MQ=255)

                                          gcCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGAcc                                             <  2:1159743/93‑1 (MQ=255)

                                           cccGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCa                                            >  1:3973940/1‑93 (MQ=255)

                                           cccGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCa                                            >  2:3609516/1‑93 (MQ=255)

                                           cccGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCa                                            >  1:3132201/1‑93 (MQ=255)

                                            ccGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAg                                           >  1:4377106/1‑93 (MQ=255)

                                            ccGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAg                                           >  1:159576/1‑93 (MQ=255)

                                             cGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGa                                          <  1:1216319/93‑1 (MQ=255)

                                                aGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGAGCAGAATg                                       >  2:961339/1‑93 (MQ=255)

                                                aGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATg                                       <  2:3312455/93‑1 (MQ=255)

                                                aGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATg                                       <  2:1930017/93‑1 (MQ=255)

                                                 gCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAAtcc                                      >  2:4350948/1‑91 (MQ=255)

                                                 gCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgc                                      >  2:1395839/1‑93 (MQ=255)

                                                 gCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTAAGCCCGCAGACCAGAATgc                                      >  2:1322129/1‑93 (MQ=255)

                                                   tCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgcgc                                    >  2:4105712/1‑93 (MQ=255)

                                                   tCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgcgc                                    >  1:4213855/1‑93 (MQ=255)

                                                    cATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGcgcg                                   <  2:420518/93‑1 (MQ=255)

                                                      ttGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCTCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGt                                 <  1:2724638/93‑1 (MQ=255)

                                                      ttGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGt                                 >  1:2897972/1‑93 (MQ=255)

                                                      ttGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGt                                 >  1:4702456/1‑93 (MQ=255)

                                                          gCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCgg                             >  2:1702803/1‑93 (MQ=255)

                                                           cACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCGTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGc                            <  1:591150/93‑1 (MQ=255)

                                                             ccGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCaa                          <  2:2297389/93‑1 (MQ=255)

                                                             ccGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCaa                          <  1:259174/93‑1 (MQ=255)

                                                             ccGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCACCCGCGTGCGGATTCAGCTCGCAGACCAGAATGCGCGGTTCGGCaa                          >  1:2223248/1‑93 (MQ=255)

                                                               ggAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATa                        >  1:2389536/1‑93 (MQ=255)

                                                                gAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATAc                       >  2:1815481/1‑93 (MQ=255)

                                                                gAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATAc                       >  2:2801985/1‑93 (MQ=255)

                                                                gAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATAc                       >  2:3013784/1‑93 (MQ=255)

                                                                gAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATAc                       >  1:157127/1‑93 (MQ=255)

                                                                 aataatGGTGTCTATCTCTTCCGTGCCCATATGGCCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATAcc                      <  2:244025/93‑1 (MQ=255)

                                                                 aataatGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATAcc                      <  1:58241/93‑1 (MQ=255)

                                                                  ataatGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCATACCAGAATGCGCGGTTCGGCAATACCa                     >  2:2572503/1‑93 (MQ=255)

                                                                  ataatGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCa                     >  2:4387709/1‑93 (MQ=255)

                                                                   taatGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCaa                    >  1:3444946/1‑93 (MQ=255)

                                                                     atGTTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAAt                  <  2:1975362/93‑1 (MQ=255)

                                                                     atGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAAt                  >  1:784449/1‑93 (MQ=255)

                                                                     atGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAAt                  >  1:1519190/1‑93 (MQ=255)

                                                                      tGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAAtt                 >  2:2710572/1‑93 (MQ=255)

                                                                      tGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAAtt                 <  1:3219663/93‑1 (MQ=255)

                                                                       ggTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAAttt                <  1:2070694/93‑1 (MQ=255)

                                                                         tgtCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTgg              >  2:3660657/1‑93 (MQ=255)

                                                                         tgtCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTgg              <  1:18718/93‑1 (MQ=255)

                                                                           tCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGCTTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTAGTc            >  1:577280/1‑93 (MQ=255)

                                                                             tATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCg          >  2:4490509/1‑93 (MQ=255)

                                                                              aTCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGc         >  1:2513848/1‑93 (MQ=255)

                                                                              aTCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGc         >  2:4466038/1‑93 (MQ=255)

                                                                               tctcTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCa        >  1:3998523/1‑93 (MQ=255)

                                                                                ctctTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCaa       >  1:1856252/1‑93 (MQ=255)

                                                                                 tctTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCaaa      >  2:2648774/1‑93 (MQ=255)

                                                                                 tctTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCaaa      >  1:686541/1‑93 (MQ=255)

                                                                                  ctTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAAt     <  2:3257275/93‑1 (MQ=255)

                                                                                   ttCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCCAATc    >  1:3331937/1‑93 (MQ=255)

                                                                                   ttCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATc    >  1:3067905/1‑93 (MQ=255)

                                                                                     ccGTGCCCCTACGGCCGCCTCCGCCCGCGTGCGGATTCAGCCCGCAGCCCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGt  <  2:1248722/93‑1 (MQ=255)

AGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGT  >  NC_000913/52671‑52848

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: `‑`