BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	112,238	A→G	V121A (GTG→GCG)	yacF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	112,238	0	A	G	96.5%	430.1 / 3.9	143	V121A (GTG→GCG)	yacF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (4/1); new base G (72/65); total (77/66)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

Only 100 of 145 total aligned reads displayed.

CGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATT  >  NC_000913/112149‑112327

cgcgCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCa                                                                                        <  1:3316874/93‑1 (MQ=255)

cgcgCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCa                                                                                        <  1:2891604/93‑1 (MQ=255)

 gcgcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCag                                                                                       >  2:3565869/1‑93 (MQ=255)

  cgcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCaga                                                                                      <  1:641973/93‑1 (MQ=255)

  cgcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCaga                                                                                      >  1:1848761/1‑93 (MQ=255)

   gcTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCagag                                                                                     >  2:1101583/1‑93 (MQ=255)

    cTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGc                                                                                    <  2:796409/93‑1 (MQ=255)

    cTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGc                                                                                    >  1:1803233/1‑93 (MQ=255)

     tGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCa                                                                                   >  2:4657702/1‑93 (MQ=255)

      gcgcCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCaa                                                                                  >  2:2420980/1‑93 (MQ=255)

       cgcCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAAt                                                                                 <  2:3812398/93‑1 (MQ=255)

       cgcCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAAt                                                                                 <  1:3976652/93‑1 (MQ=255)

          cTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCaa                                                                              >  2:3896017/1‑93 (MQ=255)

          cTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCaa                                                                              >  1:4677251/1‑93 (MQ=255)

           tGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCaaa                                                                             >  1:3701006/1‑93 (MQ=255)

            ggggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAAc                                                                            <  2:1735292/93‑1 (MQ=255)

            ggggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAAc                                                                            <  2:3068845/93‑1 (MQ=255)

             gggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACg                                                                           <  2:2125020/93‑1 (MQ=255)

              ggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGa                                                                          <  1:3972558/93‑1 (MQ=255)

              ggTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGa                                                                          >  1:550043/1‑93 (MQ=255)

                  gATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCtt                                                                      <  1:4192901/93‑1 (MQ=255)

                   aTGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTc                                                                     <  1:827015/93‑1 (MQ=255)

                     gCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCAc                                                                   <  2:1314268/93‑1 (MQ=255)

                      cAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACg                                                                  >  1:4674311/1‑93 (MQ=255)

                        gCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCa                                                                >  1:2343090/1‑93 (MQ=255)

                         ccAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAg                                                               >  2:3918329/1‑93 (MQ=255)

                          cAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGa                                                              >  2:749526/1‑93 (MQ=255)

                          cAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGa                                                              >  1:3992234/1‑93 (MQ=255)

                          cAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGa                                                              >  1:446988/1‑93 (MQ=255)

                            aaTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGaaa                                                            >  2:3989767/1‑93 (MQ=255)

                            aaTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGaaa                                                            <  2:41706/93‑1 (MQ=255)

                             aTGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAAt                                                           <  2:4073445/93‑1 (MQ=255)

                              tgtgCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAAtt                                                          <  2:397073/93‑1 (MQ=255)

                               gtgCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTg                                                         >  2:4746785/1‑93 (MQ=255)

                                 gCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGcc                                                       >  2:270175/1‑93 (MQ=255)

                                 gCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGcc                                                       <  2:2926772/93‑1 (MQ=255)

                                 gCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGcc                                                       >  1:855800/1‑93 (MQ=255)

                                    aTGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGa                                                    >  2:4625832/1‑93 (MQ=255)

                                    aTGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGa                                                    >  1:543894/1‑93 (MQ=255)

                                    aTGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGa                                                    <  1:3653717/93‑1 (MQ=255)

                                     tGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGAt                                                   >  1:3641466/1‑93 (MQ=255)

                                     tGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGAt                                                   >  2:3358673/1‑93 (MQ=255)

                                     tGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGAt                                                   >  1:1157353/1‑93 (MQ=255)

                                     tGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGAt                                                   <  2:4666943/93‑1 (MQ=255)

                                     tGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGAt                                                   <  1:320488/93‑1 (MQ=255)

                                     tGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGAt                                                   >  1:765071/1‑93 (MQ=255)

                                     tGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGAt                                                   <  1:1238253/93‑1 (MQ=255)

                                     tGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGAt                                                   <  2:2716583/93‑1 (MQ=255)

                                      gTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATa                                                  <  2:3223370/93‑1 (MQ=255)

                                      gTAGGTAAAGCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATa                                                  <  1:1712368/93‑1 (MQ=255)

                                       tAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAc                                                 >  1:1326548/1‑93 (MQ=255)

                                       tAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATaa                                                 >  1:4530361/1‑92 (MQ=255)

                                         ggTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATActc                                               >  2:4425895/1‑91 (MQ=255)

                                         ggTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgc                                               <  2:3223999/93‑1 (MQ=255)

                                           tAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcgg                                             >  1:3874612/1‑93 (MQ=255)

                                           tAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcgg                                             >  2:4600449/1‑93 (MQ=255)

                                             aaTCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggcg                                           >  1:921249/1‑93 (MQ=255)

                                               tCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggcgcg                                         >  2:3161625/1‑93 (MQ=255)

                                               tCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATAcgcggcgcg                                         <  1:2533319/93‑1 (MQ=255)

                                               tCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCCCCAGAGCAATCAAACGATATTCAAGCAGAAATTGCCAGATAcgcggcgcg                                         >  1:3746601/1‑93 (MQ=255)

                                                 aaaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGa                                       >  2:1713564/1‑93 (MQ=255)

                                                 aaaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGa                                       <  2:241386/93‑1 (MQ=255)

                                                 aaaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGa                                       >  1:3886780/1‑93 (MQ=255)

                                                  aaGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGaa                                      >  1:1183308/1‑93 (MQ=255)

                                                   aGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGaaa                                     >  1:2940122/1‑93 (MQ=255)

                                                    gCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAAt                                    >  2:4006023/1‑93 (MQ=255)

                                                     cTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAAtt                                   <  1:3064499/93‑1 (MQ=255)

                                                       gCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTaa                                 <  1:3081081/93‑1 (MQ=255)

                                                        cAACAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAAt                                >  1:434629/1‑93 (MQ=255)

                                                         aaCAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATa                               <  2:2679186/93‑1 (MQ=255)

                                                          aCAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATAc                              >  2:2299036/1‑93 (MQ=255)

                                                          aCAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATAc                              <  2:3075615/93‑1 (MQ=255)

                                                          aCAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATAc                              <  2:3132478/93‑1 (MQ=255)

                                                          aCAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATAc                              <  2:2242107/93‑1 (MQ=255)

                                                          aCAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATAc                              <  1:2744039/93‑1 (MQ=255)

                                                           cAGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACg                             >  2:3691731/1‑93 (MQ=255)

                                                            aGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGc                            <  2:134480/93‑1 (MQ=255)

                                                            aGCCGCCTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGc                            <  1:461345/93‑1 (MQ=255)

                                                             gccgccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCt                           >  1:3035351/1‑93 (MQ=255)

                                                              ccgccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTc                          >  1:3203741/1‑93 (MQ=255)

                                                              ccgccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTc                          >  2:4432628/1‑93 (MQ=255)

                                                                gccTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTccc                        >  1:674239/1‑93 (MQ=255)

                                                                  cTGGGATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccg                      >  2:2885177/1‑93 (MQ=255)

                                                                   tGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccgc                     >  2:3989633/1‑93 (MQ=255)

                                                                      gaggCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCcgccgc                  <  2:981822/90‑1 (MQ=255)

                                                                      gATGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCcgccgc                  <  2:700132/93‑1 (MQ=255)

                                                                      gATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCcgccgc                  >  1:4105579/1‑93 (MQ=255)

                                                                       aTGCTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCt                 <  2:511475/93‑1 (MQ=255)

                                                                          cTCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCtttt              <  1:4189610/93‑1 (MQ=255)

                                                                           tCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTa             >  2:2245526/1‑93 (MQ=255)

                                                                           tCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTa             <  2:3276193/93‑1 (MQ=255)

                                                                           tCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTa             >  1:1544327/1‑93 (MQ=255)

                                                                           tCAGTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTa             >  2:3574260/1‑93 (MQ=255)

                                                                             aGTCGCTGACGCGCAAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAAc           <  1:4072450/93‑1 (MQ=255)

                                                                              gTCGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAAct          <  2:4093574/93‑1 (MQ=255)

                                                                                cGCTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAActgc        <  1:1384235/93‑1 (MQ=255)

                                                                                  cTGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAActgctg      <  2:2273194/93‑1 (MQ=255)

                                                                                   tGACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGa     >  2:4589708/1‑93 (MQ=255)

                                                                                    gACGCGCCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGat    >  1:2322300/1‑92 (MQ=255)

                                                                                      cgcgCCAAAGCAATCAAACGCTCTTCACGCAGAAATTACCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCATTTAACTGCTGAAtt  >  2:3058571/1‑93 (MQ=255)

CGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATT  >  NC_000913/112149‑112327

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