BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,893,993	A→G	A82A (GCT→GCC)	yieG ←	putative membrane/transport protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,893,993	0	A	G	100.0%	787.3 / NA	244	A82A (GCT→GCC)	yieG	putative membrane/transport protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (130/114); total (130/114)

Only 100 of 249 total aligned reads displayed.

GACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAA  >  NC_000913/3893903‑3894081

gACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGc                                                                                        <  2:949032/93‑1 (MQ=255)

gACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGc                                                                                        <  1:3576839/93‑1 (MQ=255)

 aCCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCa                                                                                       >  2:3928244/1‑93 (MQ=255)

  ccTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCaa                                                                                      <  2:1149955/93‑1 (MQ=255)

     gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACa                                                                                   <  2:4440179/93‑1 (MQ=255)

     gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACa                                                                                   <  1:3910859/93‑1 (MQ=255)

     gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACa                                                                                   <  1:1635103/93‑1 (MQ=255)

        aCGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTc                                                                                <  2:3484063/93‑1 (MQ=255)

         cGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTcc                                                                               <  2:2028912/93‑1 (MQ=255)

           gCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCca                                                                             <  2:18122/93‑1 (MQ=255)

           gCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCca                                                                             <  1:455494/93‑1 (MQ=255)

           gCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCca                                                                             <  2:3759879/93‑1 (MQ=255)

           gCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGGTGGCAAACAGTCCca                                                                             >  2:3035881/1‑93 (MQ=255)

            cAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcat                                                                            <  2:428001/93‑1 (MQ=255)

             aaGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcatc                                                                           <  2:1676119/93‑1 (MQ=255)

                 ccATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataa                                                                       >  1:719496/1‑93 (MQ=255)

                 ccATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataa                                                                       <  1:2526879/93‑1 (MQ=255)

                  cATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataat                                                                      >  1:931472/1‑93 (MQ=255)

                  cATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataat                                                                      >  1:4199128/1‑93 (MQ=255)

                  cATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataat                                                                      >  1:307125/1‑93 (MQ=255)

                  cATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataat                                                                      >  1:1851565/1‑93 (MQ=255)

                    tCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATAc                                                                    <  2:2758020/93‑1 (MQ=255)

                    tCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATAc                                                                    <  2:1246572/93‑1 (MQ=255)

                    tCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATAc                                                                    >  1:3398823/1‑93 (MQ=255)

                     cGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACt                                                                   <  1:2420770/93‑1 (MQ=255)

                       ccTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGc                                                                 >  1:140892/1‑93 (MQ=255)

                         tGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCg                                                               >  2:1428941/1‑93 (MQ=255)

                         tGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCg                                                               >  1:443193/1‑93 (MQ=255)

                           tACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTTCCGaa                                                             <  2:2686280/93‑1 (MQ=255)

                           tACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGaa                                                             <  1:409806/93‑1 (MQ=255)

                           tACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGaa                                                             <  1:324929/93‑1 (MQ=255)

                           tACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGaa                                                             >  1:3262197/1‑93 (MQ=255)

                            aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAAt                                                            <  1:836710/93‑1 (MQ=255)

                            aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAAt                                                            <  1:2097780/93‑1 (MQ=255)

                            aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAAt                                                            <  2:63045/93‑1 (MQ=255)

                            aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAAt                                                            <  2:1313922/93‑1 (MQ=255)

                              aaCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGc                                                          >  2:1979466/1‑93 (MQ=255)

                              aaCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGc                                                          <  1:1032077/93‑1 (MQ=255)

                               aCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCa                                                         >  2:4607963/1‑93 (MQ=255)

                               aCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCa                                                         <  2:3095490/93‑1 (MQ=255)

                               aCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCa                                                         <  2:650290/93‑1 (MQ=255)

                                cGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAg                                                        >  2:2603773/1‑93 (MQ=255)

                                cGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGGGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAg                                                        >  1:201298/1‑93 (MQ=255)

                                 gACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGc                                                       >  2:4408512/1‑93 (MQ=255)

                                 gACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGc                                                       <  2:7457/93‑1 (MQ=255)

                                 gACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGc                                                       >  2:1481358/1‑93 (MQ=255)

                                  aCAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCg                                                      >  2:4662363/1‑93 (MQ=255)

                                  aCAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCg                                                      <  2:1083086/93‑1 (MQ=255)

                                   cAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGa                                                     >  1:3961454/1‑93 (MQ=255)

                                   cAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGa                                                     >  2:3479912/1‑93 (MQ=255)

                                   cAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGa                                                     >  1:3946585/1‑93 (MQ=255)

                                     aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                   <  1:1887289/93‑1 (MQ=255)

                                     aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                   <  2:1936705/93‑1 (MQ=255)

                                     aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                   >  2:1263115/1‑93 (MQ=255)

                                     aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                   >  1:2480554/1‑93 (MQ=255)

                                      aaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCa                                                  >  1:1481311/1‑93 (MQ=255)

                                      aaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCa                                                  <  2:2379903/93‑1 (MQ=255)

                                      aaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCa                                                  >  1:1983429/1‑93 (MQ=255)

                                      aaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCa                                                  <  1:12063/93‑1 (MQ=255)

                                       aGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAg                                                 <  1:2972798/93‑1 (MQ=255)

                                       aGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAg                                                 >  2:2172633/1‑93 (MQ=255)

                                          gaagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACa                                              >  2:3551632/1‑93 (MQ=255)

                                            agaaCTCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAgg                                            <  2:1916060/93‑1 (MQ=255)

                                            agaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAgg                                            >  2:409326/1‑93 (MQ=255)

                                                 gCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTg                                       <  2:3189718/93‑1 (MQ=255)

                                                 gCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTg                                       >  2:4473787/1‑93 (MQ=255)

                                                 gCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTg                                       >  1:2379494/1‑93 (MQ=255)

                                                  cATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGa                                      >  1:470347/1‑93 (MQ=255)

                                                  cATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGa                                      <  1:3739502/93‑1 (MQ=255)

                                                    ttCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACg                                    >  2:1559313/1‑93 (MQ=255)

                                                    ttCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACg                                    >  1:2747948/1‑93 (MQ=255)

                                                    ttCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGCGACg                                    >  1:3124735/1‑93 (MQ=255)

                                                     tCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGa                                   >  1:2810937/1‑93 (MQ=255)

                                                         gCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGAc                               <  1:795988/93‑1 (MQ=255)

                                                            cATAGCGGGGGCCAGAGCAACTGGCAGGTTGGCAAACAGTCCCATCATACTACTGCCGAATGCAGCGATCAGACAGGTAGAGACGAAGACgtg                            >  2:1366577/1‑91 (MQ=255)

                                                            cATAGCGGGAGCCAGTGCAACTGGCAGGGTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgtg                            >  2:3954612/1‑91 (MQ=255)

                                                             aTAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGgcg                           >  1:1787589/1‑93 (MQ=255)

                                                              tAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcgc                          >  2:1625845/1‑93 (MQ=255)

                                                              tAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcgc                          <  1:33141/93‑1 (MQ=255)

                                                               aGCGGGGGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATTCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCt                         <  1:1057995/93‑1 (MQ=255)

                                                                 cGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTgg                       >  1:1082110/1‑93 (MQ=255)

                                                                   ggTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTa                     <  2:2252650/93‑1 (MQ=255)

                                                                   ggTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTa                     >  1:3653724/1‑93 (MQ=255)

                                                                     tGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATc                   >  1:2269647/1‑93 (MQ=255)

                                                                      gCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATcc                  >  1:4101220/1‑93 (MQ=255)

                                                                       ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                 <  1:2224504/93‑1 (MQ=255)

                                                                         aGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATg               <  2:32394/93‑1 (MQ=255)

                                                                         aGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATg               <  2:337813/93‑1 (MQ=255)

                                                                          gTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCGATGc              >  2:873194/1‑93 (MQ=255)

                                                                          gTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGc              >  1:1747782/1‑93 (MQ=255)

                                                                            gCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCa            <  2:1183943/93‑1 (MQ=255)

                                                                              aaCTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGc          >  1:4511986/1‑93 (MQ=255)

                                                                              aaCTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGc          <  1:1360138/93‑1 (MQ=255)

                                                                                cTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCaa        <  1:1196085/93‑1 (MQ=255)

                                                                                cTGGCAGGTTGGAAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCaa        <  2:3420972/93‑1 (MQ=255)

                                                                                 tGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAAc       >  1:614178/1‑93 (MQ=255)

                                                                                  ggCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACg      >  2:4104911/1‑93 (MQ=255)

                                                                                  ggCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACg      >  1:181911/1‑93 (MQ=255)

                                                                                     aGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCa   >  2:547950/1‑93 (MQ=255)

                                                                                      ggTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCaa  <  1:2612959/93‑1 (MQ=255)

GACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAA  >  NC_000913/3893903‑3894081

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: `‑`