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	99.6%	782.9 / ‑6.8	249	W222R (TGG→CGG)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (134/113); total (135/114)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

TATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCGCAG  >  NC_000913/3822715‑3822894

   ccGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCTGACTCGCGGGct                                                                                      >  2:1370653/1‑93 (MQ=255)

     ggTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctac                                                                                    <  2:3038663/93‑1 (MQ=255)

      gTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAACCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctact                                                                                   <  1:2514555/93‑1 (MQ=255)

       tGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctacta                                                                                  >  1:348779/1‑93 (MQ=255)

       tGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGccacta                                                                                  >  1:3901975/1‑93 (MQ=255)

        gCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTAt                                                                                 <  2:345961/93‑1 (MQ=255)

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

           ggCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTAGCGGGCTACTATGCa                                                                              <  2:1662798/93‑1 (MQ=255)

              gAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGc                                                                           >  1:623296/1‑93 (MQ=255)

              gAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGc                                                                           >  1:2830466/1‑93 (MQ=255)

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

              gAAAGTCACAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAATCTGATGCCGACTGGCGGGCTACTATGCAGGc                                                                           <  2:611655/93‑1 (MQ=255)

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

               aaaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGAGGGCTACTATGCAGGcg                                                                          >  1:853447/1‑92 (MQ=255)

                 aGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGc                                                                        <  2:1078403/93‑1 (MQ=255)

                  gTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGcc                                                                       >  2:3370693/1‑92 (MQ=255)

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

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

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

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

                       aaaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTa                                                                  >  2:3096281/1‑93 (MQ=255)

                       aaaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGAAGGGTa                                                                  >  2:2928619/1‑93 (MQ=255)

                          aCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAgg                                                               >  1:1069398/1‑92 (MQ=255)

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

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

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

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

                            gCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGc                                                             <  1:224248/93‑1 (MQ=255)

                            gCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGGGc                                                             >  1:3165224/1‑93 (MQ=255)

                             cAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGGGcc                                                            >  2:340486/1‑93 (MQ=255)

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

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

                              aaaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCa                                                           <  1:187341/93‑1 (MQ=255)

                              aaaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCa                                                           >  1:3408517/1‑93 (MQ=255)

                               aaGGCCTGCCTTATCTCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCCACTATGCAGGCTGGAGGGCAAGTGCCAt                                                          >  2:120118/1‑93 (MQ=255)

                               aaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCAt                                                          >  1:3308770/1‑93 (MQ=255)

                               aaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGGCAt                                                          >  2:2008965/1‑93 (MQ=255)

                                aGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATg                                                         >  2:1784280/1‑93 (MQ=255)

                                 ggCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATTCAGGCTGCAGGGTAAGTGCCATGa                                                        >  2:2734327/1‑93 (MQ=255)

                                  gcctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGaa                                                       <  2:3294604/93‑1 (MQ=255)

                                  gcctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGaa                                                       >  1:3328910/1‑93 (MQ=255)

                                    ctgactTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAg                                                     <  2:1012142/93‑1 (MQ=255)

                                     tgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGGAGGGTAAGTGCCATGAAAgg                                                    >  1:1279822/1‑93 (MQ=255)

                                     tgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAgg                                                    >  2:3199163/1‑93 (MQ=255)

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

                                      gcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGt                                                   >  2:3648849/1‑93 (MQ=255)

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

                                         ttATCCCCAGGTCAATCAGCAAGGTGAGATTGTATCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTATGTGCCATTAATGGTctc                                                >  1:4682001/1‑91 (MQ=255)

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

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

                                          tATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTTGcc                                               >  1:2437368/1‑93 (MQ=255)

                                           aTCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCt                                              <  2:1445955/93‑1 (MQ=255)

                                            tCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGTCATGAAAGGTCGCCTg                                             >  2:4319726/1‑93 (MQ=255)

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

                                              cccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGGGCCATGAAAGGTCGCCTGtt                                           >  1:3797870/1‑93 (MQ=255)

                                              cccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTCGCCTGtt                                           >  2:3992433/1‑93 (MQ=255)

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

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

                                               ccACGGCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTa                                          <  2:2096773/93‑1 (MQ=255)

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

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

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

                                                    tCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGc                                     >  2:1411712/1‑93 (MQ=255)

                                                     cAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCt                                    >  2:2136680/1‑93 (MQ=255)

                                                     cAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCt                                    <  1:4194147/93‑1 (MQ=255)

                                                      aaTCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTg                                   >  2:1934618/1‑93 (MQ=255)

                                                      aaTCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTg                                   >  1:2206428/1‑93 (MQ=255)

                                                      aaTCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGGCATGAAAGGTCGCCTGTTAGATGcgg                                   >  1:3722560/1‑91 (MQ=255)

                                                       aTCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGt                                  <  2:2489416/93‑1 (MQ=255)

                                                        tCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTc                                 >  1:4326282/1‑93 (MQ=255)

                                                        tCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTc                                 >  2:1749317/1‑93 (MQ=255)

                                                          aGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTccc                               <  1:3571742/93‑1 (MQ=255)

                                                          aGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGCGCCATGGAAGGTCGCCTGTTAGATGCTGTccc                               >  1:582319/1‑93 (MQ=255)

                                                           gCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCa                              <  2:3777897/93‑1 (MQ=255)

                                                           gCAAGGCGAGATCGAAGCTGAGGCCGACTGGCGGGGTACTAAGCAGGCTGCAGGGTAAGTGCCATGAAAGGGCGCCCGTTTGATGCTGTCCCa                              >  2:4560540/1‑93 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGGAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCAccc                           >  1:4216144/1‑91 (MQ=255)

                                                               ggagagaTCGAAGCTGCGGCCTGCTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGCCGCCTGTTAGATGCGGTCCCACTCa                          <  2:1627516/90‑1 (MQ=255)

                                                               ggCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGGAGGGGAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCa                          >  1:1821165/1‑93 (MQ=255)

                                                                   agaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTc                      >  1:2368437/1‑93 (MQ=255)

                                                                    gaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTcc                     >  2:2595157/1‑93 (MQ=255)

                                                                     aTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTccc                    >  1:2461945/1‑93 (MQ=255)

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

                                                                         aaGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAAc                >  1:2745374/1‑93 (MQ=255)

                                                                         aaGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAAc                >  1:2137863/1‑93 (MQ=255)

                                                                          aGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGAGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACg               >  2:2029538/1‑93 (MQ=255)

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

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

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

                                                                              gATGCCGACTGGCGGGCTACTATGCAGGGTGTAGGGTAAGTGCCATGGAATGTCGCCTGTTAGATGCTGTCCCACTCAGCTCCCTCACCGGCg           >  2:253318/1‑93 (MQ=255)

                                                                              gATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCg           >  2:2216338/1‑93 (MQ=255)

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

                                                                                  ccGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTgg       >  1:2076289/1‑93 (MQ=255)

                                                                                   cGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgc      <  1:4462473/93‑1 (MQ=255)

                                                                                   cGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgc      <  1:1404470/93‑1 (MQ=255)

                                                                                   cGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTACCTAACGGGCGTTggg      >  1:2108508/1‑92 (MQ=255)

                                                                                    gACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGGCGCCTGTTAGATGCTGTCCCCCTCCGTTCCCTAACGGGCGTTGgcg     >  2:856793/1‑93 (MQ=255)

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

                                                                                      cTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgca   <  1:2603974/93‑1 (MQ=255)

                                                                                      cTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTCTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgca   <  1:3296396/93‑1 (MQ=255)

                                                                                       tggcggGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgcag  >  2:594253/1‑93 (MQ=255)

TATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCGCAG  >  NC_000913/3822715‑3822894

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