BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,623,346	T→C	L71P (CTG→CCG)	ppx →	exopolyphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,623,346	0	T	C	93.5%	481.1 / 16.4	155	L71P (CTG→CCG)	ppx	exopolyphosphatase
Reads supporting (aligned to +/- strand): ref base T (4/6); new base C (66/79); total (70/85)
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.45e-01

Only 100 of 158 total aligned reads displayed.

CTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGAC  >  NC_000913/2623255‑2623434

 tGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGt                                                                                        >  1:4989618/1‑93 (MQ=255)

   aaaCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGttt                                                                                      <  2:4563682/93‑1 (MQ=255)

   aaaCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGttt                                                                                      <  2:2503719/93‑1 (MQ=255)

   aaaCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGttt                                                                                      <  2:157361/93‑1 (MQ=255)

     aCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGc                                                                                    <  2:213522/93‑1 (MQ=255)

       aGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCg                                                                                  >  2:4898052/1‑93 (MQ=255)

       aGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCg                                                                                  <  2:4125386/93‑1 (MQ=255)

       aGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCg                                                                                  >  1:1601344/1‑93 (MQ=255)

        gCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGa                                                                                 >  2:1233947/1‑93 (MQ=255)

         cGGGTGCATCTGGCGGACGGCGTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGaa                                                                                <  2:3150779/93‑1 (MQ=255)

         cGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGaa                                                                                <  2:5984623/93‑1 (MQ=255)

          gggTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAAc                                                                               <  1:3066859/93‑1 (MQ=255)

          gggTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAAc                                                                               <  1:5408849/93‑1 (MQ=255)

           ggTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACg                                                                              <  2:3423214/93‑1 (MQ=255)

           ggTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACg                                                                              >  2:710397/1‑93 (MQ=255)

                 tCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACa                                                                        >  2:1476579/1‑93 (MQ=255)

                 tCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACa                                                                        >  1:1085169/1‑93 (MQ=255)

                  cTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACaa                                                                       <  2:1143608/93‑1 (MQ=255)

                  cTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACaa                                                                       <  1:2474175/93‑1 (MQ=255)

                   tGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAg                                                                      >  2:5374498/1‑93 (MQ=255)

                   tGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAg                                                                      <  2:1594706/93‑1 (MQ=255)

                   tGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAg                                                                      >  2:1796637/1‑93 (MQ=255)

                      cGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGt                                                                   <  1:2836607/93‑1 (MQ=255)

                       ggACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGCACGGCTACAAgggat                                                                  >  2:2090543/1‑91 (MQ=255)

                       ggACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGtt                                                                  >  2:5656546/1‑93 (MQ=255)

                        gACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGttt                                                                 >  2:2286756/1‑93 (MQ=255)

                         aCGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGtttt                                                                <  1:5352139/93‑1 (MQ=255)

                          cGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTtc                                                               <  2:4969072/93‑1 (MQ=255)

                          cGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTtc                                                               <  1:1152000/93‑1 (MQ=255)

                           ggCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTtct                                                              <  1:5813914/93‑1 (MQ=255)

                              cTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCt                                                           >  1:279615/1‑93 (MQ=255)

                              cTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCt                                                           >  1:1090427/1‑93 (MQ=255)

                                 gggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGcc                                                        >  1:3688500/1‑93 (MQ=255)

                                  ggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCa                                                       >  1:3899071/1‑93 (MQ=255)

                                  ggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCa                                                       <  2:4802877/93‑1 (MQ=255)

                                    ccAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGc                                                     >  2:4579774/1‑93 (MQ=255)

                                    ccAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGc                                                     >  1:1419442/1‑93 (MQ=255)

                                     cAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCg                                                    <  1:1654971/93‑1 (MQ=255)

                                       gATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTc                                                  >  1:5365957/1‑93 (MQ=255)

                                       gATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTc                                                  <  1:16621/93‑1 (MQ=255)

                                       gATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTc                                                  >  1:225748/1‑93 (MQ=255)

                                         taataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTCCAAGGGTTTTCTCCTGCCAGCGTCTg                                                <  1:4018353/93‑1 (MQ=255)

                                         taataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTg                                                <  1:1849567/93‑1 (MQ=255)

                                          aataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGt                                               <  1:3749371/93‑1 (MQ=255)

                                           atatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGta                                              >  2:1347795/1‑93 (MQ=255)

                                            tatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtat                                             >  2:1692634/1‑93 (MQ=255)

                                             atGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtata                                            <  1:921331/93‑1 (MQ=255)

                                              tGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAg                                           >  2:2724783/1‑93 (MQ=255)

                                              tGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAg                                           >  2:4781626/1‑93 (MQ=255)

                                                ttGAGTGAAGAGGCAATGACGCGCGGTTTAAGCTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGtt                                         >  2:4866730/1‑93 (MQ=255)

                                                ttGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGtt                                         >  1:1445821/1‑93 (MQ=255)

                                                  gAGTGACGAGGCCATGACGCGCGGTTTAAACTGTCTGTCGCTGATTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTgg                                       <  1:5492676/93‑1 (MQ=255)

                                                  gAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTgg                                       >  1:2093217/1‑93 (MQ=255)

                                                  gAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTgg                                       <  1:5047816/93‑1 (MQ=255)

                                                   aGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGt                                      <  1:4780609/93‑1 (MQ=255)

                                                       aaGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTAccc                                  <  1:3960477/93‑1 (MQ=255)

                                                        agagGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCa                                 <  1:2244133/93‑1 (MQ=255)

                                                        agagGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCa                                 <  2:2225023/93‑1 (MQ=255)

                                                         gagGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCAt                                >  1:4632840/1‑93 (MQ=255)

                                                         gagGCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCAt                                >  2:5337396/1‑93 (MQ=255)

                                                           ggCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATAc                              <  1:5613593/93‑1 (MQ=255)

                                                           ggCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATAc                              >  2:227951/1‑93 (MQ=255)

                                                            gCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACg                             <  2:202640/93‑1 (MQ=255)

                                                            gCAATGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACg                             >  2:593907/1‑93 (MQ=255)

                                                               aTGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTg                          >  2:5115155/1‑93 (MQ=255)

                                                               aTGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTGTAGTTGGTACCCATACGCTg                          <  1:5934225/93‑1 (MQ=255)

                                                               aTGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATCCGCTg                          >  2:1688917/1‑93 (MQ=255)

                                                               aTGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTg                          >  2:5660299/1‑93 (MQ=255)

                                                                tGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGc                         >  2:2506341/1‑93 (MQ=255)

                                                                tGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGc                         <  2:1157293/93‑1 (MQ=255)

                                                                tGACGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGc                         >  2:3894677/1‑93 (MQ=255)

                                                                tGACGCACGGTTTAACCTGTCGGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGc                         <  1:48910/93‑1 (MQ=255)

                                                                  aCGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGt                       <  1:1418586/93‑1 (MQ=255)

                                                                  aCGCGCGGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGt                       <  2:4229799/93‑1 (MQ=255)

                                                                   cgcgcgGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTc                      <  1:869473/93‑1 (MQ=255)

                                                                   cgcgcgGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACGAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTc                      <  1:1104969/93‑1 (MQ=255)

                                                                   cgcgcgGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTc                      <  2:2473419/93‑1 (MQ=255)

                                                                   cgcgcgGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTc                      <  2:3225060/93‑1 (MQ=255)

                                                                   cgcgcgGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTc                      <  2:4990552/93‑1 (MQ=255)

                                                                   cgcgcgGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTc                      <  1:1758565/93‑1 (MQ=255)

                                                                    gcgcgGTTTAAACTGTCTGTCGCCGTTTGCGGAACGGCTACAAGGGTGTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCa                     <  1:3636426/93‑1 (MQ=255)

                                                                    gcgcgGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCa                     <  2:2674713/93‑1 (MQ=255)

                                                                    ccggggTTTAAACTGTCTGTCGCCGTGTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATCCGCTGCGTCa                     <  1:2043179/89‑1 (MQ=255)

                                                                      ccgGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAgg                   <  2:1766186/92‑1 (MQ=255)

                                                                       cgGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgc                  <  2:4076967/93‑1 (MQ=255)

                                                                       cgGTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgc                  <  1:331518/93‑1 (MQ=255)

                                                                        ggTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgcg                 >  2:2241560/1‑93 (MQ=255)

                                                                         gTTTAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgcgc                >  2:3906831/1‑93 (MQ=255)

                                                                         gTTTAAACTGTCTGTCGCCGGTTGCCGAACGGCCACAAGGGCTTTCTCCTGCCAGCGTCTGGATAGTTGGGACCCATACGCCGCGTCAGgcgc                >  1:4772835/1‑93 (MQ=255)

                                                                          tttAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCt               >  1:6022988/1‑93 (MQ=255)

                                                                          tttAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCt               >  1:2562337/1‑93 (MQ=255)

                                                                           ttAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTg              >  2:1022845/1‑93 (MQ=255)

                                                                            tAAACTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGa             >  1:175900/1‑93 (MQ=255)

                                                                             aaaaTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGaa            >  1:150679/1‑93 (MQ=255)

                                                                             aaaCTGTCTGTCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGaa            <  2:4242553/93‑1 (MQ=255)

                                                                                ctgtctgtCGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGc         <  1:4374061/93‑1 (MQ=255)

                                                                                  gtctgtcGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCa       <  2:5319749/93‑1 (MQ=255)

                                                                                   tctgtcGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCAc      <  1:769728/93‑1 (MQ=255)

                                                                                      gtcGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGa   >  1:3546578/1‑93 (MQ=255)

                                                                                       tcGCCGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGAc  <  2:513482/93‑1 (MQ=255)

CTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGAC  >  NC_000913/2623255‑2623434

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