BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,473,878	G→A	P116P (CCG→CCA)	b2362 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,473,878	0	G	A	85.1%	578.9 / 70.5	201	P116P (CCG→CCA)	b2362	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (15/15); new base A (75/96); total (90/111)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 205 total aligned reads displayed.

AAAAGAAACCGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGAT  >  NC_000913/2473788‑2473967

aaaaGAAACCGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGAACCAGa                                                                                         >  2:4377742/1‑93 (MQ=255)

   aGAAACCGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGc                                                                                      >  1:4099546/1‑93 (MQ=255)

    gAAACCGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCt                                                                                     >  2:5637070/1‑93 (MQ=255)

      aaCCGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTgg                                                                                   >  1:4226164/1‑93 (MQ=255)

       aCCGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGa                                                                                  <  1:1372919/93‑1 (MQ=255)

       aCCGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGa                                                                                  >  1:2463730/1‑93 (MQ=255)

         cGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGc                                                                                <  1:3625246/93‑1 (MQ=255)

         cGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGc                                                                                <  1:3327202/93‑1 (MQ=255)

         cGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGc                                                                                >  1:2858627/1‑93 (MQ=255)

          ggAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCa                                                                               >  1:6655409/1‑93 (MQ=255)

          ggAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCa                                                                               <  2:5392476/93‑1 (MQ=255)

          ggAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCa                                                                               >  1:1476259/1‑93 (MQ=255)

              aaTGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAAt                                                                           <  1:1903579/93‑1 (MQ=255)

              aaTGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAAt                                                                           >  2:3924635/1‑93 (MQ=255)

              aaTGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAAt                                                                           <  2:538354/93‑1 (MQ=255)

               aTGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATg                                                                          >  2:3189899/1‑93 (MQ=255)

                tGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGa                                                                         >  2:5173699/1‑93 (MQ=255)

                tGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGa                                                                         >  2:5592240/1‑93 (MQ=255)

                 ggCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGAc                                                                        <  1:2238836/93‑1 (MQ=255)

                 ggCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGAc                                                                        >  1:6553407/1‑93 (MQ=255)

                   cAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACca                                                                      >  1:4532147/1‑93 (MQ=255)

                     gTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACcaca                                                                    >  2:3610286/1‑93 (MQ=255)

                      tGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAACGACCACAt                                                                   <  1:6501360/93‑1 (MQ=255)

                       ggACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACAtt                                                                  >  2:1208117/1‑93 (MQ=255)

                       ggACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACAtt                                                                  <  1:1285234/93‑1 (MQ=255)

                         aCATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTa                                                                >  2:4585338/1‑93 (MQ=255)

                           aTTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTAcc                                                              >  2:2230203/1‑93 (MQ=255)

                             ttCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCaa                                                            <  1:5058615/93‑1 (MQ=255)

                              tCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAAt                                                           <  1:3882151/93‑1 (MQ=255)

                              tCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAAt                                                           >  2:1867121/1‑93 (MQ=255)

                               cAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATa                                                          <  2:5726848/93‑1 (MQ=255)

                                 ggTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATaaa                                                        <  2:1816462/93‑1 (MQ=255)

                                 ggTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATaaa                                                        >  2:797825/1‑93 (MQ=255)

                                  gTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAg                                                       <  2:5571029/93‑1 (MQ=255)

                                    tGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGaa                                                     >  2:6017548/1‑93 (MQ=255)

                                    tGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGaa                                                     <  2:4271047/93‑1 (MQ=255)

                                     ggggACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAAc                                                    >  1:1563918/1‑93 (MQ=255)

                                      gggACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACt                                                   <  2:5117502/93‑1 (MQ=255)

                                       ggACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTg                                                  <  1:5879154/93‑1 (MQ=255)

                                       ggACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTg                                                  <  2:6354575/93‑1 (MQ=255)

                                       ggACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTg                                                  <  2:411564/93‑1 (MQ=255)

                                        gACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGa                                                 <  1:785863/93‑1 (MQ=255)

                                         aCCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGAt                                                <  2:1039287/93‑1 (MQ=255)

                                          ccTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGAtt                                               <  2:463592/93‑1 (MQ=255)

                                           cTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTa                                              >  2:2282309/1‑93 (MQ=255)

                                           cTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTa                                              >  1:2727410/1‑93 (MQ=255)

                                             tGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATATAGAACTGATTaaa                                            <  2:1953639/93‑1 (MQ=255)

                                             tGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTaaa                                            <  1:3965822/93‑1 (MQ=255)

                                              gAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAg                                           <  1:4757012/93‑1 (MQ=255)

                                               aaaaCGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGa                                          <  2:4306752/93‑1 (MQ=255)

                                                aaaCGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGaa                                         >  1:2409641/1‑93 (MQ=255)

                                                aaaCGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGaa                                         <  1:5999404/93‑1 (MQ=255)

                                                   cGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATc                                      <  2:5893923/93‑1 (MQ=255)

                                                     tCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCaa                                    <  2:4040216/93‑1 (MQ=255)

                                                     tCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCaa                                    >  1:2182876/1‑93 (MQ=255)

                                                      cAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCaaa                                   <  1:2966498/93‑1 (MQ=255)

                                                      cAATATCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCaaa                                   <  1:1079555/93‑1 (MQ=255)

                                                        atatCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGa                                 <  2:14517/93‑1 (MQ=255)

                                                         tatCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGaa                                <  1:6233354/93‑1 (MQ=255)

                                                         tatCACTCACTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGaa                                <  1:3644895/93‑1 (MQ=255)

                                                           tcactcacTGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACg                              <  2:5612931/93‑1 (MQ=255)

                                                            cactcactGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGa                             <  2:6635376/93‑1 (MQ=255)

                                                             actcactGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGaa                            <  1:1427519/93‑1 (MQ=255)

                                                              ctcactGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAAt                           >  2:603709/1‑93 (MQ=255)

                                                              ctcactGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAAt                           >  2:2823630/1‑93 (MQ=255)

                                                              atcactGGGCAGTAAATGCAAAAGCACCGGAGGCTGGAGCATAATGGCCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAAt                           <  2:5072300/92‑1 (MQ=255)

                                                               tcactGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATc                          <  2:3669513/93‑1 (MQ=255)

                                                               tcactGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATc                          <  1:199203/93‑1 (MQ=255)

                                                                cactGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATca                         >  1:564542/1‑93 (MQ=255)

                                                                cactGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATca                         >  2:5303973/1‑93 (MQ=255)

                                                                cactGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATca                         <  1:4288929/93‑1 (MQ=255)

                                                                cactGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAAAGAATca                         >  1:2630289/1‑93 (MQ=255)

                                                                  ctGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATcagc                       >  1:927351/1‑93 (MQ=255)

                                                                   tGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATcagca                      >  1:3944577/1‑93 (MQ=255)

                                                                   tGGGCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATcagca                      >  1:1871871/1‑93 (MQ=255)

                                                                    gggCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATcagcag                     >  1:3502467/1‑93 (MQ=255)

                                                                    gggCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATcagcag                     >  1:4988777/1‑93 (MQ=255)

                                                                     ggCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCagcagc                    >  2:5505072/1‑93 (MQ=255)

                                                                     ggCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCagcagc                    >  2:6404608/1‑93 (MQ=255)

                                                                      gCAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGcc                   >  1:6342495/1‑93 (MQ=255)

                                                                       cAGTAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCt                  <  2:2991771/93‑1 (MQ=255)

                                                                        aGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTa                 >  2:1627407/1‑93 (MQ=255)

                                                                          tAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTaga               >  2:1727614/1‑93 (MQ=255)

                                                                          tAATGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTaga               >  2:571053/1‑93 (MQ=255)

                                                                           aaTGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTagag              <  2:3986935/93‑1 (MQ=255)

                                                                           aaTGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTagag              <  1:5480700/93‑1 (MQ=255)

                                                                            aTGCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAgg             >  2:2485072/1‑93 (MQ=255)

                                                                              gCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGtt           <  1:385241/93‑1 (MQ=255)

                                                                              gCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGtt           <  2:4826995/93‑1 (MQ=255)

                                                                              gCCAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGtt           <  2:1473410/93‑1 (MQ=255)

                                                                               ccAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTc          <  1:1532026/93‑1 (MQ=255)

                                                                                cAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCg         >  1:1820922/1‑93 (MQ=255)

                                                                                cAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCg         >  1:3597988/1‑93 (MQ=255)

                                                                                cAAAAGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCg         >  2:3455923/1‑93 (MQ=255)

                                                                                 aaaaGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCga        <  1:5706726/93‑1 (MQ=255)

                                                                                 aaaaGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCga        >  2:4684305/1‑93 (MQ=255)

                                                                                   aaGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCgaga      <  1:3032336/93‑1 (MQ=255)

                                                                                    aGCACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGAc     <  1:3607491/93‑1 (MQ=255)

                                                                                      cACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGa   <  2:4516643/93‑1 (MQ=255)

                                                                                      cACCAGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGa   >  2:2793577/1‑93 (MQ=255)

AAAAGAAACCGGAGAATGGCAGTGGACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGAT  >  NC_000913/2473788‑2473967

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