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.5%	605.7 / ‑4.9	185	H68H (CAT→CAC)	uspA	universal stress protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (93/91); total (94/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 195 total aligned reads displayed.

GTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGC  >  NC_000913/3637855‑3638033

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

 tAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCAc                                                                                       >  1:3002671/1‑93 (MQ=255)

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

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

  aGATGTAATCTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACt                                                                                      <  2:1062645/93‑1 (MQ=255)

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

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

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

       tAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCg                                                                                 <  2:4461389/93‑1 (MQ=255)

        aaaCTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGa                                                                                >  2:4452764/1‑93 (MQ=255)

        aaaCTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGa                                                                                >  2:3290662/1‑93 (MQ=255)

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

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

          actactCTGACCTAGACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGc                                                                              <  2:2771531/93‑1 (MQ=255)

           ctactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCt                                                                             <  2:201137/93‑1 (MQ=255)

           ctactCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCt                                                                             >  1:4337091/1‑93 (MQ=255)

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

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

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

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

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

                 tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACCGACCGAGCTTTCCAc                                                                       <  1:3441371/93‑1 (MQ=255)

                  gACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACt                                                                      <  1:785521/93‑1 (MQ=255)

                  gACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACt                                                                      <  1:2036988/93‑1 (MQ=255)

                   aCCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTa                                                                     >  1:287585/1‑93 (MQ=255)

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

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

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

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

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

                        taCACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCa                                                                <  2:3932544/93‑1 (MQ=255)

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

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

                           acCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGc                                                             <  1:2762832/93‑1 (MQ=255)

                           acCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGc                                                             <  1:4092586/93‑1 (MQ=255)

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

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

                              tgtcTTATTGATTGGAATCTGGGGGATATGCAGAAAAGCTTCTCCGAAGAGACACATCACGCACCGACCGAGCCTTCCCCTAATGCAGGCTAc                                                          >  2:4199767/4‑93 (MQ=255)

                              gggCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTAc                                                          >  2:1255659/1‑93 (MQ=255)

                                 cTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCa                                                       >  1:4639133/1‑93 (MQ=255)

                                 cTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCa                                                       >  1:2525470/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                 tGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGa                                       <  2:1428445/93‑1 (MQ=255)

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

                                                   ggTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagc                                     <  2:3432425/93‑1 (MQ=255)

                                                   ggTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagc                                     <  1:4153005/93‑1 (MQ=255)

                                                    gTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcg                                    <  1:9305/93‑1 (MQ=255)

                                                    gTGATATGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcg                                    <  2:779821/93‑1 (MQ=255)

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

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

                                                       atatGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggca                                 <  2:2986604/93‑1 (MQ=255)

                                                         atGCAGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcagc                               <  2:183343/93‑1 (MQ=255)

                                                           gCAGAAACGCATCTCTGAAGAGACACATCACACACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcagcgg                             <  1:3881983/93‑1 (MQ=255)

                                                             aGAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAgcggcggcggcg                           >  1:4245958/1‑93 (MQ=255)

                                                              gAAACGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAgcggcgc                          >  2:3977763/1‑92 (MQ=255)

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

                                                              aaaaCGCATCTCAGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGa                          <  1:2638642/92‑1 (MQ=255)

                                                               aaaCGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGAc                         >  1:4698927/1‑93 (MQ=255)

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

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

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

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

                                                                  cGCATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTg                      >  1:3045406/1‑93 (MQ=255)

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

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

                                                                    cATCTCTGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTggg                    >  1:2881314/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                                                                         ctGAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAgg               >  2:2857180/1‑93 (MQ=255)

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

                                                                           gAAGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGtt             <  1:4068972/93‑1 (MQ=255)

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

                                                                            aaGAGACACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTc            <  2:4701094/93‑1 (MQ=255)

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

                                                                             agagACACATCACGCACTGACCGAGATTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCt           <  2:283397/93‑1 (MQ=255)

                                                                              gagaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTg          <  2:3771479/93‑1 (MQ=255)

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

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

                                                                                gaCACATCACGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTATggc        >  2:1192487/1‑92 (MQ=255)

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

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

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

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

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

                                                                                     atcacGCACTTACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATg   <  2:4621220/93‑1 (MQ=255)

                                                                                      tcacGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGc  >  2:1578210/1‑93 (MQ=255)

GTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGC  >  NC_000913/3637855‑3638033

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