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.8%	1314.7 / ‑5.5	416	W222R (TGG→CGG)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (192/223); total (192/224)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 448 total aligned reads displayed.

TATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCGCAG  >  NC_000913/3822715‑3822894

  tCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGc                                                                                       >  1:4602935/1‑93 (MQ=255)

  tCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGc                                                                                       >  1:2291196/1‑93 (MQ=255)

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

       tGCTGTCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctacta                                                                                  <  1:3960723/93‑1 (MQ=255)

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

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

        gCTGGCGAAAGCCGCAAAACGCAAAGGCCTGCCTTATCCCGACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTAt                                                                                 <  2:2330114/93‑1 (MQ=255)

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

              gAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGCCTGGCGGGCTACTATGCAGGc                                                                           <  1:1414128/93‑1 (MQ=255)

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

                aaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCgggggg                                                                         >  1:5821206/1‑90 (MQ=255)

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

                aaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTg                                                                         <  2:6763759/93‑1 (MQ=255)

                 agccGCAAAACGCAAAGGCCTGCCTTCTCCCCACGTCAATCAGCAAGGCGAGATCTAAGCTGATGACGCCTGGCGGGCTACTATGCAGGCTGc                                                                        <  1:5690548/90‑1 (MQ=255)

                 aGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGc                                                                        >  2:2181547/1‑93 (MQ=255)

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

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

                  gTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCa                                                                       >  1:6937425/1‑93 (MQ=255)

                   tCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAg                                                                      >  2:2040538/1‑93 (MQ=255)

                   tCGCAAAACGCAAACGCCTGCTTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAg                                                                      <  1:2033795/93‑1 (MQ=255)

                    cGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAgg                                                                     >  2:1362334/1‑93 (MQ=255)

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

                     gCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAggg                                                                    >  2:4154157/1‑93 (MQ=255)

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

                         aaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAg                                                                <  1:1814279/93‑1 (MQ=255)

                         aaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAg                                                                <  2:7416842/93‑1 (MQ=255)

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

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

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

                             cAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAATCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGcc                                                            <  2:4678229/93‑1 (MQ=255)

                             cAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTggc                                                            >  2:7109489/1‑91 (MQ=255)

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

                              aaaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTCGCGTGCTACCATGCCGGCTGCATCGTCAGTGACa                                                           >  2:567053/1‑93 (MQ=255)

                               aaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCAt                                                          <  2:4080733/93‑1 (MQ=255)

                                aGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATg                                                         <  1:3705344/93‑1 (MQ=255)

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

                                   cctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGaaa                                                      >  1:4499931/1‑93 (MQ=255)

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

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

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

                                    ctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCGACTATGCAGGCTGAAGGGGAAGCGCCCTGgaaa                                                     >  1:3113331/1‑92 (MQ=255)

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

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

                                     tgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAgg                                                    >  1:86025/1‑93 (MQ=255)

                                     tgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAgg                                                    >  1:1185206/1‑93 (MQ=255)

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

                                        ctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCg                                                 <  1:6480620/93‑1 (MQ=255)

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

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

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

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

                                             ccccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGt                                            >  2:6813720/1‑93 (MQ=255)

                                              cccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGtt                                           <  2:1494568/93‑1 (MQ=255)

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

                                                cACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAg                                         >  1:304056/1‑93 (MQ=255)

                                                 gcgacAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                        <  2:5862543/89‑1 (MQ=255)

                                                 ccGTCAATCAGCAAGGCGAGATTGAAGCTGATGCCGACGGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                        <  1:2517515/92‑1 (MQ=255)

                                                 aCGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                        >  2:764103/1‑93 (MQ=255)

                                                 aCGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGGGCCATGCAAAGTCGCCTGTTAGa                                        >  2:4983117/1‑93 (MQ=255)

                                                   gTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGGGCCATGAAAGGTCGCCTGTTAGATg                                      >  2:5591814/1‑93 (MQ=255)

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

                                                          aGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGCAAGGTCGCCTGTTAGATGCTGTccc                               >  1:2863959/1‑93 (MQ=255)

                                                           gCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTTGAAGCTGTcccc                              >  2:7043436/1‑92 (MQ=255)

                                                           gCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCa                              >  1:900974/1‑93 (MQ=255)

                                                           gCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCa                              >  2:7347765/1‑93 (MQ=255)

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

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

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

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGGAAGGTGCGCTGTTAGATGGTGTCCCagta                           >  2:1082553/1‑90 (MQ=255)

                                                                gCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAg                         <  2:6139356/93‑1 (MQ=255)

                                                                gCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAg                         <  1:4351793/93‑1 (MQ=255)

                                                                  gagaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGtt                       <  1:970416/93‑1 (MQ=255)

                                                                  gagaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGtt                       <  1:2659665/93‑1 (MQ=255)

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

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

                                                                    gaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTcc                     <  2:3563733/93‑1 (MQ=255)

                                                                     aTCGAAGCTGGTGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTccc                    <  2:1852552/93‑1 (MQ=255)

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

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

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

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

                                                                        gAAGATGATGCCGACTGGCGGGCTACTATGCAGGCTGCTGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTaa                 >  1:7837587/1‑93 (MQ=255)

                                                                         aaGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAAc                <  1:3601496/93‑1 (MQ=255)

                                                                         aaGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAAc                <  2:5286754/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                               aTGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGt          <  2:2666733/93‑1 (MQ=255)

                                                                                  ccGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTgg       <  1:5419131/93‑1 (MQ=255)

                                                                                  ccGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTgg       <  1:2619857/93‑1 (MQ=255)

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

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

                                                                                   cGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgc      >  2:490656/1‑93 (MQ=255)

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

                                                                                      cTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgca   >  2:830292/1‑93 (MQ=255)

                                                                                       tggcggGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgcag  >  1:7412116/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` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`