BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,935,727	C→T	P270S (CCA→TCA)	rbsK →	ribokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,935,727	0	C	T	98.7%	1329.1 / ‑1.4	374	P270S (CCA→TCA)	rbsK	ribokinase
Reads supporting (aligned to +/- strand): ref base C (1/4); new base T (189/180); total (190/184)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 382 total aligned reads displayed.

GTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAA  >  NC_000913/3935637‑3935816

gTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCa                                                                                         >  1:3739260/1‑93 (MQ=255)

 ttCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCAt                                                                                        <  1:5390543/93‑1 (MQ=255)

 ttCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCAt                                                                                        >  1:157057/1‑93 (MQ=255)

 ttCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCAt                                                                                        >  1:4464993/1‑93 (MQ=255)

 ttCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCAt                                                                                        >  1:2900738/1‑93 (MQ=255)

 ttCCTGGATTCCGGGTGCAGGCTGTCAATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCAt                                                                                        >  1:7013633/1‑93 (MQ=255)

  tCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCAtt                                                                                       <  1:1001528/93‑1 (MQ=255)

  tCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCAtt                                                                                       >  1:5630263/1‑93 (MQ=255)

   ccTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTg                                                                                      <  1:6669602/93‑1 (MQ=255)

    cTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGc                                                                                     <  2:2884553/93‑1 (MQ=255)

     tGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGcc                                                                                    <  2:5742475/93‑1 (MQ=255)

      ggATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCa                                                                                   <  2:2484481/93‑1 (MQ=255)

       gATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAACACTCATTGCCag                                                                                  >  1:3974851/1‑93 (MQ=255)

        aTTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCTCGGCATTGCTGGAAGAAAAATCATTGCCaga                                                                                 <  1:5754678/93‑1 (MQ=255)

            cGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCg                                                                             <  2:883687/93‑1 (MQ=255)

              ggTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGAt                                                                           >  2:6453664/1‑93 (MQ=255)

                tGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTc                                                                         <  2:4711284/93‑1 (MQ=255)

                tGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTc                                                                         >  2:6836091/1‑93 (MQ=255)

                  cAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGt                                                                       <  1:1365749/93‑1 (MQ=255)

                   aGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGtt                                                                      <  1:7692534/93‑1 (MQ=255)

                   aGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGtt                                                                      >  1:1641990/1‑93 (MQ=255)

                    ggCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGttt                                                                     >  1:4683434/1‑93 (MQ=255)

                    ggCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGttt                                                                     >  1:2650255/1‑93 (MQ=255)

                    ggCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTCAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGttt                                                                     >  2:4240790/1‑93 (MQ=255)

                       tGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGc                                                                  >  2:5219268/1‑93 (MQ=255)

                         tCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGccc                                                                <  2:5520795/93‑1 (MQ=255)

                          cGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCa                                                               <  2:4480191/93‑1 (MQ=255)

                             ttccATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCAtgc                                                            <  1:1261780/91‑1 (MQ=255)

                             tACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCAtgc                                                            >  1:2019330/1‑93 (MQ=255)

                             tACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCAtgc                                                            >  1:3530084/1‑93 (MQ=255)

                             tACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCAtgc                                                            >  1:1789861/1‑93 (MQ=255)

                               ccATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                                          <  2:5564698/93‑1 (MQ=255)

                               ccATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCAtgctg                                                          >  1:5877280/1‑93 (MQ=255)

                                cATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCAtgctgc                                                         >  1:2708688/1‑93 (MQ=255)

                                 aTTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGcc                                                        <  2:2789786/93‑1 (MQ=255)

                                 aTTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGcc                                                        >  2:5329471/1‑93 (MQ=255)

                                 aTTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGcc                                                        >  1:6975497/1‑93 (MQ=255)

                                 aTTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGcc                                                        >  1:3234915/1‑93 (MQ=255)

                                    gctgcCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCt                                                     <  2:115498/93‑1 (MQ=255)

                                    gctgcCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCt                                                     <  1:5264146/93‑1 (MQ=255)

                                     ctgcCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTg                                                    >  1:950024/1‑93 (MQ=255)

                                     ctgcCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTg                                                    >  2:777866/1‑93 (MQ=255)

                                      tgcCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGc                                                   >  2:7532128/1‑93 (MQ=255)

                                        ccGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGa                                                 >  2:7246699/1‑93 (MQ=255)

                                         cGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGAt                                                <  2:3555680/93‑1 (MQ=255)

                                         cGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGAt                                                >  1:2147564/1‑93 (MQ=255)

                                           gagaTACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTg                                              >  1:3032734/1‑93 (MQ=255)

                                              aTACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCg                                           >  2:7936665/1‑93 (MQ=255)

                                              aTACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCg                                           <  1:438763/93‑1 (MQ=255)

                                               tACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGt                                          <  1:559307/93‑1 (MQ=255)

                                                   tttAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAaca                                      >  2:4788399/1‑93 (MQ=255)

                                                   tttAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAaca                                      <  1:651989/93‑1 (MQ=255)

                                                   tttAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAaca                                      >  2:4188514/1‑93 (MQ=255)

                                                   tttAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAaca                                      <  2:7091039/93‑1 (MQ=255)

                                                   tttAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAaca                                      >  1:4623901/1‑93 (MQ=255)

                                                   tttAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAaca                                      <  2:16449/93‑1 (MQ=255)

                                                    ttAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAacac                                     >  2:1753651/1‑93 (MQ=255)

                                                     tAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACg                                    <  1:1751897/93‑1 (MQ=255)

                                                     tAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACg                                    >  1:3123696/1‑93 (MQ=255)

                                                      aaCGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGt                                   >  2:2399677/1‑93 (MQ=255)

                                                      aaCGGTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGt                                   <  1:2200178/93‑1 (MQ=255)

                                                         ggTGCGTTAATCACGGCATTGCTTGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaaa                                <  2:3786314/93‑1 (MQ=255)

                                                         ggTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaaa                                <  2:2613993/93‑1 (MQ=255)

                                                          gTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAg                               >  2:1040906/1‑93 (MQ=255)

                                                          gTGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAg                               >  1:868290/1‑93 (MQ=255)

                                                           tGCGTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTCCCGCTGCGATTGCCGTAACACGTAAAgg                              >  1:4944260/1‑93 (MQ=255)

                                                              gTTAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGgcgc                           >  1:5126796/1‑93 (MQ=255)

                                                               ttAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTACAGGCGca                          >  2:3733371/1‑93 (MQ=255)

                                                               ttAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGca                          >  1:2427986/1‑93 (MQ=255)

                                                                tAATCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGcac                         <  1:3347982/93‑1 (MQ=255)

                                                                  aTCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACaa                       >  2:2970168/1‑93 (MQ=255)

                                                                   tCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAc                      <  1:480510/93‑1 (MQ=255)

                                                                   tCACGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAc                      <  2:7595843/93‑1 (MQ=255)

                                                                     aCGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCt                    >  2:1835575/1‑93 (MQ=255)

                                                                      cGGCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCtt                   <  2:178679/93‑1 (MQ=255)

                                                                        gCATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTcc                 <  1:6042743/93‑1 (MQ=255)

                                                                         cATTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCg                <  2:321144/93‑1 (MQ=255)

                                                                          aTTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGt               >  1:6860161/1‑93 (MQ=255)

                                                                          aTTGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGt               >  2:6007734/1‑93 (MQ=255)

                                                                           ttGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTa              <  1:4838205/93‑1 (MQ=255)

                                                                           ttGCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTa              >  2:3656468/1‑93 (MQ=255)

                                                                             gCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTAcc            <  2:5021350/93‑1 (MQ=255)

                                                                             gCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTAcc            >  2:5431309/1‑93 (MQ=255)

                                                                             gCTGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTAcc            <  2:4080940/93‑1 (MQ=255)

                                                                               tGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGt          >  1:7482065/1‑93 (MQ=255)

                                                                               tGGAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGt          >  2:1333578/1‑93 (MQ=255)

                                                                                ggAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTg         >  2:1922166/1‑93 (MQ=255)

                                                                                ggAAGAAAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTg         <  2:6224267/93‑1 (MQ=255)

                                                                                 gaagaaAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTgg        >  1:7573796/1‑93 (MQ=255)

                                                                                  aagaaAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGc       >  2:5120993/1‑93 (MQ=255)

                                                                                   agaaAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCg      >  2:741905/1‑93 (MQ=255)

                                                                                   agaaAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCg      >  1:2347113/1‑93 (MQ=255)

                                                                                   agaaAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCg      <  2:2290162/93‑1 (MQ=255)

                                                                                    gaaAAATCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGt     >  2:5931676/1‑93 (MQ=255)

                                                                                      aaaaCCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGa   <  1:7274833/93‑1 (MQ=255)

                                                                                       aaaTCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGaa  <  1:6766442/93‑1 (MQ=255)

                                                                                       aaaTCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGaa  <  1:6346593/93‑1 (MQ=255)

                                                                                       aaaTCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGaa  >  2:6332116/1‑93 (MQ=255)

                                                                                       aaaTCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGaa  >  1:5475998/1‑93 (MQ=255)

                                                                                       aaaTCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGaa  >  1:2515843/1‑93 (MQ=255)

GTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAA  >  NC_000913/3935637‑3935816

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