BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,682,619	A→G	R199R (CGT→CGC)	yhjK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,682,619	0	A	G	99.2%	1303.9 / ‑1.5	387	R199R (CGT→CGC)	yhjK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (170/214); total (171/216)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

Only 100 of 398 total aligned reads displayed.

TTTTCGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATCAG  >  NC_000913/3682527‑3682708

   tCGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGgcgc                                                                                        >  1:628606/1‑93 (MQ=255)

     gttctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCgg                                                                                      >  1:8756842/1‑93 (MQ=255)

     gttctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCgg                                                                                      <  1:3618939/93‑1 (MQ=255)

      ttctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCgga                                                                                     >  2:7613990/1‑92 (MQ=255)

      ttctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGt                                                                                     >  2:8422549/1‑93 (MQ=255)

      ttctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGt                                                                                     >  1:6787915/1‑93 (MQ=255)

      ttctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGt                                                                                     >  2:3432971/1‑93 (MQ=255)

      ttctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGt                                                                                     >  1:2661377/1‑93 (MQ=255)

      ttctgttctTCATAATGGCGCTGCAGCAATTGCTGCTTGAGGTTGTAACTGCGCACCAACATCCCGATTTCATCGTCCTGATGCAGGCGCGGt                                                                                     <  1:9643/93‑1 (MQ=255)

       tctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTa                                                                                    >  1:3224406/1‑93 (MQ=255)

       tctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTa                                                                                    >  1:5446923/1‑93 (MQ=255)

        ctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTaa                                                                                   <  1:6938661/93‑1 (MQ=255)

        ctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTaa                                                                                   <  2:801747/93‑1 (MQ=255)

         tgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTACCTGGGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAAt                                                                                  <  2:5437821/93‑1 (MQ=255)

         tgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAAt                                                                                  <  1:2218662/93‑1 (MQ=255)

           ttcttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGc                                                                                >  2:6316536/1‑93 (MQ=255)

             cttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCa                                                                              >  2:212967/1‑93 (MQ=255)

              ttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAg                                                                             >  1:2434550/1‑93 (MQ=255)

               tcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGt                                                                            <  1:3673143/93‑1 (MQ=255)

                 aTAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTg                                                                          <  1:7727440/93‑1 (MQ=255)

                 aTAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTg                                                                          <  1:3982514/93‑1 (MQ=255)

                 aTAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTg                                                                          <  2:6764849/93‑1 (MQ=255)

                 aTAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTg                                                                          <  1:8181043/93‑1 (MQ=255)

                  tAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCTTCCTGATGCAGGCGCGGTAATGCCAGTTgg                                                                         <  2:1647834/93‑1 (MQ=255)

                  tAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTgg                                                                         >  1:7292885/1‑93 (MQ=255)

                  tAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTgg                                                                         >  1:4195042/1‑93 (MQ=255)

                  tAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTgg                                                                         <  1:1335990/93‑1 (MQ=255)

                   aaTGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGt                                                                        <  1:5647624/93‑1 (MQ=255)

                    aTGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTg                                                                       >  2:5962316/1‑93 (MQ=255)

                     tGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTgc                                                                      >  2:1760279/1‑92 (MQ=255)

                      ggCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGAc                                                                     <  1:7386546/93‑1 (MQ=255)

                       gcgcTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGAcc                                                                    >  2:316091/1‑93 (MQ=255)

                         gcTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGGAATGCCAGTTGGTGACcaa                                                                  >  1:1203261/1‑93 (MQ=255)

                          cTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGCCcaac                                                                 >  2:2690549/1‑93 (MQ=255)

                           tGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACcaaca                                                                <  2:222618/93‑1 (MQ=255)

                            gcagcaATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACcaacaa                                                               >  1:3188556/1‑93 (MQ=255)

                            gcagcaATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACcaacaa                                                               <  2:2559010/93‑1 (MQ=255)

                            gcagcaATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACcaacaa                                                               >  2:6964888/1‑93 (MQ=255)

                            gcagcaATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACcaacaa                                                               <  1:8332034/93‑1 (MQ=255)

                             cagcaATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAg                                                              >  1:7621589/1‑93 (MQ=255)

                              agcaATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGc                                                             >  2:6575418/1‑93 (MQ=255)

                               gcaATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCt                                                            >  2:5708541/1‑93 (MQ=255)

                                caATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTc                                                           <  2:4978035/93‑1 (MQ=255)

                                  aTTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCt                                                         <  2:6097197/93‑1 (MQ=255)

                                  aTTGCTGGTTGAGGTTGTAAATGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCt                                                         >  2:7387803/1‑93 (MQ=255)

                                   ttGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCtt                                                        <  1:8559673/93‑1 (MQ=255)

                                    tGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTg                                                       <  2:7076160/93‑1 (MQ=255)

                                    tGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTg                                                       <  2:3517325/93‑1 (MQ=255)

                                      cTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGc                                                     >  1:1986977/1‑93 (MQ=255)

                                      cTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGc                                                     <  1:3377749/93‑1 (MQ=255)

                                        ggTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTg                                                   >  1:2900864/1‑93 (MQ=255)

                                         gTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCATGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTgg                                                  >  1:8451851/1‑93 (MQ=255)

                                         gTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTgg                                                  <  1:300800/93‑1 (MQ=255)

                                           tGAGGTTGTAACTGCGCTCCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGa                                                <  1:7810268/93‑1 (MQ=255)

                                           tGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTTATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGa                                                >  2:1917456/1‑93 (MQ=255)

                                           tGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGa                                                >  2:4685655/1‑93 (MQ=255)

                                           tGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGa                                                <  1:7507703/93‑1 (MQ=255)

                                            gAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGAt                                               <  1:5022952/93‑1 (MQ=255)

                                            gAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGAGACCAACAAGCTCCTTGGCTGGGAt                                               <  2:5159797/93‑1 (MQ=255)

                                             aGGTTGTAACTGCGCACAAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATg                                              >  1:7153930/1‑93 (MQ=255)

                                               gTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGc                                            >  1:4265953/1‑93 (MQ=255)

                                                ttGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCg                                           >  2:1135210/1‑93 (MQ=255)

                                                  gTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGCGACCAACAAGCTCCTTGGCTGGGATGGCGtt                                         >  1:5747297/1‑93 (MQ=255)

                                                   tAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTa                                        <  2:2450111/93‑1 (MQ=255)

                                                      cTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGt                                     <  2:935080/93‑1 (MQ=255)

                                                       tGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGtt                                    >  1:8627886/1‑93 (MQ=255)

                                                        gcgcACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTc                                   >  1:2303657/1‑93 (MQ=255)

                                                        gcgcACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTc                                   >  2:647458/1‑93 (MQ=255)

                                                             ccAACATCCCGATTTCATCGTCCTGATGCAGGCTCTGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAg                              <  1:6707474/93‑1 (MQ=255)

                                                               accATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCAAGTTGGTGACCAACAAGCTCCTTGGCTGCGATGGCGTTAAGCTCGCGAGCa                            >  2:7171218/3‑93 (MQ=255)

                                                                aCATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCaa                           <  1:465855/93‑1 (MQ=255)

                                                                 cATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAat                          <  2:3664794/93‑1 (MQ=255)

                                                                  aTACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAata                         >  1:611585/1‑93 (MQ=255)

                                                                  aTACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAata                         <  2:4807948/93‑1 (MQ=255)

                                                                    aCCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATAtt                       >  1:1494294/1‑93 (MQ=255)

                                                                    aCCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATAtt                       <  2:1296366/93‑1 (MQ=255)

                                                                     ccGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTg                      <  1:1011875/93‑1 (MQ=255)

                                                                      cGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGc                     <  1:6381015/93‑1 (MQ=255)

                                                                      cGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGc                     <  1:5579214/93‑1 (MQ=255)

                                                                       gATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCg                    >  2:3917741/1‑93 (MQ=255)

                                                                        aTTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGt                   >  2:5346324/1‑93 (MQ=255)

                                                                        aTTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGt                   <  2:1417085/93‑1 (MQ=255)

                                                                         tttCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTa                  <  2:1547157/93‑1 (MQ=255)

                                                                           tCATCGTCCTGATGCATGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGGTAAGTTCGCGAGCAATATTGCGTAAc                >  2:7807524/1‑93 (MQ=255)

                                                                           tCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAAc                <  1:5933476/93‑1 (MQ=255)

                                                                           tCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAAc                <  2:8288891/93‑1 (MQ=255)

                                                                             aTCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGCGCAATATTGCGTAACgg              >  1:4074445/1‑93 (MQ=255)

                                                                             aTCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACgg              >  1:4789688/1‑93 (MQ=255)

                                                                             aTCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACgg              <  1:7942602/93‑1 (MQ=255)

                                                                               cGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGAt            >  2:3005253/1‑93 (MQ=255)

                                                                               cGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGAt            <  1:4914672/93‑1 (MQ=255)

                                                                               cGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGAt            <  1:5320951/93‑1 (MQ=255)

                                                                                  ccTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCa         >  2:8726714/1‑93 (MQ=255)

                                                                                   cTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCaa        >  2:3080816/1‑93 (MQ=255)

                                                                                     gATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCaaaa      <  2:3296538/93‑1 (MQ=255)

                                                                                     gATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCaaaa      >  2:4339157/1‑93 (MQ=255)

                                                                                      aTGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAAt     >  2:4546079/1‑93 (MQ=255)

                                                                                       tGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATc    >  1:1381001/1‑93 (MQ=255)

                                                                                       tGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATc    <  2:4131598/93‑1 (MQ=255)

                                                                                        gCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATCa   <  2:2829427/93‑1 (MQ=255)

TTTTCGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATCAG  >  NC_000913/3682527‑3682708

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