BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,637,944	T→C	H68H (CAT→CAC)	uspA →	universal stress protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,637,944	0	T	C	99.8%	1568.0 / ‑5.4	478	H68H (CAT→CAC)	uspA	universal stress protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (249/228); total (250/228)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

Only 100 of 509 total aligned reads displayed.

CGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGC  >  NC_000913/3637854‑3638033

cGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGc                                                                                         <  2:1186154/93‑1 (MQ=255)

 gTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCa                                                                                        >  1:9306335/1‑93 (MQ=255)

 gTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCa                                                                                        >  1:9911050/1‑93 (MQ=255)

  tAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCAc                                                                                       <  1:3922421/93‑1 (MQ=255)

   aGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACt                                                                                      >  2:6795285/1‑93 (MQ=255)

   aGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACt                                                                                      <  1:3671875/93‑1 (MQ=255)

    gATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTg                                                                                     >  1:6593962/1‑93 (MQ=255)

    gATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTg                                                                                     <  1:6774215/93‑1 (MQ=255)

     aTGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGa                                                                                    <  1:8731254/93‑1 (MQ=255)

      tGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGAc                                                                                   >  1:7265309/1‑93 (MQ=255)

      tGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGAc                                                                                   <  2:7455870/93‑1 (MQ=255)

         aaaCTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGa                                                                                <  1:8885060/93‑1 (MQ=255)

         aaaCTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGa                                                                                <  1:7935349/93‑1 (MQ=255)

          aaCTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGac                                                                               >  1:2699469/1‑92 (MQ=255)

          aaCTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAg                                                                               >  1:7603888/1‑93 (MQ=255)

           actactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGc                                                                              >  1:301512/1‑93 (MQ=255)

             tactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCtt                                                                            >  1:2087940/1‑93 (MQ=255)

              actCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCttt                                                                           >  1:4307317/1‑93 (MQ=255)

              actCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCttt                                                                           <  1:6432495/93‑1 (MQ=255)

               ctctGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTc                                                                          <  1:8341194/93‑1 (MQ=255)

                 ctGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCa                                                                        >  1:2952722/1‑93 (MQ=255)

                  tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCAc                                                                       >  1:4165913/1‑93 (MQ=255)

                  tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCAc                                                                       <  1:1328250/93‑1 (MQ=255)

                    aCCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTa                                                                     <  2:5400908/93‑1 (MQ=255)

                     ccTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTaa                                                                    <  1:3529608/93‑1 (MQ=255)

                       tataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATg                                                                  >  2:2944820/1‑93 (MQ=255)

                       tataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATg                                                                  >  2:4985684/1‑93 (MQ=255)

                       tataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATg                                                                  >  2:6362735/1‑93 (MQ=255)

                        ataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGc                                                                 >  2:7278145/1‑93 (MQ=255)

                         taCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCa                                                                >  1:8059645/1‑93 (MQ=255)

                         taCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCAGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCa                                                                <  2:4446263/93‑1 (MQ=255)

                              cGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTa                                                           <  1:222041/93‑1 (MQ=255)

                               gggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAc                                                          <  2:3445864/93‑1 (MQ=255)

                                ggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAcc                                                         <  1:9407330/93‑1 (MQ=255)

                                   ttATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCaa                                                      >  1:5198336/1‑93 (MQ=255)

                                   ttATTGATGTGAATCTGGGGGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCaa                                                      >  2:3744271/1‑93 (MQ=255)

                                    tATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAAt                                                     <  2:8506517/93‑1 (MQ=255)

                                    tATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAAt                                                     >  2:8102249/1‑93 (MQ=255)

                                    tATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAAt                                                     >  2:2343315/1‑93 (MQ=255)

                                    tATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAAt                                                     <  1:5842970/93‑1 (MQ=255)

                                     aTTGATGTGAATCTGGGTTATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATc                                                    >  2:6169327/1‑93 (MQ=255)

                                     aTTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAAta                                                    >  2:2740400/1‑92 (MQ=255)

                                     aTTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATc                                                    >  1:7007840/1‑93 (MQ=255)

                                      ttGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCa                                                   <  1:4993995/93‑1 (MQ=255)

                                      ttGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCa                                                   <  2:4208529/93‑1 (MQ=255)

                                        gATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACt                                                 >  1:2523981/1‑93 (MQ=255)

                                         aTGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTg                                                <  1:5923333/93‑1 (MQ=255)

                                            tgAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGaaa                                             >  2:8201547/1‑93 (MQ=255)

                                            tgAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCACGCTAACGAATGACTGcaa                                             >  1:737221/1‑93 (MQ=255)

                                               aTCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAccc                                          >  2:8742354/1‑93 (MQ=255)

                                               aTCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAccc                                          <  1:1232685/93‑1 (MQ=255)

                                                tCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCt                                         <  2:3343400/93‑1 (MQ=255)

                                                tCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGCCTACCCAATCACTGAAACCCt                                         <  2:6447400/93‑1 (MQ=255)

                                                 cTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTg                                        <  1:5548003/93‑1 (MQ=255)

                                                 cTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTg                                        >  2:8834197/1‑93 (MQ=255)

                                                 cTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTg                                        <  2:9210276/93‑1 (MQ=255)

                                                  tGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACT‑ATGCAGGCTACCCAATCACTGAAACCCTGag                                      <  1:641218/93‑1 (MQ=255)

                                                   tggTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGag                                      <  1:6227032/92‑1 (MQ=255)

                                                   gggTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGag                                      <  1:6974982/93‑1 (MQ=255)

                                                    ggTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagc                                     >  1:1467390/1‑93 (MQ=255)

                                                       gATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggc                                  >  2:5225837/1‑93 (MQ=255)

                                                           tGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcagcg                              >  1:1919579/1‑93 (MQ=255)

                                                           tGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcgacagcg                              >  1:6334527/1‑93 (MQ=255)

                                                             cAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcagcggc                            <  2:4817765/93‑1 (MQ=255)

                                                             cAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcagcggc                            >  1:4418768/1‑93 (MQ=255)

                                                               gAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGa                          >  2:6665569/1‑93 (MQ=255)

                                                               gAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGa                          <  2:7701379/93‑1 (MQ=255)

                                                                 aaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGAcc                        <  1:9031334/93‑1 (MQ=255)

                                                                  aCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCt                       <  1:1979829/93‑1 (MQ=255)

                                                                   cGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTg                      <  2:7135249/93‑1 (MQ=255)

                                                                   cGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTGTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTg                      <  1:5325490/93‑1 (MQ=255)

                                                                    gCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTgg                     <  1:4093023/93‑1 (MQ=255)

                                                                    gCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTgg                     <  1:4731899/93‑1 (MQ=255)

                                                                     cATCTCTGAAGAGACACATCACGCACTGACCGAGCGTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTggg                    >  2:2660221/1‑93 (MQ=255)

                                                                      ctctcTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGc                   <  2:6037982/92‑1 (MQ=255)

                                                                      aTCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGc                   <  2:1548952/93‑1 (MQ=255)

                                                                       tctcTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGcc                  >  1:6059192/1‑93 (MQ=255)

                                                                       tctcTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGcc                  >  1:10032121/1‑93 (MQ=255)

                                                                         tctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAg                <  2:3042102/93‑1 (MQ=255)

                                                                          ctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAgg               <  1:9798226/93‑1 (MQ=255)

                                                                           tGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGt              >  2:177070/1‑93 (MQ=255)

                                                                            gAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGtt             >  2:5619796/1‑93 (MQ=255)

                                                                            gAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGtt             >  2:9436194/1‑93 (MQ=255)

                                                                            gAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGtt             >  2:1646631/1‑93 (MQ=255)

                                                                             aaGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGGTc            >  2:9732855/1‑93 (MQ=255)

                                                                              agagACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCt           <  2:3867743/93‑1 (MQ=255)

                                                                              agagACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCt           >  2:2266694/1‑93 (MQ=255)

                                                                               gagaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTg          >  2:4626542/1‑93 (MQ=255)

                                                                               gagaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTg          >  2:7020866/1‑93 (MQ=255)

                                                                                agaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTCAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTgg         >  1:4877856/1‑93 (MQ=255)

                                                                                 gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt        <  2:3629133/93‑1 (MQ=255)

                                                                                 gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt        <  2:1960137/93‑1 (MQ=255)

                                                                                 gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt        <  1:6900384/93‑1 (MQ=255)

                                                                                  acacATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTc       >  1:5470887/1‑93 (MQ=255)

                                                                                  acacATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTc       <  1:2289696/93‑1 (MQ=255)

                                                                                  acacATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTc       >  2:9541696/1‑93 (MQ=255)

                                                                                   cacaTCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCg      >  1:4483644/1‑93 (MQ=255)

                                                                                    acaTCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGa     >  1:5704695/1‑93 (MQ=255)

                                                                                     catcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGAt    >  1:3747149/1‑93 (MQ=255)

                                                                                     catcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGAt    <  1:1100743/93‑1 (MQ=255)

CGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGC  >  NC_000913/3637854‑3638033

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