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.5%	1213.0 / ‑3.6	377	W222R (TGG→CGG)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (171/204); total (173/204)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 395 total aligned reads displayed.

TATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCGCAG  >  NC_000913/3822715‑3822894

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

  tCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGc                                                                                       <  2:5358261/93‑1 (MQ=255)

   ccGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGct                                                                                      >  1:4119669/1‑93 (MQ=255)

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

         cTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATg                                                                                >  1:6222047/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                   tCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAg                                                                      <  2:1551278/93‑1 (MQ=255)

                    cGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAgg                                                                     <  1:6985097/93‑1 (MQ=255)

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

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

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

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

                         aaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAg                                                                >  1:5852745/1‑93 (MQ=255)

                         aaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGAAGGGAAAg                                                                >  2:3176807/1‑93 (MQ=255)

                          aCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGt                                                               >  2:4302491/1‑93 (MQ=255)

                            gCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGc                                                             >  1:3531342/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                    ctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAg                                                     >  2:2337998/1‑93 (MQ=255)

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

                                       cctTATCCCCACGTCAATCAGCAAGGCGATATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTc                                                  <  2:9826623/93‑1 (MQ=255)

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

                                       cctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTc                                                  >  1:7668246/1‑93 (MQ=255)

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

                                        ctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTCg                                                 >  2:2956683/1‑93 (MQ=255)

                                         ttATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGc                                                <  1:3157331/93‑1 (MQ=255)

                                         ttATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTCGc                                                >  2:118131/1‑93 (MQ=255)

                                         ttATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAAGGTAAGTGCCATGAAAGGTAGc                                                >  1:6770/1‑93 (MQ=255)

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

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

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

                                           ttCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGATACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCt                                              <  1:5795074/92‑1 (MQ=255)

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

                                           aTCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTCGCCt                                              >  1:2710675/1‑93 (MQ=255)

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

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

                                                cACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCGAGGTAAGTGCCCTGAAAGGTCGCCTGTTAg                                         >  1:9616062/1‑93 (MQ=255)

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

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

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

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

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

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

                                                    tCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGTCATGAAAGGGCGCCTGTTAGATgt                                     >  2:864419/1‑92 (MQ=255)

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

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

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

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

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

                                                         cAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGTTGCTGTcc                                >  2:4532568/1‑93 (MQ=255)

                                                         cAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTcc                                >  2:8074150/1‑93 (MQ=255)

                                                         cAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTcc                                >  1:1871720/1‑93 (MQ=255)

                                                           gCAAGGCGAGTTCGAAGCTGATGCTGACTGGCGGGCTACTATGCAGGCTGCTGGGGAAGTGTCTAGTATGGTCGCTTGTTAGATGCTTTCCCa                              >  2:8933797/1‑93 (MQ=255)

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

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

                                                             aaGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACt                            >  1:1744763/1‑93 (MQ=255)

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

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

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

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                           <  1:7370405/93‑1 (MQ=255)

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

                                                                 cGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGt                        <  2:1295134/93‑1 (MQ=255)

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

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

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

                                                                     aTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGGGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTccc                    >  1:6846281/1‑93 (MQ=255)

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

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

                                                                        gAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTaa                 <  2:4951154/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                              gATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCg           <  2:2179642/93‑1 (MQ=255)

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

                                                                              gATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCg           <  1:2963019/93‑1 (MQ=255)

                                                                               aTGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGt          >  1:2606518/1‑93 (MQ=255)

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

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

                                                                                tGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGtt         >  1:6547488/1‑93 (MQ=255)

                                                                                 gCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTg        <  2:2684051/93‑1 (MQ=255)

                                                                                 gCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTg        <  2:7721922/93‑1 (MQ=255)

                                                                                 gCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTg        <  1:1342908/93‑1 (MQ=255)

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

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

                                                                                    gACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgcg     <  2:1803085/93‑1 (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` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`