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.6%	888.8 / ‑5.2	271	H68H (CAT→CAC)	uspA	universal stress protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (160/110); total (160/111)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 282 total aligned reads displayed.

CGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGC  >  NC_000913/3637854‑3638033

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

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

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

  tAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCAc                                                                                       >  2:3870009/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

       gTAAACTACTCTGACCTATACACCGGGCTTATTGTTGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGAcc                                                                                  <  1:1175593/93‑1 (MQ=255)

       gTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGAcc                                                                                  >  1:2824386/1‑93 (MQ=255)

        tAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCg                                                                                 >  1:6631067/1‑93 (MQ=255)

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

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

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

           actactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAgg                                                                              >  2:2148895/1‑92 (MQ=255)

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

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

                tctGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTcc                                                                         <  2:6313868/93‑1 (MQ=255)

                tctGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTcc                                                                         <  1:4597169/93‑1 (MQ=255)

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

                   gACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACt                                                                      >  2:1075225/1‑93 (MQ=255)

                   gACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACt                                                                      >  1:692346/1‑93 (MQ=255)

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

                     ccTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTaa                                                                    >  1:3486809/1‑93 (MQ=255)

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

                      cTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAAt                                                                   >  1:4026601/1‑93 (MQ=255)

                      cTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAAt                                                                   <  1:3960469/93‑1 (MQ=255)

                      cTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAAt                                                                   <  1:6311773/93‑1 (MQ=255)

                      cTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAAt                                                                   >  1:2407817/1‑93 (MQ=255)

                        ataCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGc                                                                 <  1:5039617/93‑1 (MQ=255)

                         taCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATCCa                                                                >  2:3353227/1‑93 (MQ=255)

                          acacCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAg                                                               <  1:5882054/93‑1 (MQ=255)

                             ccGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCt                                                            >  2:4792604/1‑93 (MQ=255)

                              cGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTa                                                           >  1:5251350/1‑93 (MQ=255)

                              cGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTa                                                           >  1:3888770/1‑93 (MQ=255)

                                ggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAcc                                                         >  1:838189/1‑93 (MQ=255)

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

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

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

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

                                       tGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCAc                                                  <  2:4735047/93‑1 (MQ=255)

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

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

                                          tgtgAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGa                                               <  1:4387354/93‑1 (MQ=255)

                                          tgtgAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGa                                               >  2:1853098/1‑93 (MQ=255)

                                           gtgAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGaa                                              <  2:538765/93‑1 (MQ=255)

                                              aaTCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAcc                                           >  1:5748663/1‑93 (MQ=255)

                                              aaTCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAcc                                           >  2:6263049/1‑93 (MQ=255)

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

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

                                                tCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCt                                         >  1:6012064/1‑93 (MQ=255)

                                                  tGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGa                                       >  2:565924/1‑93 (MQ=255)

                                                  tGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGa                                       >  2:3828997/1‑93 (MQ=255)

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

                                                      tGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGACACCCTGagagg                                   >  1:899977/1‑93 (MQ=255)

                                                       gATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggc                                  <  1:3214197/93‑1 (MQ=255)

                                                         tatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcag                                >  2:4657421/1‑93 (MQ=255)

                                                          atGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcagc                               >  2:4860976/1‑93 (MQ=255)

                                                            gCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcagcgg                             <  1:777718/93‑1 (MQ=255)

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

                                                                 aaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGAcc                        >  1:2978635/1‑93 (MQ=255)

                                                                  actcATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCt                       <  2:3704949/90‑1 (MQ=255)

                                                                  aCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCt                       >  1:1027433/1‑93 (MQ=255)

                                                                  aCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCt                       >  2:917391/1‑93 (MQ=255)

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

                                                                  aCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCt                       >  2:6038679/1‑93 (MQ=255)

                                                                  aCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCt                       >  1:4935098/1‑93 (MQ=255)

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

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

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

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

                                                                    gCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTgg                     >  2:1194205/1‑93 (MQ=255)

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

                                                                      aTCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGc                   >  1:1097415/1‑93 (MQ=255)

                                                                       tctcTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGcc                  <  2:2593506/93‑1 (MQ=255)

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

                                                                        ctctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCa                 <  1:1991117/93‑1 (MQ=255)

                                                                        ctctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCa                 <  1:4838634/93‑1 (MQ=255)

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

                                                                         tctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAg                >  1:3375225/1‑93 (MQ=255)

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

                                                                           tGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGt              <  1:2766431/93‑1 (MQ=255)

                                                                             aaGAGACACATCACGCGCTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTc            >  2:6136706/1‑93 (MQ=255)

                                                                             aaGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTc            >  1:2687916/1‑93 (MQ=255)

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

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

                                                                                agaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTgg         >  1:5360395/1‑93 (MQ=255)

                                                                                agaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTgg         <  2:4966875/93‑1 (MQ=255)

                                                                                 gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt        >  1:3803014/1‑93 (MQ=255)

                                                                                 gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt        >  2:4126204/1‑93 (MQ=255)

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

                                                                                    acaTCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGa     <  2:1051392/93‑1 (MQ=255)

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

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

                                                                                      atcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATg   >  1:3549214/1‑93 (MQ=255)

CGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGC  >  NC_000913/3637854‑3638033

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