BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,246,351	T→C	S68G (AGT→GGT)	lysP ←	lysine‑specific permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,246,351	0	T	C	99.7%	928.0 / NA	297	S68G (AGT→GGT)	lysP	lysine‑specific permease
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (162/134); total (163/134)
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.99e-01

Only 100 of 299 total aligned reads displayed.

GCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGG  >  NC_000913/2246261‑2246440

   cgAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtca                                                                                      >  2:6505857/1‑93 (MQ=255)

   cgAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtca                                                                                      >  2:5568215/1‑93 (MQ=255)

      aGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtcatca                                                                                   >  1:5053787/1‑93 (MQ=255)

       gCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAg                                                                                  <  2:816003/93‑1 (MQ=255)

         cAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGa                                                                                <  2:7318509/93‑1 (MQ=255)

         cAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGa                                                                                >  1:4079590/1‑93 (MQ=255)

          aaaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGaa                                                                               <  1:2221085/93‑1 (MQ=255)

          aaaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGaa                                                                               >  2:4589203/1‑93 (MQ=255)

           aaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAg                                                                              >  2:859863/1‑93 (MQ=255)

           aaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAg                                                                              >  1:6278489/1‑93 (MQ=255)

             gCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTa                                                                            >  2:6655309/1‑93 (MQ=255)

             gCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTa                                                                            >  2:4985514/1‑93 (MQ=255)

               cttcttCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTaaa                                                                          >  2:1153605/1‑93 (MQ=255)

                ttcttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAc                                                                         <  2:1665979/93‑1 (MQ=255)

                ttcttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAc                                                                         <  1:3390248/93‑1 (MQ=255)

                 tcttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAcc                                                                        >  2:2654351/1‑93 (MQ=255)

                  cttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCa                                                                       <  2:2927544/93‑1 (MQ=255)

                  cttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCa                                                                       <  1:6663890/93‑1 (MQ=255)

                   ttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCAt                                                                      <  1:3507508/93‑1 (MQ=255)

                   ttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCAt                                                                      >  1:7377792/1‑93 (MQ=255)

                   ttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCAt                                                                      >  1:4141843/1‑93 (MQ=255)

                    tcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATc                                                                     <  1:894826/93‑1 (MQ=255)

                     cAACATAGTTCTGACCGTAAGTGGCAAACGACCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCa                                                                    <  1:7790965/93‑1 (MQ=255)

                      caCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAg                                                                   <  2:5347336/92‑1 (MQ=255)

                      aaCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAg                                                                   <  1:3790711/93‑1 (MQ=255)

                      aaCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAg                                                                   <  2:670714/93‑1 (MQ=255)

                        cATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGc                                                                 <  1:2630847/93‑1 (MQ=255)

                        cATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGc                                                                 >  1:1336689/1‑93 (MQ=255)

                         aTAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGcc                                                                <  1:5665033/93‑1 (MQ=255)

                            gTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAAt                                                             <  1:2438738/93‑1 (MQ=255)

                             ttCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATc                                                            >  2:7766739/1‑93 (MQ=255)

                                tGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGc                                                         >  2:6838898/1‑93 (MQ=255)

                                 gACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCa                                                        <  1:1108610/93‑1 (MQ=255)

                                 gACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCa                                                        <  1:7190230/93‑1 (MQ=255)

                                   ccGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCatc                                                      <  2:4794397/93‑2 (MQ=255)

                                   ccGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCata                                                      <  1:652135/93‑1 (MQ=255)

                                   ccGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCata                                                      <  1:2898987/93‑1 (MQ=255)

                                    cGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCatat                                                     <  1:3956024/93‑1 (MQ=255)

                                    cGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCatat                                                     <  2:7603864/93‑1 (MQ=255)

                                     gTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCAtata                                                    >  2:3581046/1‑93 (MQ=255)

                                     gTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCAtata                                                    <  2:7390100/93‑1 (MQ=255)

                                     gTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATAAGCAtata                                                    >  2:2083786/1‑93 (MQ=255)

                                      tAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATAc                                                   >  2:7191496/1‑93 (MQ=255)

                                      tAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATAc                                                   <  1:5404651/93‑1 (MQ=255)

                                      tAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATAc                                                   >  1:2819569/1‑93 (MQ=255)

                                       aaGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACg                                                  <  1:5507456/93‑1 (MQ=255)

                                        aGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACga                                                 >  1:3339303/1‑93 (MQ=255)

                                         gTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACgag                                                >  1:927910/1‑93 (MQ=255)

                                         gTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACgag                                                >  1:3103102/1‑93 (MQ=255)

                                           ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGagag                                              >  2:3264865/1‑93 (MQ=255)

                                           ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGagag                                              <  1:859773/93‑1 (MQ=255)

                                            gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGc                                             <  1:3536469/93‑1 (MQ=255)

                                              aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCaa                                           <  2:1440202/93‑1 (MQ=255)

                                              aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCaa                                           >  2:317518/1‑93 (MQ=255)

                                               atctaaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAAt                                          <  1:493013/89‑1 (MQ=255)

                                               aaCGACCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAAt                                          <  2:4311252/93‑1 (MQ=255)

                                                aCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATg                                         >  2:3086236/1‑93 (MQ=255)

                                                 cGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGc                                        <  2:5243175/93‑1 (MQ=255)

                                                  gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGcc                                       >  1:2967806/1‑93 (MQ=255)

                                                  gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGcc                                       >  1:4404829/1‑93 (MQ=255)

                                                  gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGcc                                       >  2:2509169/1‑93 (MQ=255)

                                                   aaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGccc                                      >  2:3973894/1‑93 (MQ=255)

                                                   aaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGccc                                      >  2:5731259/1‑93 (MQ=255)

                                                   aaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGccc                                      <  1:416106/93‑1 (MQ=255)

                                                       ggAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgcc                                  >  1:4469596/1‑93 (MQ=255)

                                                       ggAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgcc                                  <  1:1187403/93‑1 (MQ=255)

                                                       ggAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgcc                                  <  2:4698098/93‑1 (MQ=255)

                                                        gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGATAGCAATGCccggccc                                 >  2:5545891/1‑93 (MQ=255)

                                                         aaaCCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgcccg                                >  1:1409534/1‑93 (MQ=255)

                                                          aaCCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgcccgc                               >  1:1827605/1‑92 (MQ=255)

                                                          aaCCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCgg                               >  1:151974/1‑93 (MQ=255)

                                                           aCCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGa                              <  2:2897480/93‑1 (MQ=255)

                                                            ccGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAc                             <  2:2346948/93‑1 (MQ=255)

                                                            ccGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAc                             >  2:6338535/1‑93 (MQ=255)

                                                            ccGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAc                             >  2:2726613/1‑93 (MQ=255)

                                                             cGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcc                            >  2:5975852/1‑93 (MQ=255)

                                                              ggCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcct                           <  1:7732523/93‑1 (MQ=255)

                                                               gCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctg                          >  2:6202666/1‑93 (MQ=255)

                                                                 atatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcc                        >  2:2996916/1‑93 (MQ=255)

                                                                  tatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcct                       >  1:1009238/1‑93 (MQ=255)

                                                                  tatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcct                       >  2:6596042/1‑93 (MQ=255)

                                                                  tatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcct                       >  1:6143828/1‑93 (MQ=255)

                                                                  tatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcct                       >  2:2959535/1‑93 (MQ=255)

                                                                    tatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTga                     >  1:3287183/1‑93 (MQ=255)

                                                                    tatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTga                     >  2:3323864/1‑93 (MQ=255)

                                                                        caaccAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGaaa                 <  1:7825223/90‑1 (MQ=255)

                                                                         aGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAAt                <  2:1779092/93‑1 (MQ=255)

                                                                          gCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATc               <  1:1972200/93‑1 (MQ=255)

                                                                            cAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGt             <  1:7646802/93‑1 (MQ=255)

                                                                             aGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGtt            >  1:6075629/1‑93 (MQ=255)

                                                                             aGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGtt            >  2:3982752/1‑93 (MQ=255)

                                                                              gTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTg           >  1:4063431/1‑93 (MQ=255)

                                                                              gTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTg           >  1:4818325/1‑93 (MQ=255)

                                                                               ttCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTgc          >  2:4414490/1‑93 (MQ=255)

                                                                                 cACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTgcgc        >  2:6986230/1‑93 (MQ=255)

                                                                                 cACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTgcgc        <  2:7568209/93‑1 (MQ=255)

                                                                                  aCCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGcc       <  1:42095/93‑1 (MQ=255)

                                                                                     gagaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCaga    >  2:7879967/1‑93 (MQ=255)

                                                                                      agaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCagag   <  2:3858472/93‑1 (MQ=255)

                                                                                       gaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAgg  <  1:3593414/93‑1 (MQ=255)

GCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGG  >  NC_000913/2246261‑2246440

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