BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,097,623	T→C	M162V (ATG→GTG)	rhaT ←	RhaT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,097,623	0	T	C	99.7%	1262.1 / ‑6.1	391	M162V (ATG→GTG)	rhaT	RhaT
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (222/168); total (223/168)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 400 total aligned reads displayed.

TCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCGCCTT  >  NC_000913/4097533‑4097712

tCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCt                                                                                         <  2:5951675/93‑1 (MQ=255)

 cGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCtt                                                                                        >  1:2798261/1‑93 (MQ=255)

 cGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCtt                                                                                        >  2:2451659/1‑93 (MQ=255)

  gCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTg                                                                                       <  1:9143297/93‑1 (MQ=255)

    aaaGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTgcg                                                                                     >  2:343060/1‑93 (MQ=255)

    aaaGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTgcg                                                                                     >  1:5411364/1‑93 (MQ=255)

    aaaGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTgcg                                                                                     >  1:52300/1‑93 (MQ=255)

     aaGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTgcgc                                                                                    >  1:7874541/1‑93 (MQ=255)

      aGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGct                                                                                   >  2:4537630/1‑93 (MQ=255)

      aGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGct                                                                                   <  2:8591863/93‑1 (MQ=255)

       ggACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGctc                                                                                  <  2:5367103/93‑1 (MQ=255)

        gccATCCCGGCAGAGAAAATGCCGCACCTCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGctct                                                                                 <  1:2540054/91‑1 (MQ=255)

        gACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGctct                                                                                 <  2:7981524/93‑1 (MQ=255)

         aCATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCtt                                                                                <  1:7588004/93‑1 (MQ=255)

         aCATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCtt                                                                                >  2:4824645/1‑93 (MQ=255)

           aTCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTc                                                                              <  1:8192303/93‑1 (MQ=255)

           aTCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTc                                                                              >  1:9449746/1‑93 (MQ=255)

           aTCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTc                                                                              <  2:950234/93‑1 (MQ=255)

              ccGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAAc                                                                           <  1:6487861/93‑1 (MQ=255)

               cGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACt                                                                          >  1:8893101/1‑93 (MQ=255)

                   agagAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGccc                                                                      >  2:8337187/1‑93 (MQ=255)

                   agagAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGccc                                                                      <  2:705101/93‑1 (MQ=255)

                    gagaAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCg                                                                     <  1:8013222/93‑1 (MQ=255)

                     agaAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCgg                                                                    <  2:7353162/93‑1 (MQ=255)

                     agaAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCgg                                                                    >  1:66200/1‑93 (MQ=255)

                     agaAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCgg                                                                    <  2:9172003/93‑1 (MQ=255)

                       aaaaTGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGgcg                                                                  <  1:1919233/93‑1 (MQ=255)

                          aTGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGa                                                               <  2:9430900/93‑1 (MQ=255)

                          aTGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGa                                                               <  2:7616573/93‑1 (MQ=255)

                          aTGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGa                                                               >  1:2289990/1‑93 (MQ=255)

                            gCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGt                                                             >  2:390896/1‑93 (MQ=255)

                               gCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTAc                                                          >  2:9037245/1‑93 (MQ=255)

                                cacaTCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACa                                                         >  2:373099/1‑93 (MQ=255)

                                cacaTCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACa                                                         >  2:6896250/1‑93 (MQ=255)

                                 acaTCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACac                                                        >  1:8653285/1‑92 (MQ=255)

                                   aTCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATc                                                      >  2:3073573/1‑93 (MQ=255)

                                    tCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATcc                                                     >  1:3606516/1‑93 (MQ=255)

                                    tCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATcc                                                     >  2:6625200/1‑93 (MQ=255)

                                      cccGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCt                                                   <  2:84872/92‑1 (MQ=255)

                                       ccGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTa                                                  >  2:1905731/1‑93 (MQ=255)

                                          ccAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGc                                               >  1:5135995/1‑93 (MQ=255)

                                          ccAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGc                                               <  1:4201791/93‑1 (MQ=255)

                                            aGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCa                                             >  1:4398004/1‑93 (MQ=255)

                                            aGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCa                                             <  1:5661434/93‑1 (MQ=255)

                                             gCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCaa                                            >  1:1875837/1‑93 (MQ=255)

                                             gCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCaa                                            >  2:4697153/1‑93 (MQ=255)

                                             gCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCaa                                            >  2:7481172/1‑93 (MQ=255)

                                             gCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCaa                                            >  2:1679956/1‑93 (MQ=255)

                                               aCCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATc                                          <  1:8860952/93‑1 (MQ=255)

                                                ccAGCCCTTTTTTCAGATTGAGCTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCACCTGCCCGGCGCGAGTTACAATCCCTACGCCAATCa                                         <  2:1495077/93‑1 (MQ=255)

                                                ccAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCa                                         >  2:9866771/1‑93 (MQ=255)

                                                ccAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCa                                         >  2:5730096/1‑93 (MQ=255)

                                                ccAGCCCTGTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCa                                         <  1:3136726/93‑1 (MQ=255)

                                                 cAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAg                                        >  1:7813149/1‑93 (MQ=255)

                                                   gCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAgcg                                      <  1:1339755/93‑1 (MQ=255)

                                                   gCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAgcg                                      <  2:816721/93‑1 (MQ=255)

                                                    cccTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAgcgc                                     >  2:9693725/1‑93 (MQ=255)

                                                      cTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGcca                                   <  2:1774026/93‑1 (MQ=255)

                                                      cTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGcca                                   <  2:6388331/93‑1 (MQ=255)

                                                       ttttttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCaa                                                                            >  2:1133105/1‑51 (MQ=255)

                                                       ttttttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccac                                  >  2:4414233/1‑93 (MQ=255)

                                                       ttttttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccac                                  >  2:5551443/1‑93 (MQ=255)

                                                        tttttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccacc                                 <  2:1338164/93‑1 (MQ=255)

                                                        tttttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccacc                                 >  2:5215339/1‑93 (MQ=255)

                                                        tttttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccacc                                 >  1:8523157/1‑93 (MQ=255)

                                                          tttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAg                               >  1:4779618/1‑93 (MQ=255)

                                                          tttCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAg                               >  1:6031435/1‑93 (MQ=255)

                                                           ttCAGATTGACCTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGa                              <  1:5906488/93‑1 (MQ=255)

                                                            tCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGaa                             <  2:7779967/93‑1 (MQ=255)

                                                            tCAGATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGaa                             <  1:6118897/93‑1 (MQ=255)

                                                               gATTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGc                          >  1:7191212/1‑93 (MQ=255)

                                                                aTTGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGcc                         <  1:9234065/93‑1 (MQ=255)

                                                                 ttGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCg                        >  2:5427220/1‑93 (MQ=255)

                                                                 ttGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCg                        >  2:6141982/1‑93 (MQ=255)

                                                                  tGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGa                       <  2:5854173/93‑1 (MQ=255)

                                                                  tGAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGa                       <  1:2879806/93‑1 (MQ=255)

                                                                   gAACTCTTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAg                      <  2:2753735/93‑1 (MQ=255)

                                                                      ctctTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCaa                   <  1:8432847/93‑1 (MQ=255)

                                                                      ctctTCGGCTTTAATGCCCACCTTGCACTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCaa                   <  1:7018498/93‑1 (MQ=255)

                                                                       tctTCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAAc                  >  2:1214721/1‑93 (MQ=255)

                                                                         ttCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGt                >  2:4218927/1‑93 (MQ=255)

                                                                          tCGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTc               >  1:6321409/1‑93 (MQ=255)

                                                                           cGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCa              <  1:5626219/93‑1 (MQ=255)

                                                                           cGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCa              <  2:2003740/93‑1 (MQ=255)

                                                                           cGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCa              <  1:4129472/93‑1 (MQ=255)

                                                                           cGGCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCa              <  1:10038769/93‑1 (MQ=255)

                                                                            ggCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCAt             >  1:9995019/1‑93 (MQ=255)

                                                                             gCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATg            >  1:4562778/1‑93 (MQ=255)

                                                                             gCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATg            >  1:3979743/1‑93 (MQ=255)

                                                                             gCTTTAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATg            <  1:3514877/93‑1 (MQ=255)

                                                                               tttAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCg          <  2:319263/93‑1 (MQ=255)

                                                                               tttAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCg          >  1:2064879/1‑93 (MQ=255)

                                                                                ttAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGt         <  2:2338547/93‑1 (MQ=255)

                                                                                ttAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGt         >  1:3824307/1‑93 (MQ=255)

                                                                                 tAATGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTc        >  1:9682067/1‑93 (MQ=255)

                                                                                  aaTGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAACCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTcc       >  1:9072167/1‑93 (MQ=255)

                                                                                   aTGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCg      >  2:6771426/1‑93 (MQ=255)

                                                                                   aTGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCg      >  1:3804624/1‑93 (MQ=255)

                                                                                    tGCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCGc     <  1:1170251/93‑1 (MQ=255)

                                                                                     gCCCACCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCGcc    <  1:6855564/93‑1 (MQ=255)

TCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCACCAGAACGCCGAGCAACGTCATGCGTCCGCCTT  >  NC_000913/4097533‑4097712

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