BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,822,805	T→C	W222R (TGG→CGG)	spoU →	putative RNA methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,822,805	0	T	C	100.0%	837.8 / NA	257	W222R (TGG→CGG)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (135/122); total (135/122)

Only 100 of 267 total aligned reads displayed.

TATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCGCAG  >  NC_000913/3822715‑3822894

 aTCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCggg                                                                                        <  2:1716232/93‑1 (MQ=255)

      gTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctact                                                                                   <  2:931820/93‑1 (MQ=255)

      gTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctact                                                                                   >  2:5726801/1‑93 (MQ=255)

      gTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctact                                                                                   >  1:3402479/1‑93 (MQ=255)

      gTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctact                                                                                   >  2:4204607/1‑93 (MQ=255)

      gTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctact                                                                                   <  1:2392468/93‑1 (MQ=255)

       tGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctacta                                                                                  <  1:2161205/93‑1 (MQ=255)

          tGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGc                                                                               <  2:5113998/93‑1 (MQ=255)

          tGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGc                                                                               <  1:3815855/93‑1 (MQ=255)

          tGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGc                                                                               >  2:2555118/1‑93 (MQ=255)

          tGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGc                                                                               <  1:6273269/93‑1 (MQ=255)

          tGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGc                                                                               <  1:4091063/93‑1 (MQ=255)

           ggCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCa                                                                              >  1:1195795/1‑93 (MQ=255)

           ggCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCa                                                                              <  1:5809693/93‑1 (MQ=255)

              gAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGc                                                                           <  1:305468/93‑1 (MQ=255)

              gAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCACGc                                                                           >  2:3311352/1‑93 (MQ=255)

               aaaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGGAGGCt                                                                          >  2:4785614/1‑93 (MQ=255)

               aaaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCt                                                                          >  2:2237445/1‑93 (MQ=255)

                aaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTg                                                                         >  1:997000/1‑93 (MQ=255)

                  gTCGCAAAACGCAAGGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAAGCTGCa                                                                       >  1:2132007/1‑93 (MQ=255)

                  gTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCa                                                                       <  1:2755975/93‑1 (MQ=255)

                  gTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCa                                                                       <  1:1822768/93‑1 (MQ=255)

                     gCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAggg                                                                    <  2:6153464/93‑1 (MQ=255)

                      cAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGt                                                                   <  2:6134475/93‑1 (MQ=255)

                      cAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGt                                                                   <  1:4425613/93‑1 (MQ=255)

                       aaaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTa                                                                  <  2:3230616/93‑1 (MQ=255)

                       aaaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTa                                                                  <  1:1365376/93‑1 (MQ=255)

                        aaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTaa                                                                 <  1:2209918/93‑1 (MQ=255)

                        aaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTaa                                                                 <  1:4140017/93‑1 (MQ=255)

                          aCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTTGTATGCAGGCTGCAGGGTAAgc                                                               >  1:2262357/1‑92 (MQ=255)

                          aCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGt                                                               <  2:6018562/93‑1 (MQ=255)

                          aCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGATACTATGCAGGATTCAGGATAAga                                                               >  1:109361/1‑92 (MQ=255)

                           cGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTg                                                              <  1:3224909/93‑1 (MQ=255)

                              aaaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCa                                                           >  2:5879980/1‑93 (MQ=255)

                              aaaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCa                                                           >  1:5164767/1‑93 (MQ=255)

                                   cctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGaaa                                                      <  1:4809436/93‑1 (MQ=255)

                                    ctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAg                                                     <  1:276806/93‑1 (MQ=255)

                                    ctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAg                                                     >  1:5707230/1‑93 (MQ=255)

                                     tgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAgg                                                    <  2:4365552/93‑1 (MQ=255)

                                     tgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGGGCCATGAAAgg                                                    >  2:4555067/1‑93 (MQ=255)

                                      gcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTGAGCGGCATGAAAGGt                                                   >  2:3424705/1‑93 (MQ=255)

                                      gcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGt                                                   <  2:5362453/93‑1 (MQ=255)

                                      gcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGt                                                   >  1:2330049/1‑93 (MQ=255)

                                       cctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTc                                                  <  2:3061801/93‑1 (MQ=255)

                                       cctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTc                                                  <  1:3547764/93‑1 (MQ=255)

                                        ctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCg                                                 >  2:4761818/1‑93 (MQ=255)

                                        ctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCg                                                 >  1:611057/1‑93 (MQ=255)

                                         ttATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGc                                                >  2:3871553/1‑93 (MQ=255)

                                          tagccccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGcc                                               <  2:2388130/90‑1 (MQ=255)

                                          tATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGcc                                               >  1:2901794/1‑93 (MQ=255)

                                          tATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGcc                                               >  2:1526138/1‑93 (MQ=255)

                                          tATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGcc                                               <  1:5537207/93‑1 (MQ=255)

                                          tATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGcc                                               >  1:6059156/1‑93 (MQ=255)

                                            tCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTg                                             <  2:3668971/93‑1 (MQ=255)

                                              cccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGGCTGtt                                           >  1:4266210/1‑93 (MQ=255)

                                              cccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGtt                                           >  1:1524188/1‑93 (MQ=255)

                                              cccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGtt                                           >  2:3606375/1‑93 (MQ=255)

                                               ccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTa                                          <  2:693295/93‑1 (MQ=255)

                                               ccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTa                                          >  1:2511871/1‑93 (MQ=255)

                                                 aCGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                        <  2:53018/93‑1 (MQ=255)

                                                 aCGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                        <  1:5268437/93‑1 (MQ=255)

                                                  cGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGAt                                       >  1:4625610/1‑93 (MQ=255)

                                                  cGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGAt                                       <  1:682042/93‑1 (MQ=255)

                                                   gTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATg                                      >  2:4987676/1‑93 (MQ=255)

                                                    tCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGGCGCCTGTTAGATGc                                     >  2:3923071/1‑93 (MQ=255)

                                                        tCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTc                                 <  2:3947163/93‑1 (MQ=255)

                                                            cAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCAc                             <  1:2679906/93‑1 (MQ=255)

                                                             aaGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACt                            <  1:4050704/93‑1 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                           >  2:3016231/1‑93 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                           >  2:484588/1‑93 (MQ=255)

                                                               ggCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTGAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCa                          >  2:6011336/1‑93 (MQ=255)

                                                               ggCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCa                          <  2:4096551/93‑1 (MQ=255)

                                                                  gagaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGtt                       >  1:5728895/1‑93 (MQ=255)

                                                                  gagaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTCTCCCACTCAGtt                       >  2:2066799/1‑93 (MQ=255)

                                                                     aTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTccc                    <  2:1331070/93‑1 (MQ=255)

                                                                      tCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCt                   <  1:1009420/93‑1 (MQ=255)

                                                                      tCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCt                   >  1:1598305/1‑93 (MQ=255)

                                                                       cGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTa                  <  2:5382203/93‑1 (MQ=255)

                                                                       cGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCACTa                  >  2:5750092/1‑93 (MQ=255)

                                                                          aGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACg               >  2:3138374/1‑93 (MQ=255)

                                                                          aGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACg               >  2:183385/1‑93 (MQ=255)

                                                                           gCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACgg              <  1:4200543/93‑1 (MQ=255)

                                                                           gCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACgg              <  1:799912/93‑1 (MQ=255)

                                                                           gCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACgg              >  2:2199944/1‑93 (MQ=255)

                                                                            cTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACggg             <  2:1184162/93‑1 (MQ=255)

                                                                            cTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACggg             >  2:1956867/1‑93 (MQ=255)

                                                                            cTGATGCCGACTGGCGGGCTACTATGCACGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACggg             >  1:4952213/1‑93 (MQ=255)

                                                                             tGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGc            >  1:2025422/1‑93 (MQ=255)

                                                                              gATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCCCTCAGTTCCCTAACGGGCg           >  1:3584461/1‑93 (MQ=255)

                                                                               aTGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGCAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTTCCTAACCGGCGt          >  2:2488175/1‑93 (MQ=255)

                                                                                tGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGGCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGtt         <  2:2777641/93‑1 (MQ=255)

                                                                                tGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGtt         <  1:1995479/93‑1 (MQ=255)

                                                                                 gCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTg        >  1:4763826/1‑93 (MQ=255)

                                                                                 gCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTg        >  1:1803538/1‑93 (MQ=255)

                                                                                 gCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAAGTCGCCTGTTAGATGCTGTCCCACTCTGTTCCCTAACGGGCGGTg        >  2:5581341/1‑93 (MQ=255)

                                                                                    gACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgcg     >  1:543738/1‑93 (MQ=255)

                                                                                     aCTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgcgc    >  2:1663051/1‑93 (MQ=255)

                                                                                     aCTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgcgc    <  1:5014938/93‑1 (MQ=255)

                                                                                       tggcggGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgcag  <  2:1155447/93‑1 (MQ=255)

                                                                                       tggcggGCAACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgcag  >  2:6115355/1‑93 (MQ=255)

TATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCGCAG  >  NC_000913/3822715‑3822894

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