BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,472,649	G→A	R107R (CGC→CGT)	yheM ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,472,649	0	G	A	96.7%	816.7 / 5.1	242	R107R (CGC→CGT)	yheM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (5/3); new base A (129/105); total (134/108)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

Only 100 of 255 total aligned reads displayed.

GCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGCACAAA  >  NC_000913/3472559‑3472737

gcaAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCgcgacg                                                                                        <  2:1329636/93‑1 (MQ=255)

gcaAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCgcgacg                                                                                        >  2:2005915/1‑93 (MQ=255)

gcaAAATCCGTCAGCCAGGGCGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCgcgacg                                                                                        >  1:2499613/1‑93 (MQ=255)

  aaaaTCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTa                                                                                      <  1:1046007/93‑1 (MQ=255)

   aaaTCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTaa                                                                                     <  1:1568642/93‑1 (MQ=255)

   aaaTCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTaa                                                                                     <  2:1204845/93‑1 (MQ=255)

   aaaTCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTaa                                                                                     <  1:383289/93‑1 (MQ=255)

    aaTCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAg                                                                                    <  2:4229635/93‑1 (MQ=255)

     aTCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAgg                                                                                   >  2:2779286/1‑93 (MQ=255)

      tCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGc                                                                                  >  1:4237030/1‑93 (MQ=255)

      tCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGc                                                                                  <  2:4483159/93‑1 (MQ=255)

      tCCGTCAGCCAGGGAGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGc                                                                                  >  2:1830998/1‑93 (MQ=255)

        cGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCAt                                                                                >  2:2064139/1‑93 (MQ=255)

        cGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCAt                                                                                <  1:4263835/93‑1 (MQ=255)

         gTCAGCCCGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATc                                                                               <  2:1258190/93‑1 (MQ=255)

         gTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATc                                                                               <  2:2716220/93‑1 (MQ=255)

          tCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCt                                                                              >  2:3554171/1‑93 (MQ=255)

            aGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGc                                                                            <  1:3140910/93‑1 (MQ=255)

                aGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCg                                                                        <  1:3494026/93‑1 (MQ=255)

                aGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCg                                                                        <  1:4252847/93‑1 (MQ=255)

                 gggTTAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGa                                                                       >  1:1629509/1‑93 (MQ=255)

                 gggTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGa                                                                       >  2:638464/1‑93 (MQ=255)

                 gggTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGa                                                                       >  1:5468/1‑93 (MQ=255)

                     gAGCGATGTAATGTGTGCAGCATAAACGCCTAAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAgcgg                                                                   <  1:2176630/93‑1 (MQ=255)

                       gCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAgcggcg                                                                 <  1:4328282/93‑1 (MQ=255)

                        cGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGt                                                                <  1:3708751/93‑1 (MQ=255)

                         gATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTg                                                               >  1:2400307/1‑93 (MQ=255)

                            gTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCt                                                            <  2:3413342/93‑1 (MQ=255)

                              aaTGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTc                                                          >  1:2661690/1‑93 (MQ=255)

                                tgtgtgCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGa                                                        >  2:1776376/1‑93 (MQ=255)

                                tgtgtgCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGa                                                        >  2:327677/1‑93 (MQ=255)

                                tgtgtgCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCCGCGTGGCTTCgc                                                        >  1:3884630/1‑92 (MQ=255)

                                  tgtgCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGaca                                                      >  2:2617102/1‑93 (MQ=255)

                                  tgtgCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGaca                                                      <  2:1958000/93‑1 (MQ=255)

                                   gtgCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGacaa                                                     >  1:3639965/1‑93 (MQ=255)

                                   gtgCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGacaa                                                     >  1:594781/1‑93 (MQ=255)

                                       agcaTAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACaaa                                                 <  1:3590292/93‑1 (MQ=255)

                                       agcaTAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACaaa                                                 <  2:4507638/93‑1 (MQ=255)

                                       agcaTAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACaaa                                                 >  1:2925886/1‑93 (MQ=255)

                                        gcaTAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAg                                                >  2:2472153/1‑93 (MQ=255)

                                         caTAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAgg                                               <  1:4477066/93‑1 (MQ=255)

                                          aTAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAggg                                              >  1:420102/1‑93 (MQ=255)

                                          aTAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAggg                                              >  2:3993765/1‑93 (MQ=255)

                                          aTAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAggg                                              >  1:1441519/1‑93 (MQ=255)

                                          aTAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAggg                                              >  2:4583277/1‑93 (MQ=255)

                                          aTAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAggg                                              >  2:4118524/1‑93 (MQ=255)

                                           tAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGgt                                             <  2:544737/93‑1 (MQ=255)

                                           tAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGgt                                             >  1:4628139/1‑93 (MQ=255)

                                           tAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGgt                                             <  1:886322/93‑1 (MQ=255)

                                             aaCGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGgtgt                                           >  1:2831241/1‑93 (MQ=255)

                                              aCGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTc                                          <  1:4168741/93‑1 (MQ=255)

                                              aCGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTc                                          <  1:4585076/93‑1 (MQ=255)

                                              aCGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTc                                          <  1:2620567/93‑1 (MQ=255)

                                                gCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTCGCTTCGACAACAAAGGGTGTCTg                                        >  1:2279501/1‑93 (MQ=255)

                                                 cctcaaaacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGc                                       >  1:3207699/1‑93 (MQ=255)

                                                  ctcaaaacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCg                                      >  2:4164872/1‑93 (MQ=255)

                                                  ctcaaaacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCg                                      <  1:3093633/93‑1 (MQ=255)

                                                  ctcaaaacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCg                                      >  2:4659490/1‑93 (MQ=255)

                                                  ctcaaaacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCg                                      <  1:2225343/93‑1 (MQ=255)

                                                   tcaaaacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCgg                                     >  1:4117371/1‑93 (MQ=255)

                                                   tcaaaacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGAGTCTGCgg                                     >  2:4428024/1‑93 (MQ=255)

                                                    caaaacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGa                                    <  1:2572495/93‑1 (MQ=255)

                                                     aaaacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGAt                                   >  1:3322418/1‑93 (MQ=255)

                                                      aaacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATc                                  <  2:2952517/93‑1 (MQ=255)

                                                       aacctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCt                                 <  1:2728229/93‑1 (MQ=255)

                                                        acctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTa                                >  2:372419/1‑93 (MQ=255)

                                                          ctcaacaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAAc                              <  2:4136893/93‑1 (MQ=255)

                                                          ctcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAAc                              <  2:2356684/93‑1 (MQ=255)

                                                           tcaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAAcc                             <  2:3090707/93‑1 (MQ=255)

                                                           taaaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAAcc                             <  2:1023244/91‑1 (MQ=255)

                                                            caaaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAAccc                            >  2:1457117/1‑93 (MQ=255)

                                                              aaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCcgc                          <  1:1164744/93‑1 (MQ=255)

                                                              aaaTAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCcgc                          >  2:677817/1‑93 (MQ=255)

                                                                 tAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGtt                       >  2:111900/1‑93 (MQ=255)

                                                                 tAACATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGtt                       >  1:3738619/1‑93 (MQ=255)

                                                                   aCATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTcg                     >  2:2113060/1‑93 (MQ=255)

                                                                   aCATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTcg                     >  2:3847427/1‑93 (MQ=255)

                                                                   aCATCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTcg                     <  2:4323950/93‑1 (MQ=255)

                                                                     aTCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTcgcg                   >  1:3559788/1‑93 (MQ=255)

                                                                     aTCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGCCTGCGGCTCTAACCCGCGTTcgcg                   >  1:2347453/1‑93 (MQ=255)

                                                                      tCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCgcgc                  <  2:2634637/93‑1 (MQ=255)

                                                                      tCGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCgcgc                  <  1:2013290/93‑1 (MQ=255)

                                                                       cGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCa                 <  1:138499/93‑1 (MQ=255)

                                                                       cGTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCa                 <  1:1746181/93‑1 (MQ=255)

                                                                        gTAGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAg                >  2:800246/1‑93 (MQ=255)

                                                                          aGTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTg              <  1:1813256/93‑1 (MQ=255)

                                                                           gTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGa             >  1:23058/1‑93 (MQ=255)

                                                                           gTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGa             >  1:3949572/1‑93 (MQ=255)

                                                                           gTTGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGa             >  2:859096/1‑93 (MQ=255)

                                                                            ttGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGaa            <  2:932043/93‑1 (MQ=255)

                                                                            ttGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGaa            <  1:3283006/93‑1 (MQ=255)

                                                                            ttGGCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGaa            <  2:1164708/93‑1 (MQ=255)

                                                                              ggCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGc          <  1:2789768/93‑1 (MQ=255)

                                                                              agCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCCGTCGCGCAGTGAAGc          >  2:1861708/2‑93 (MQ=255)

                                                                               gCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGcc         <  2:1673138/93‑1 (MQ=255)

                                                                               gCGAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGcc         >  2:710197/1‑93 (MQ=255)

                                                                                 gAGTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGc       >  2:4626454/1‑93 (MQ=255)

                                                                                   gTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGcac     <  1:152468/93‑1 (MQ=255)

                                                                                   gTTCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGcac     >  2:2930759/1‑93 (MQ=255)

                                                                                    ttCGCGACGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGcaca    >  2:301000/1‑93 (MQ=255)

GCAAAATCCGTCAGCCAGGGTGAGCGATGTAATGTGTGCAGCATAAACGCCTCAAAACCTCAAAATAACATCGTAGTTGGCGAGTTCGCGGCGTAAGGCATCTGCTTCGAGCGGCGTGGCTTCGACAACAAAGGGTGTCTGCGGATCTAACCCGCGTTCGCGCAGTGAAGCCGCACAAA  >  NC_000913/3472559‑3472737

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