BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,351,962	C→T	R194H (CGT→CAT)	sapC ←	SapC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,351,962	0	C	T	99.6%	840.1 / ‑5.7	233	R194H (CGT→CAT)	sapC	SapC
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (121/111); total (121/112)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 237 total aligned reads displayed.

AGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACA  >  NC_000913/1351872‑1352051

aGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATgg                                                                                         <  2:928894/93‑1 (MQ=255)

aGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATgg                                                                                         <  2:2465918/93‑1 (MQ=255)

aGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATgg                                                                                         >  2:5169356/1‑93 (MQ=255)

 gTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATggg                                                                                        <  1:7858640/93‑1 (MQ=255)

 gTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATggg                                                                                        >  2:4883871/1‑93 (MQ=255)

 gTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATggg                                                                                        >  1:283007/1‑93 (MQ=255)

  tGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGggc                                                                                       >  1:933428/1‑93 (MQ=255)

   gcgcGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGggcg                                                                                      <  2:1129300/93‑1 (MQ=255)

       gggTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGa                                                                                  <  2:1780569/93‑1 (MQ=255)

        ggTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGAt                                                                                 >  1:6902432/1‑93 (MQ=255)

        ggTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGAt                                                                                 >  2:4654470/1‑93 (MQ=255)

         gTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATa                                                                                <  1:551180/93‑1 (MQ=255)

         gTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATa                                                                                <  2:1502492/93‑1 (MQ=255)

          tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATaa                                                                               >  1:5899163/1‑93 (MQ=255)

          tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATaa                                                                               <  2:6428662/93‑1 (MQ=255)

          tGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATaa                                                                               >  2:3897819/1‑93 (MQ=255)

          tGATCTCGGTGACCAGCCCGGCGCTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATaa                                                                               <  2:4304430/93‑1 (MQ=255)

             tctcGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGt                                                                            <  1:2508075/93‑1 (MQ=255)

              ctcGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTa                                                                           <  1:2940661/93‑1 (MQ=255)

              ctcGGGCACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCTCCATCCAGATGGGCGGCGATAACGTa                                                                           <  2:789272/93‑1 (MQ=255)

               tcGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTAc                                                                          >  1:4473365/1‑93 (MQ=255)

                  gTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTActct                                                                       >  2:684732/1‑93 (MQ=255)

                  gTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTActct                                                                       >  2:7064341/1‑93 (MQ=255)

                   ggACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCtt                                                                      <  2:1845099/92‑1 (MQ=255)

                    gACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCttt                                                                     <  2:4027948/93‑1 (MQ=255)

                     aCCAGCCCGGCGGTGATGTTTGTCATCACAGAAAACAAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTTCTCtttt                                                                    <  2:7059614/93‑1 (MQ=255)

                       cAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTc                                                                  >  1:6058631/1‑93 (MQ=255)

                         gCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCa                                                                >  2:2660047/1‑93 (MQ=255)

                         gCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCa                                                                >  1:5715276/1‑93 (MQ=255)

                         gCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCa                                                                >  1:1140101/1‑93 (MQ=255)

                         gCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCa                                                                >  2:6219599/1‑93 (MQ=255)

                          cccGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAg                                                               <  2:1324940/93‑1 (MQ=255)

                           ccGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGt                                                              >  2:4861965/1‑93 (MQ=255)

                           ccGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGt                                                              >  2:1473982/1‑93 (MQ=255)

                            cggcggTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGtt                                                             >  2:399393/1‑93 (MQ=255)

                              gcggTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCg                                                           >  2:6101859/1‑93 (MQ=255)

                              gcggTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCg                                                           >  1:7621787/1‑93 (MQ=255)

                               cggTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGt                                                          <  2:2619246/93‑1 (MQ=255)

                                 gTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCa                                                        >  2:586491/1‑93 (MQ=255)

                                  tGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCAt                                                       >  2:1223378/1‑93 (MQ=255)

                                   gATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATg                                                      >  1:1562660/1‑93 (MQ=255)

                                     tGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCa                                                    <  2:3694914/93‑1 (MQ=255)

                                     tGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCa                                                    >  2:6797545/1‑93 (MQ=255)

                                      gTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCAc                                                   <  1:3972863/93‑1 (MQ=255)

                                            cATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCt                                             <  1:3220009/93‑1 (MQ=255)

                                            cATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCt                                             >  2:6690181/1‑93 (MQ=255)

                                             aTCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTg                                            <  1:6793394/93‑1 (MQ=255)

                                             aTCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTg                                            >  1:738088/1‑93 (MQ=255)

                                              tCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGt                                           >  2:573462/1‑93 (MQ=255)

                                              tCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGt                                           >  2:2951908/1‑93 (MQ=255)

                                              tCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGt                                           >  1:1002275/1‑93 (MQ=255)

                                              tCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGt                                           >  2:1905566/1‑93 (MQ=255)

                                              tCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGt                                           >  2:7351221/1‑93 (MQ=255)

                                               cacaGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTa                                          <  2:2851990/93‑1 (MQ=255)

                                                    cAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCg                                     >  2:3129083/1‑93 (MQ=255)

                                                    cAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCg                                     >  1:2604470/1‑93 (MQ=255)

                                                     aaaCCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGa                                    >  1:4240506/1‑93 (MQ=255)

                                                     aaaCCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGa                                    >  1:3826651/1‑93 (MQ=255)

                                                      aaCCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGaa                                   <  1:3780186/93‑1 (MQ=255)

                                                      aaCCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGaa                                   <  1:6442215/93‑1 (MQ=255)

                                                      aaCCAGAGAATATGCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGaa                                   <  2:4525650/93‑1 (MQ=255)

                                                       aCCAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAAc                                  <  2:5414575/93‑1 (MQ=255)

                                                        ccAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACg                                 >  1:5420267/1‑93 (MQ=255)

                                                         cAGAGAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGt                                >  1:2709805/1‑93 (MQ=255)

                                                          agagAATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTa                               >  1:558849/1‑93 (MQ=255)

                                                           gagaATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTAc                              >  2:1537095/1‑93 (MQ=255)

                                                           gagaATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTAc                              <  1:7235428/93‑1 (MQ=255)

                                                             taaTATTCAGCGTTGACGCGTCATCCTGATGGGCGTCGATAACATACTCTTTTCCCAGTTCATCATGCACCATGCATCAATTCGAACGTGCCa                            <  2:2508537/92‑1 (MQ=255)

                                                             gaATATTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCa                            >  1:1958139/1‑93 (MQ=255)

                                                               atatTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATa                          >  1:1983353/1‑93 (MQ=255)

                                                               atatTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATa                          <  1:4977672/93‑1 (MQ=255)

                                                               atatTCAGCGTTGATGCGCCATCCAGATGGACGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATa                          <  2:3581517/93‑1 (MQ=255)

                                                                tatTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATAc                         >  1:3088281/1‑93 (MQ=255)

                                                                tatTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATAc                         <  1:6642011/93‑1 (MQ=255)

                                                                tatTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATAc                         >  2:1955067/1‑93 (MQ=255)

                                                                 atTCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATAcg                        >  2:3206635/1‑93 (MQ=255)

                                                                  ttCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATAcgc                       >  1:3058207/1‑93 (MQ=255)

                                                                   tCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATAcgcg                      >  1:509063/1‑93 (MQ=255)

                                                                   tCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATAcgcg                      <  1:6371822/93‑1 (MQ=255)

                                                                   tCAGCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATAcgcg                      >  2:4914743/1‑93 (MQ=255)

                                                                      gCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACTATACGCGgcc                   >  2:7501116/1‑92 (MQ=255)

                                                                      gCGTTGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGgca                   <  1:4605867/93‑1 (MQ=255)

                                                                         ttGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGgcagca                <  1:1766584/93‑1 (MQ=255)

                                                                         ttGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGgcagca                >  1:4174816/1‑93 (MQ=255)

                                                                          tGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCaa               <  1:4412040/93‑1 (MQ=255)

                                                                          tGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCaa               <  2:7179211/93‑1 (MQ=255)

                                                                          tGATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCaa               <  2:7075908/93‑1 (MQ=255)

                                                                           gATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAAt              <  1:322965/93‑1 (MQ=255)

                                                                           gATGCGCCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAAt              >  1:5311804/1‑93 (MQ=255)

                                                                              gcgcCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATgcc           >  2:6712114/1‑93 (MQ=255)

                                                                              gcgcCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATgcc           <  1:5369128/93‑1 (MQ=255)

                                                                              gcgcCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATgcc           <  1:4789424/93‑1 (MQ=255)

                                                                               cgcCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATgcca          <  1:1131499/93‑1 (MQ=255)

                                                                                gcCATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATgccag         >  1:5976713/1‑93 (MQ=255)

                                                                                 ccATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATgccagc        <  1:2198217/93‑1 (MQ=255)

                                                                                 ccATCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATgccagc        >  1:1210116/1‑93 (MQ=255)

                                                                                   aTCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAAAgccagcca      >  2:5502345/1‑93 (MQ=255)

                                                                                    tCCAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGccagccag     >  1:1420878/1‑93 (MQ=255)

                                                                                      cAGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGAc   >  2:6195432/1‑93 (MQ=255)

                                                                                       aGATGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATgccagc        >  2:3085911/1‑87 (MQ=255)

AGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACA  >  NC_000913/1351872‑1352051

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