BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,655,609	G→C	R449R (CGC→CGG)	hscA ←	HscA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,655,609	0	G	C	97.3%	665.6 / 5.7	188	R449R (CGC→CGG)	hscA	HscA
Reads supporting (aligned to +/- strand): ref base G (2/3); new base C (77/106); total (79/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

Only 100 of 193 total aligned reads displayed.

ACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGC  >  NC_000913/2655519‑2655698

aCCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGAAGCGCAAAACGCGCCAGTGACCg                                                                                         >  1:2848494/1‑93 (MQ=255)

   gTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCa                                                                                      <  2:1823447/93‑1 (MQ=255)

      ggCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTc                                                                                   >  2:406734/1‑93 (MQ=255)

       gCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTcc                                                                                  <  1:479233/93‑1 (MQ=255)

        cATCGACGTGGAACGTCACGCGAATCGGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCt                                                                                 <  1:860931/93‑1 (MQ=255)

         aTCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTg                                                                                >  2:4568496/1‑93 (MQ=255)

          tCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGc                                                                               >  2:497313/1‑93 (MQ=255)

             aCCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCAcc                                                                            <  2:5433423/93‑1 (MQ=255)

              ccTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCa                                                                           >  1:6075783/1‑93 (MQ=255)

              ccTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCa                                                                           >  1:1963571/1‑93 (MQ=255)

               cTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAg                                                                          >  1:971484/1‑93 (MQ=255)

               cTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAg                                                                          >  2:5820083/1‑93 (MQ=255)

                tGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGt                                                                         >  2:2038142/1‑93 (MQ=255)

                tGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGt                                                                         >  2:1739006/1‑93 (MQ=255)

                   accGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTcg                                                                      <  2:1057842/91‑1 (MQ=255)

                    aCGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTcgc                                                                     <  2:3630365/93‑1 (MQ=255)

                    aCGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTcgc                                                                     >  2:5405139/1‑93 (MQ=255)

                      gTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCgcgc                                                                   >  2:161420/1‑93 (MQ=255)

                      gTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCgcgc                                                                   >  2:4149393/1‑93 (MQ=255)

                       tCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCt                                                                  <  1:1028858/93‑1 (MQ=255)

                        cACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTc                                                                 >  2:3172517/1‑93 (MQ=255)

                           gcgAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCAcc                                                              >  2:4332124/1‑93 (MQ=255)

                            cgAATATGCGCCCCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCAccc                                                             <  1:2130298/93‑1 (MQ=255)

                            cgAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCAccc                                                             <  2:1499182/93‑1 (MQ=255)

                            cgAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCAccc                                                             <  2:4696898/93‑1 (MQ=255)

                             gAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTAGCGCTCACCCt                                                            >  2:3948719/1‑93 (MQ=255)

                              aaTATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTg                                                           <  2:1805173/93‑1 (MQ=255)

                                tatGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCa                                                         <  1:312680/93‑1 (MQ=255)

                                 atGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCAt                                                        <  1:3413869/93‑1 (MQ=255)

                                   gcgcACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTa                                                      >  2:4112149/1‑93 (MQ=255)

                                    cgcACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTAc                                                     <  1:5714731/93‑1 (MQ=255)

                                     gcACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACa                                                    >  1:4249214/1‑93 (MQ=255)

                                     gcACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACa                                                    >  2:4458830/1‑93 (MQ=255)

                                     gcACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACa                                                    <  2:515096/93‑1 (MQ=255)

                                      cACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACAt                                                   <  2:14177/93‑1 (MQ=255)

                                      cACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACAt                                                   <  1:3367011/93‑1 (MQ=255)

                                         cGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGa                                                >  2:4460474/1‑93 (MQ=255)

                                          gCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGAt                                               >  1:5528612/1‑93 (MQ=255)

                                           ccAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATa                                              <  1:729667/93‑1 (MQ=255)

                                           ccAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATa                                              <  1:1551570/93‑1 (MQ=255)

                                           caaGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATa                                              <  2:4846094/91‑1 (MQ=255)

                                              gCCGGTAGCGCCGGAATACCAAGCAGCGCAAAACGCGCCAGTGACAGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGAc                                           >  1:4219164/1‑93 (MQ=255)

                                               ccGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACa                                          <  2:2341855/93‑1 (MQ=255)

                                                cGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACAt                                         >  1:5895393/1‑93 (MQ=255)

                                                 ggTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATc                                        >  1:472778/1‑93 (MQ=255)

                                                  gTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCg                                       >  1:3371205/1‑93 (MQ=255)

                                                  gTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCg                                       >  1:4294898/1‑93 (MQ=255)

                                                    aGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGcc                                     <  1:3991183/93‑1 (MQ=255)

                                                    aGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGcc                                     <  1:2671711/93‑1 (MQ=255)

                                                       gcCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTc                                  <  1:799305/93‑1 (MQ=255)

                                                       gcCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTc                                  <  2:729552/93‑1 (MQ=255)

                                                       gcCGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTc                                  >  1:1600297/1‑93 (MQ=255)

                                                         cGGAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTg                                <  1:3581741/93‑1 (MQ=255)

                                                          ggAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGa                               >  1:3938455/1‑93 (MQ=255)

                                                          ggAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGa                               <  2:3735017/93‑1 (MQ=255)

                                                           gAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGAc                              >  1:1318071/1‑93 (MQ=255)

                                                           gAATACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGAc                              >  1:2904936/1‑93 (MQ=255)

                                                             aTACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCa                            <  1:3188130/93‑1 (MQ=255)

                                                              tACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCAt                           >  1:897097/1‑93 (MQ=255)

                                                              tACCACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCAt                           >  1:828411/1‑93 (MQ=255)

                                                              tACCACGCAGCGCAAAACGCGCCAGTGACCGGAAGTCCTGCACCAGTTCGCGGTCACACTGCATTACATGGATAGACATCGCCGTCTGACCAt                           >  2:261804/1‑93 (MQ=255)

                                                                ccACGCAGCGCAAAACGCGCCAGTGACCGGCCGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCt                         >  2:4224704/1‑93 (MQ=255)

                                                                ccACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCt                         <  2:392345/93‑1 (MQ=255)

                                                                 cACGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCtt                        <  2:3876798/93‑1 (MQ=255)

                                                                  aCGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCttt                       >  1:3796126/1‑93 (MQ=255)

                                                                  aCGCAGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCttt                       <  2:4187013/93‑1 (MQ=255)

                                                                      aGCGCAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAg                   <  2:1878621/93‑1 (MQ=255)

                                                                       gcgcAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAgg                  >  1:2065091/1‑93 (MQ=255)

                                                                       gcgcAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAgg                  >  1:5006954/1‑93 (MQ=255)

                                                                        cgcAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAggt                 <  1:3232163/93‑1 (MQ=255)

                                                                        cgcAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCACTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAggt                 <  2:2285726/93‑1 (MQ=255)

                                                                          cAAAACGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAggtgg               <  1:4652155/93‑1 (MQ=255)

                                                                           aaaaCGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAggtggt              <  2:2789446/93‑1 (MQ=255)

                                                                            aaaCGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGgtggtg             <  2:441122/93‑1 (MQ=255)

                                                                            aaaCGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGgtggtg             <  2:2211171/93‑1 (MQ=255)

                                                                             aaCGCGCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGa            >  2:5935039/1‑93 (MQ=255)

                                                                              aCGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGAGaa           <  1:4141141/93‑1 (MQ=255)

                                                                               cgcgCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGaaa          <  2:2179650/93‑1 (MQ=255)

                                                                               cgcgCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGaaa          <  2:4520914/93‑1 (MQ=255)

                                                                               cgcgCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGaaa          <  2:1781826/93‑1 (MQ=255)

                                                                               cgcgCCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGaaa          <  1:880044/93‑1 (MQ=255)

                                                                                 cgcCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATc        <  2:5155195/93‑1 (MQ=255)

                                                                                 cgcCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATc        <  1:1155646/93‑1 (MQ=255)

                                                                                 cgcCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATc        <  2:2366170/93‑1 (MQ=255)

                                                                                 cgcCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATc        <  1:3782382/93‑1 (MQ=255)

                                                                                 caccAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATc        <  1:795109/91‑1 (MQ=255)

                                                                                  gcCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATcc       <  2:2591367/93‑1 (MQ=255)

                                                                                  gcCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATcc       >  1:409742/1‑93 (MQ=255)

                                                                                  gcCAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATcc       <  2:3985073/93‑1 (MQ=255)

                                                                                   ccAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCt      <  1:2288260/93‑1 (MQ=255)

                                                                                    cAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTg     <  1:725371/93‑1 (MQ=255)

                                                                                    cAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTg     <  1:2712872/93‑1 (MQ=255)

                                                                                    cAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTg     >  2:5248432/1‑93 (MQ=255)

                                                                                    cAGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTg     >  1:4715968/1‑93 (MQ=255)

                                                                                     aGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGa    <  2:1130716/93‑1 (MQ=255)

                                                                                     aGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGa    >  2:1398221/1‑93 (MQ=255)

                                                                                     aGTGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGa    <  2:1966081/93‑1 (MQ=255)

                                                                                       tGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAgc  <  1:3206745/93‑1 (MQ=255)

                                                                                       tGACCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAgc  <  2:2893586/93‑1 (MQ=255)

                                                                                       tGACCGGCAGCCCTGCACCAGTTTGCGCTAACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAgc  >  1:4317150/1‑93 (MQ=255)

ACCGTCGGCATCGACCTGGAACGTCACGCGAATATGCGCACCGCCAGCCGGTAGCGCCGGAATACCACGCAGCGCAAAACGCGCCAGTGAGCGGCAGTCCTGCACCAGTTCGCGCTCACCCTGCATTACATGGATAGACATCGCCGTCTGACCATCTTTAAAGGTGGTGAAATCCTGAGC  >  NC_000913/2655519‑2655698

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