BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,443,236	A→G	P516P (CCA→CCG)	ytfN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,443,236	0	A	G	99.6%	702.3 / NA	225	P516P (CCA→CCG)	ytfN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (121/103); minor base T (1/0); total (122/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 234 total aligned reads displayed.

GGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGC  >  NC_000913/4443146‑4443325

ggAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCggg                                                                                         >  1:594820/1‑93 (MQ=255)

 gAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCgggt                                                                                        <  1:167539/93‑2 (MQ=255)

 gAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCgggg                                                                                        <  1:3729609/93‑1 (MQ=255)

   gTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCt                                                                                      >  1:1250170/1‑93 (MQ=255)

   gTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCt                                                                                      >  1:4048692/1‑93 (MQ=255)

      gAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCa                                                                                   >  2:4084085/1‑93 (MQ=255)

        aGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATc                                                                                 <  1:891972/93‑1 (MQ=255)

        aGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATc                                                                                 <  1:1953508/93‑1 (MQ=255)

             aCCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGaa                                                                            >  1:3126703/1‑93 (MQ=255)

             aCCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGaa                                                                            >  2:914815/1‑93 (MQ=255)

              ccGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAAc                                                                           >  2:4440798/1‑93 (MQ=255)

              ccGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAAc                                                                           >  2:4280186/1‑93 (MQ=255)

                ggTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTc                                                                         >  1:4627224/1‑93 (MQ=255)

                 gTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATTATCCCGGGGCTTCATCTGGAACTCg                                                                        >  1:4335018/1‑93 (MQ=255)

                 gTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCg                                                                        <  1:115606/93‑1 (MQ=255)

                  tAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCgg                                                                       <  1:1057390/93‑1 (MQ=255)

                  tAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCgg                                                                       <  2:2600039/93‑1 (MQ=255)

                  tAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCgg                                                                       >  1:160496/1‑93 (MQ=255)

                   aaCGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCggg                                                                      <  2:3858439/93‑1 (MQ=255)

                    aagTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGc                                                                     >  2:2947441/3‑93 (MQ=255)

                    aCGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGc                                                                     <  2:101192/93‑1 (MQ=255)

                      gTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCa                                                                   >  2:3558052/1‑93 (MQ=255)

                       ttAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCaa                                                                  >  1:3069714/1‑93 (MQ=255)

                         aaaCAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAAc                                                                >  2:4376618/1‑93 (MQ=255)

                          aaCAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACa                                                               <  1:2390618/93‑1 (MQ=255)

                          aaCAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACa                                                               <  1:3099397/93‑1 (MQ=255)

                           aCAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAg                                                              >  2:1259506/1‑93 (MQ=255)

                           aCAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAg                                                              >  1:154996/1‑93 (MQ=255)

                            cAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGt                                                             <  1:2598139/93‑1 (MQ=255)

                            cAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGt                                                             <  1:1813314/93‑1 (MQ=255)

                            cAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGt                                                             <  2:329828/93‑1 (MQ=255)

                              gAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGc                                                           <  2:3277920/93‑1 (MQ=255)

                              gAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGc                                                           <  2:301225/93‑1 (MQ=255)

                              gAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGc                                                           >  2:3441655/1‑93 (MQ=255)

                               aaCAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGcc                                                          >  1:3988322/1‑93 (MQ=255)

                                   aaGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAg                                                      >  2:2929753/1‑93 (MQ=255)

                                   aaGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAg                                                      <  1:1790008/93‑1 (MQ=255)

                                   aaGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAg                                                      >  1:3795525/1‑93 (MQ=255)

                                    aGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGt                                                     >  1:647743/1‑93 (MQ=255)

                                    aGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGt                                                     <  1:1027418/93‑1 (MQ=255)

                                     gTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTg                                                    <  2:1605348/93‑1 (MQ=255)

                                       gAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGTTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGaa                                                  >  1:2591365/1‑93 (MQ=255)

                                       gAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGaa                                                  <  2:2046245/93‑1 (MQ=255)

                                       gAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGaa                                                  >  2:506879/1‑93 (MQ=255)

                                        aaCGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGaaa                                                 >  1:3932007/1‑93 (MQ=255)

                                         aCGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAg                                                >  2:1910772/1‑93 (MQ=255)

                                         aCGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAg                                                <  2:2670413/93‑1 (MQ=255)

                                          cGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAgg                                               <  1:2678884/93‑1 (MQ=255)

                                           gTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGc                                              >  2:132293/1‑93 (MQ=255)

                                           gTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGc                                              <  2:4562782/93‑1 (MQ=255)

                                            ttGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCg                                             >  1:1288873/1‑93 (MQ=255)

                                             tGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATAATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGa                                            <  1:1708404/93‑1 (MQ=255)

                                              gACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAg                                           >  2:2489786/1‑93 (MQ=255)

                                               aCGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGGGAGc                                          >  2:1542400/1‑93 (MQ=255)

                                               aCGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGc                                          <  1:3510498/93‑1 (MQ=255)

                                                 ggCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTg                                        <  2:4508293/93‑1 (MQ=255)

                                                 ggCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTg                                        <  1:3679825/93‑1 (MQ=255)

                                                   cACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTggg                                      <  2:1828279/93‑1 (MQ=255)

                                                    aCGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACCGTGCCGAAGTGAAAGGCGAGCTgggg                                     >  2:3607547/1‑93 (MQ=255)

                                                     cGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTggggg                                    >  1:1398356/1‑93 (MQ=255)

                                                      gCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTggggga                                   >  2:1743704/1‑92 (MQ=255)

                                                       cTGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTggcggga                                  >  1:2908509/1‑91 (MQ=255)

                                                        tGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGCAGTGAAAGGCGAGCTGGGGGTaa                                 >  1:2791321/1‑93 (MQ=255)

                                                        tGAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTaa                                 >  1:1071560/1‑93 (MQ=255)

                                                         gAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTaaa                                >  2:3954953/1‑93 (MQ=255)

                                                         gAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTaaa                                <  1:3213229/93‑1 (MQ=255)

                                                         gAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTaaa                                >  2:1914664/1‑93 (MQ=255)

                                                         gAAAGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTaaa                                <  1:2108297/93‑1 (MQ=255)

                                                          aaaGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTGATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTaaaa                               <  1:4273749/93‑1 (MQ=255)

                                                          aaaGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTaaaa                               <  2:4501106/93‑1 (MQ=255)

                                                            aGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGa                             >  1:505414/1‑93 (MQ=255)

                                                            aGGCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGa                             >  2:2977865/1‑93 (MQ=255)

                                                             ggCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGAt                            >  1:3417331/1‑93 (MQ=255)

                                                             ggCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGAt                            >  2:3591268/1‑93 (MQ=255)

                                                             ggCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGAt                            <  2:3396962/93‑1 (MQ=255)

                                                             agCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGAt                            <  1:1916724/92‑1 (MQ=255)

                                                              gCAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGAtc                           <  2:543081/93‑1 (MQ=255)

                                                               cAACAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGAtct                          <  2:3721468/93‑1 (MQ=255)

                                                                aaCAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGAtctc                         >  1:754533/1‑93 (MQ=255)

                                                                aaCAGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGAtctc                         <  2:4115104/93‑1 (MQ=255)

                                                                   aGTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAAt                      >  1:363583/1‑93 (MQ=255)

                                                                    gTTATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTACAAGATCTCAATc                     <  2:3757580/93‑1 (MQ=255)

                                                                     ttATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCt                    >  1:1986860/1‑93 (MQ=255)

                                                                     ttATATGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGAGAGCTGGGGGTAAAAGATCTCAATCt                    >  1:819939/1‑93 (MQ=255)

                                                                       atatGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTg                  <  2:2761205/93‑1 (MQ=255)

                                                                        tatGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGa                 <  2:4057071/93‑1 (MQ=255)

                                                                        tatGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGCGAAAAAGATCTCAATCTTGa                 <  2:4310320/93‑1 (MQ=255)

                                                                          tGCAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATg               >  2:2127012/1‑93 (MQ=255)

                                                                          tGCAGTGGATGATCCCGGGGCTTCATCCGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATg               <  1:3269697/93‑1 (MQ=255)

                                                                            cAGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGcc             >  2:1972062/1‑93 (MQ=255)

                                                                             aGTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGcca            <  1:3023357/93‑1 (MQ=255)

                                                                              gTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGccac           >  2:2566163/1‑93 (MQ=255)

                                                                              gTGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGccac           >  1:3428046/1‑93 (MQ=255)

                                                                               tGGATGATCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGccacc          >  1:1331560/1‑93 (MQ=255)

                                                                                  atgatCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATc       >  2:600300/1‑93 (MQ=255)

                                                                                    gatCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCaa     <  2:3218673/93‑1 (MQ=255)

                                                                                     atCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAAc    <  2:3689454/93‑1 (MQ=255)

                                                                                     atCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAAc    <  2:2682269/93‑1 (MQ=255)

                                                                                     atCCCGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAAc    <  2:935344/93‑1 (MQ=255)

                                                                                       cccGGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAAcgc  <  2:1401482/93‑1 (MQ=255)

GGAGTTGAAGCTGACCGGTAACGTTAAACAGAACAAAGTGAACGTTGACGGCACGCTGAAAGGCAACAGTTATATGCAGTGGATGATCCCAGGGCTTCATCTGGAACTCGGGCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGC  >  NC_000913/4443146‑4443325

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