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%	916.5 / NA	280	A82A (GCT→GCC)	yieG	putative membrane/transport protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (141/139); total (141/139)

Only 100 of 287 total aligned reads displayed.

GACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAG  >  NC_000913/3893903‑3894082

gACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGc                                                                                         >  1:4802107/1‑93 (MQ=255)

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

    tGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAAc                                                                                     <  1:4342827/93‑1 (MQ=255)

     gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACa                                                                                    >  1:1246410/1‑93 (MQ=255)

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

      ccACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAg                                                                                   >  1:3072728/1‑93 (MQ=255)

       cACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGt                                                                                  >  2:1431621/1‑93 (MQ=255)

         cGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTcc                                                                                >  2:1033671/1‑93 (MQ=255)

          ggCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTccc                                                                               >  1:2029353/1‑93 (MQ=255)

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

            cAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcat                                                                             >  1:4149828/1‑93 (MQ=255)

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

              aGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcatca                                                                           <  1:2514429/93‑1 (MQ=255)

               gCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcatcat                                                                          <  2:4092658/93‑1 (MQ=255)

               gCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcatcat                                                                          >  1:2605028/1‑93 (MQ=255)

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

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

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

                   aTCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCATGTTGGCAAACAGTCCCATCataata                                                                      <  1:3480674/93‑1 (MQ=255)

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

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

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

                      gCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTg                                                                   >  1:5788773/1‑93 (MQ=255)

                        cTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGcc                                                                 >  2:4479965/1‑93 (MQ=255)

                        cTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGcc                                                                 >  1:5041236/1‑93 (MQ=255)

                        cTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGcc                                                                 >  1:4652440/1‑93 (MQ=255)

                          gTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGa                                                               <  2:1195576/93‑1 (MQ=255)

                          gTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGa                                                               >  2:5612468/1‑93 (MQ=255)

                          gTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGa                                                               <  2:5365216/93‑1 (MQ=255)

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

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

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

                            aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAAt                                                             >  1:3943180/1‑93 (MQ=255)

                             cAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATg                                                            >  2:6275813/1‑93 (MQ=255)

                             cAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATg                                                            <  1:752066/93‑1 (MQ=255)

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

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

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

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

                                   cAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGa                                                      <  1:6212513/93‑1 (MQ=255)

                                    aaaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGAt                                                     >  2:3199514/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                                          gaagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACa                                               <  1:1065607/93‑1 (MQ=255)

                                           aagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAg                                              <  1:676414/93‑1 (MQ=255)

                                           aagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAg                                              <  2:683338/93‑1 (MQ=255)

                                           aagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAg                                              <  2:3987928/93‑1 (MQ=255)

                                             gaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGt                                            <  1:2101651/93‑1 (MQ=255)

                                              aaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTa                                           <  1:5243802/93‑1 (MQ=255)

                                               aCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAg                                          <  1:2986814/93‑1 (MQ=255)

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

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

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

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

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

                                                   aTTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGAc                                      <  1:2745043/93‑1 (MQ=255)

                                                   aTTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGAc                                      >  1:3400216/1‑93 (MQ=255)

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

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

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

                                                    ttCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATACAGCGATCAGACAGGTAGTGACg                                     >  2:3091710/1‑93 (MQ=255)

                                                       aGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAg                                  <  1:148118/93‑1 (MQ=255)

                                                       aGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAg                                  >  1:6264928/1‑93 (MQ=255)

                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 >  1:1625427/1‑93 (MQ=255)

                                                         gCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCTCATCCTAATTTTGCCGAGTGCAGCGATCAGACAGGTGGTGACGAAGAc                                >  1:320208/1‑93 (MQ=255)

                                                           ccATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgt                              <  2:3111321/93‑2 (MQ=255)

                                                            cATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgtc                             <  2:3892169/93‑3 (MQ=255)

                                                             aTAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcg                            >  1:5451344/1‑93 (MQ=255)

                                                             aTAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcg                            >  1:2447824/1‑93 (MQ=255)

                                                             aTAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcg                            >  2:1630292/1‑93 (MQ=255)

                                                             aTAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcg                            <  2:5544791/93‑1 (MQ=255)

                                                             aTAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcg                            >  1:324586/1‑93 (MQ=255)

                                                             aTAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcg                            >  1:4391109/1‑93 (MQ=255)

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

                                                               aGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCt                          <  2:3607796/93‑1 (MQ=255)

                                                                gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                         <  2:2103841/93‑1 (MQ=255)

                                                                gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                         >  1:6128816/1‑93 (MQ=255)

                                                                gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                         >  2:4983273/1‑93 (MQ=255)

                                                                    gTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTAt                     >  2:3433591/1‑93 (MQ=255)

                                                                       ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                  >  2:1954045/1‑93 (MQ=255)

                                                                       ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                  >  2:4199372/1‑93 (MQ=255)

                                                                       ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                  >  2:2484701/1‑93 (MQ=255)

                                                                        cAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCAt                 >  2:1764384/1‑93 (MQ=255)

                                                                        cAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCAt                 >  2:4039736/1‑93 (MQ=255)

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

                                                                                  ggCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGGAGTGACGAAGACGTCGCTGGGATCCACGCCAGCAACg       >  1:841573/1‑93 (MQ=255)

                                                                                   gCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGc      <  1:197042/93‑1 (MQ=255)

                                                                                   gCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGc      >  1:3719146/1‑93 (MQ=255)

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

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

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

                                                                                      ggTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCaa   >  2:45589/1‑93 (MQ=255)

                                                                                       gTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAg  <  1:2292888/93‑1 (MQ=255)

                                                                                       gTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAg  >  2:573346/1‑93 (MQ=255)

GACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAG  >  NC_000913/3893903‑3894082

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