BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,259,418	G→A	G194G (GGC→GGT)	ychM ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,259,418	0	G	A	84.5%	475.9 / 62.7	193	G194G (GGC→GGT)	ychM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (11/19); new base A (72/91); total (83/110)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-04

Only 100 of 195 total aligned reads displayed.

GATGGTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGCCATAAATAATGC  >  NC_000913/1259328‑1259510

gATGGTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCATCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGcc                                                                                            <  2:2410403/93‑1 (MQ=255)

gATGGTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGcc                                                                                            <  1:7299934/93‑1 (MQ=255)

 aTGGTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATAccc                                                                                           >  1:5620378/1‑93 (MQ=255)

 aTGGTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATAccc                                                                                           <  2:512868/93‑1 (MQ=255)

  tGGTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCa                                                                                          >  1:7114503/1‑93 (MQ=255)

   ggTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAg                                                                                         <  2:5070833/93‑1 (MQ=255)

    gTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGCCCAGa                                                                                        <  2:3317926/93‑1 (MQ=255)

    gTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGa                                                                                        <  2:3995850/93‑1 (MQ=255)

    gTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGa                                                                                        <  2:4138279/93‑1 (MQ=255)

    gTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGa                                                                                        <  2:4471215/93‑1 (MQ=255)

    gTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGa                                                                                        >  2:2782030/1‑93 (MQ=255)

    gTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGa                                                                                        <  1:3869494/93‑1 (MQ=255)

    gTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGa                                                                                        >  2:5049186/1‑93 (MQ=255)

      aGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGAcg                                                                                      <  2:3498973/93‑1 (MQ=255)

          aCATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGc                                                                                  >  1:4982736/1‑93 (MQ=255)

             tGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCaa                                                                               >  1:6399432/1‑93 (MQ=255)

             tGTCCGCCGAGCAGGTTAA‑AATTCCCATTACCGCGCAACCAGCCAGCAAGGCCCGGAAGG‑GACCGGGGAAACCGATAACCAAGCGCGGCCcaa                                                                              >  1:3986001/1‑93 (MQ=255)

              gTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCaaa                                                                              >  2:5262500/1‑93 (MQ=255)

              gTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCaaa                                                                              >  1:5885185/1‑93 (MQ=255)

               tCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCaaaa                                                                             <  1:2084148/93‑1 (MQ=255)

               tCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCaaaa                                                                             <  1:6265005/93‑1 (MQ=255)

                 cgccgAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACAGGGTAAACGAATACCCAGACGCGGCCaaaaaa                                                                           >  1:5986717/1‑93 (MQ=255)

                  gccgAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAAc                                                                          >  1:826876/1‑93 (MQ=255)

                  gccgAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAAc                                                                          <  2:5248485/93‑1 (MQ=255)

                  gccgAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAAc                                                                          <  1:3430334/93‑1 (MQ=255)

                  gccgAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAAc                                                                          >  1:5150961/1‑93 (MQ=255)

                    cgAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACaa                                                                        >  1:5067409/1‑93 (MQ=255)

                    cgAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCGAAAAACaa                                                                        >  2:2705881/1‑93 (MQ=255)

                    cgAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACaa                                                                        >  2:7082796/1‑93 (MQ=255)

                     gAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAg                                                                       >  2:4709890/1‑93 (MQ=255)

                       gCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGaa                                                                     <  1:945962/93‑1 (MQ=255)

                       gCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGaa                                                                     <  1:4817394/93‑1 (MQ=255)

                       gCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGaa                                                                     >  2:6333332/1‑93 (MQ=255)

                        cAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAAt                                                                    >  2:1475703/1‑93 (MQ=255)

                         aGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATa                                                                   >  1:4961/1‑93 (MQ=255)

                         aGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATa                                                                   >  2:4254359/1‑93 (MQ=255)

                          ggTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATAc                                                                  <  2:6016205/93‑1 (MQ=255)

                          ggTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATAc                                                                  <  2:5705514/93‑1 (MQ=255)

                           gTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATAcc                                                                 >  2:1412180/1‑93 (MQ=255)

                           gTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATAcc                                                                 >  2:7324808/1‑93 (MQ=255)

                           gTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATAcc                                                                 <  1:7430084/93‑1 (MQ=255)

                           gTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATAcc                                                                 >  2:81903/1‑93 (MQ=255)

                            ttAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCt                                                                <  1:1757563/93‑1 (MQ=255)

                             tAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTa                                                               <  1:6718090/93‑1 (MQ=255)

                             tAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTa                                                               <  1:4021978/93‑1 (MQ=255)

                              aaCAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAg                                                              <  1:4217622/93‑1 (MQ=255)

                                  aTCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTc                                                          >  1:5497965/1‑93 (MQ=255)

                                  aTCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTc                                                          >  1:4205054/1‑93 (MQ=255)

                                     cccATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTcaca                                                       >  1:1045196/1‑93 (MQ=255)

                                      ccATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACaa                                                      <  1:1739894/93‑1 (MQ=255)

                                       cATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAAt                                                     >  2:4225718/1‑93 (MQ=255)

                                       cATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAAt                                                     >  2:6474547/1‑93 (MQ=255)

                                        aTCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATg                                                    >  1:1144817/1‑93 (MQ=255)

                                         tCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGc                                                   >  1:5298726/1‑93 (MQ=255)

                                           aCCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCa                                                 >  1:3285212/1‑93 (MQ=255)

                                             cgcgCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAAt                                               >  2:2331676/1‑93 (MQ=255)

                                               cgcAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATgg                                             <  1:7409861/93‑1 (MQ=255)

                                                gcAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                            <  1:6204480/93‑1 (MQ=255)

                                                 cAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgca                                           >  1:459843/1‑93 (MQ=255)

                                                 cAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgca                                           <  2:582132/93‑1 (MQ=255)

                                                  aaCCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgcag                                          <  2:3665/93‑1 (MQ=255)

                                                  aaCCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgcag                                          <  1:258764/93‑1 (MQ=255)

                                                   aCCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgcagc                                         >  1:1305329/1‑93 (MQ=255)

                                                     cagccagcAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCAt                                       >  1:3090859/1‑93 (MQ=255)

                                                      agccagcAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATc                                      >  2:3855741/1‑93 (MQ=255)

                                                      agccagcAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATc                                      <  2:4934212/93‑1 (MQ=255)

                                                       gccagcAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCa                                     <  2:2619227/93‑1 (MQ=255)

                                                         cagcAAGG‑CCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCAcc                                   <  2:5579294/93‑1 (MQ=255)

                                                         cagcAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCAcc                                   <  1:6107510/93‑1 (MQ=255)

                                                         cagcAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCAcc                                   <  2:1056669/93‑1 (MQ=255)

                                                           gcAAGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCca                                 <  2:6532365/93‑1 (MQ=255)

                                                             aaGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCcaca                               <  2:745911/93‑1 (MQ=255)

                                                              aGG‑CCGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAAT‑GCAGCATCACCCACAtt                             >  1:1617237/1‑93 (MQ=255)

                                                                  ccGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTaa                           <  2:5497050/93‑1 (MQ=255)

                                                                  ccGGAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTaa                           <  1:7194763/93‑1 (MQ=255)

                                                                    ggAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATg                         >  1:2619254/1‑93 (MQ=255)

                                                                     gAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATgg                        <  2:7374130/93‑1 (MQ=255)

                                                                     gAAGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATgg                        <  2:1377029/93‑1 (MQ=255)

                                                                      aaGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGt                       <  1:3713175/93‑1 (MQ=255)

                                                                      aaGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGt                       <  2:827113/93‑1 (MQ=255)

                                                                      aaGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGt                       <  1:2184834/93‑1 (MQ=255)

                                                                      aaGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGt                       <  1:6515199/93‑1 (MQ=255)

                                                                       aGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTc                      <  2:4184046/93‑1 (MQ=255)

                                                                       aGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTc                      <  1:4717063/93‑1 (MQ=255)

                                                                       aGGTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTc                      <  2:2002953/93‑1 (MQ=255)

                                                                        ggTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCg                     >  2:1125644/1‑93 (MQ=255)

                                                                         gTGACCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCgg                    >  1:1986006/1‑93 (MQ=255)

                                                                            aCCGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAg                 >  1:3546316/1‑93 (MQ=255)

                                                                             ccGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAgc                <  1:5238994/93‑1 (MQ=255)

                                                                             ccGGGTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAgc                <  2:1932566/93‑1 (MQ=255)

                                                                               gggTAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCAAATTAATGGTCGGCAgcgc              <  2:3025411/93‑1 (MQ=255)

                                                                                ggTAAACGGATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGcc             >  2:4857180/1‑93 (MQ=255)

                                                                                  tAAACGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGCCat           >  1:299573/1‑93 (MQ=255)

                                                                                    aaCGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGCCataa         <  1:6782035/93‑1 (MQ=255)

                                                                                     aCGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGCCataaa        >  2:2819568/1‑93 (MQ=255)

                                                                                      cGAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGCCataaat       <  2:7249605/93‑1 (MQ=255)

                                                                                       gAATACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGCCataaata      <  1:2529671/93‑1 (MQ=255)

                                                                                        aaTACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGCCataaataa     <  1:1750371/93‑1 (MQ=255)

                                                                                        aaTACCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGCCataaataa     <  1:6374893/93‑1 (MQ=255)

                                                                                           gCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGCCATAAATAATgc  <  2:1765082/93‑1 (MQ=255)

GATGGTAGCAACATGTCCGCCGAGCAGGTTAACAATCCCCATCACCGCGCAACCAGCCAGCAAGG‑CCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGCCATAAATAATGC  >  NC_000913/1259328‑1259510

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