BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,188,890	A→G	L37P (CTG→CCG)	phoQ ←	PhoQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,188,890	0	A	G	85.0%	544.8 / 91.2	220	L37P (CTG→CCG)	phoQ	PhoQ
Reads supporting (aligned to +/- strand): ref base A (12/21); new base G (112/75); total (124/96)
Fisher's exact test for biased strand distribution p-value = 1.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

Only 100 of 232 total aligned reads displayed.

TGTTTTCCCACTTCGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCGGGAAAAAA  >  NC_000913/1188798‑1188979

tGTTTTCCCCCTTCGCAAGGGTAAAGAACAGGTTGGCCTCCCCCCGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCa                                                                                           <  2:2374032/93‑1 (MQ=255)

 gTTTTCCCACTTCGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCAg                                                                                          <  1:5533109/93‑1 (MQ=255)

   tttCCCACTTCGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGgcg                                                                                        <  1:5163133/93‑1 (MQ=255)

    ttCCCACTTCGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGgcgc                                                                                       <  1:3821834/93‑1 (MQ=255)

      cccACTTCGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGa                                                                                     >  1:2879576/1‑93 (MQ=255)

         aCTTCGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCa                                                                                  >  2:4321765/1‑93 (MQ=255)

           ttCGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCAtt                                                                                >  1:4539311/1‑93 (MQ=255)

            tCGCAAGGTTATAGAACAGATTGCTCTCCCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTc                                                                               <  1:1979425/93‑1 (MQ=255)

            tCGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTc                                                                               <  1:5427081/93‑1 (MQ=255)

            tCGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTc                                                                               >  1:4234417/1‑93 (MQ=255)

             cGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTcc                                                                              >  1:788520/1‑93 (MQ=255)

             cGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTcc                                                                              >  1:433014/1‑93 (MQ=255)

             cGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTcc                                                                              >  1:229954/1‑93 (MQ=255)

              gCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCg                                                                             >  1:5833613/1‑93 (MQ=255)

              gCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCg                                                                             >  1:3047277/1‑93 (MQ=255)

               cAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGt                                                                            <  1:1335142/93‑1 (MQ=255)

               cAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGt                                                                            >  2:3164588/1‑93 (MQ=255)

               cAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGt                                                                            >  2:4824273/1‑93 (MQ=255)

               cAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACAATTCCGt                                                                            <  1:6538147/93‑1 (MQ=255)

                aaGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTa                                                                           <  2:6280402/93‑1 (MQ=255)

                aaGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTa                                                                           <  1:4740055/93‑1 (MQ=255)

                 aGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCAGCGCGACCATTCCGTAg                                                                          >  2:5205010/1‑93 (MQ=255)

                  gggTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCAGCGCGACCATTCCGTAgg                                                                         >  2:5344812/1‑93 (MQ=255)

                   ggTATAGAAAAGATTGCGCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGc                                                                        <  2:5166576/93‑1 (MQ=255)

                    gTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCa                                                                       <  1:5791673/93‑1 (MQ=255)

                    gTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGCCa                                                                       >  2:2016719/1‑93 (MQ=255)

                    gTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCa                                                                       <  1:4904173/93‑1 (MQ=255)

                     tataGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCaa                                                                      <  2:682793/93‑1 (MQ=255)

                         gAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGa                                                                  >  1:1807168/1‑93 (MQ=255)

                           aCAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGaaa                                                                <  1:6629315/93‑1 (MQ=255)

                               aTTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCAc                                                            >  2:271090/1‑93 (MQ=255)

                                   ctctCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGt                                                        >  1:1450754/1‑93 (MQ=255)

                                    tctcGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGta                                                       >  1:2662469/1‑93 (MQ=255)

                                      tcGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGtact                                                     >  2:5899554/1‑93 (MQ=255)

                                        gCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGtactac                                                   <  2:4516983/93‑1 (MQ=255)

                                         ccACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTAcc                                                  <  1:1053170/93‑1 (MQ=255)

                                          cACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCg                                                 >  2:5882779/1‑93 (MQ=255)

                                           aCGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGc                                                <  2:6313726/93‑1 (MQ=255)

                                           aCGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGc                                                <  1:657258/93‑1 (MQ=255)

                                            cGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCt                                               >  1:1679211/1‑93 (MQ=255)

                                            cGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCt                                               >  1:3348453/1‑93 (MQ=255)

                                              tAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGc                                             >  1:5112805/1‑93 (MQ=255)

                                               aaCAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGcc                                            >  1:4042477/1‑93 (MQ=255)

                                                aCAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCg                                           >  2:1054259/1‑93 (MQ=255)

                                                aCAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCg                                           >  1:645636/1‑93 (MQ=255)

                                                  aGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGtt                                         >  2:2081944/1‑93 (MQ=255)

                                                   gCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTg                                        >  2:4375172/1‑93 (MQ=255)

                                                   gCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTg                                        >  2:6340519/1‑93 (MQ=255)

                                                   gCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTg                                        >  1:5680139/1‑93 (MQ=255)

                                                   gCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTg                                        >  2:2254036/1‑93 (MQ=255)

                                                    ccGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGc                                       <  1:4187309/93‑1 (MQ=255)

                                                      gAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGTAAGCACCAGTACTACCGCTGCCGTTGCCa                                     >  1:1852044/1‑93 (MQ=255)

                                                      gAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCa                                     >  2:2094619/1‑93 (MQ=255)

                                                       aaaCGTAGTTTTATCGAAACTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCaa                                    <  2:1947948/93‑1 (MQ=255)

                                                        accGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAAc                                   <  2:6166429/91‑1 (MQ=255)

                                                        aaCGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAAc                                   <  1:658456/93‑1 (MQ=255)

                                                          cGTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAg                                 >  1:5560237/1‑93 (MQ=255)

                                                           gTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGa                                <  2:4707222/93‑1 (MQ=255)

                                                           gTAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGa                                >  1:5076670/1‑93 (MQ=255)

                                                            tAGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGaa                               >  2:6426987/1‑93 (MQ=255)

                                                             aGTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGaaa                              <  2:1865815/93‑1 (MQ=255)

                                                              gTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGaaaa                             >  1:600452/1‑93 (MQ=255)

                                                              gTTTTATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGaaaa                             <  1:1539860/93‑1 (MQ=255)

                                                               ttttATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAAc                            <  1:6031497/93‑1 (MQ=255)

                                                               ttttATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAAc                            <  1:3269038/93‑1 (MQ=255)

                                                                tttATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACg                           <  1:2528788/93‑1 (MQ=255)

                                                                  tATCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTa                         >  1:3916993/1‑93 (MQ=255)

                                                                   aTCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTAc                        >  2:5034687/1‑93 (MQ=255)

                                                                   aTCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTAc                        >  1:3483100/1‑93 (MQ=255)

                                                                    tCGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTAcc                       <  2:4122769/93‑1 (MQ=255)

                                                                     cGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTAccc                      >  2:5489080/1‑93 (MQ=255)

                                                                     cGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTAccc                      <  2:5079658/93‑1 (MQ=255)

                                                                     cGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTAccc                      <  1:2613745/93‑1 (MQ=255)

                                                                     cGAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTAccc                      <  2:4758510/93‑1 (MQ=255)

                                                                     cGAAACTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTAccc                      <  2:2633806/93‑1 (MQ=255)

                                                                      gAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCg                     >  1:2792270/1‑93 (MQ=255)

                                                                      gAAACTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCg                     <  2:1802262/93‑1 (MQ=255)

                                                                      gAAACTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCg                     <  2:376759/93‑1 (MQ=255)

                                                                        aaCTGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCa                   >  2:3993125/1‑93 (MQ=255)

                                                                        aaCTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCa                   <  2:3323285/93‑1 (MQ=255)

                                                                          cTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGc                 >  2:1541744/1‑93 (MQ=255)

                                                                          cTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGc                 <  1:1229802/93‑1 (MQ=255)

                                                                           tGACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCg                >  2:5008733/1‑93 (MQ=255)

                                                                           cgACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCg                <  2:5769383/92‑1 (MQ=255)

                                                                            gACGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCga               >  2:4211498/1‑93 (MQ=255)

                                                                            gACGCTAGAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCga               >  1:5373790/1‑93 (MQ=255)

                                                                             aCGCTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCgag              <  2:878865/93‑1 (MQ=255)

                                                                                cTATAACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGc           <  1:3510904/93‑1 (MQ=255)

                                                                                cTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGc           >  2:1266225/1‑93 (MQ=255)

                                                                                cTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGc           <  1:4323611/93‑1 (MQ=255)

                                                                                 tataACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCg          <  2:563627/93‑1 (MQ=255)

                                                                                 tataACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCg          >  2:6367212/1‑93 (MQ=255)

                                                                                  ataACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCgg         <  1:3441874/93‑1 (MQ=255)

                                                                                   taACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCggg        <  1:497414/93‑1 (MQ=255)

                                                                                   taACCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGAACCAGTACTACAGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCggg        <  2:6655156/93‑1 (MQ=255)

                                                                                    aaCCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCGGGa       >  1:685294/1‑93 (MQ=255)

                                                                                     aCCGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCGGGaa      >  2:5386863/1‑93 (MQ=255)

                                                                                       cGATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCGGGaaaa    >  1:1186488/1‑93 (MQ=255)

                                                                                        gATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCGGGaaaaa   >  1:6267953/1‑93 (MQ=255)

                                                                                        gATCGGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCGGGaaaaa   >  1:6224044/1‑93 (MQ=255)

TGTTTTCCCACTTCGCAAGGGTATAGAACAGATTGCTCTCGCCACGTAACAGCCGAAACGTAGTTTTATCGAAACTGACGCTATAACCGATCAGCGCGACCATTCCGTAGGCAAGCGAAAGCACCAGTACTACCGCTGCCGTTGCCAACAGAAAACGTACCCGCAGCGAGAGCGGGAAAAAA  >  NC_000913/1188798‑1188979

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