BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,215,945	T→C	S420P (TCT→CCT)	aceA →	isocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,215,945	0	T	C	99.6%	779.3 / ‑4.5	248	S420P (TCT→CCT)	aceA	isocitrate lyase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (148/99); minor base A (1/0); total (149/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 254 total aligned reads displayed.

GCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTC  >  NC_000913/4215855‑4216034

  gAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTc                                                                                       >  1:413417/1‑93 (MQ=255)

   aaaGATGGCTATACCTTCGTATCTCACCAGAAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCa                                                                                      <  2:1789715/93‑1 (MQ=255)

     aGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGt                                                                                    >  2:3320969/1‑93 (MQ=255)

     aGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGt                                                                                    <  1:2017076/93‑1 (MQ=255)

       aTGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCa                                                                                  >  1:1136320/1‑93 (MQ=255)

        tGGCTATCCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCAc                                                                                 <  2:4318621/93‑1 (MQ=255)

        tGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCAc                                                                                 >  2:1628363/1‑93 (MQ=255)

         ggCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCAcc                                                                                >  2:1299273/1‑93 (MQ=255)

            tataCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACcgcg                                                                             >  1:1367348/1‑93 (MQ=255)

             ataCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCgcgc                                                                            <  1:4393951/93‑1 (MQ=255)

             ataCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCgcgc                                                                            >  1:1253456/1‑93 (MQ=255)

              taCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCt                                                                           >  2:2753766/1‑93 (MQ=255)

              taCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCt                                                                           >  2:1228536/1‑93 (MQ=255)

               aCCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTg                                                                          >  1:2982653/1‑93 (MQ=255)

                ccTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGa                                                                         >  2:2221082/1‑93 (MQ=255)

                 cTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGAc                                                                        >  1:1719224/1‑93 (MQ=255)

                  ttCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGAcc                                                                       >  2:2841548/1‑93 (MQ=255)

                     gTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGc                                                                    >  2:4373210/1‑93 (MQ=255)

                       aTCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTc                                                                  >  2:3257499/1‑93 (MQ=255)

                       aTCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTc                                                                  >  2:3903401/1‑93 (MQ=255)

                        tctcACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTcc                                                                 >  1:1890845/1‑93 (MQ=255)

                            aCCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTg                                                             >  1:2246397/1‑93 (MQ=255)

                            aCCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTg                                                             >  2:2584915/1‑93 (MQ=255)

                             ccAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTga                                                            >  1:228951/1‑93 (MQ=255)

                              cagcagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTgaa                                                           <  1:1086127/93‑1 (MQ=255)

                               agcagGAAGTGGGTACAGGTTACTTCGATAAGGTGACGACTATTATTCAGGGAGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTgaag                                                          >  1:2748204/1‑93 (MQ=255)

                               agcagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTgaag                                                          <  1:2484281/93‑1 (MQ=255)

                                gcagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTgaaga                                                         <  2:9407/93‑1 (MQ=255)

                                gcagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTgaaga                                                         <  2:3708995/93‑1 (MQ=255)

                                gcagGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTgaaga                                                         >  1:3076878/1‑93 (MQ=255)

                                  agGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAAt                                                       >  2:207419/1‑93 (MQ=255)

                                   ggAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATc                                                      >  2:4485163/1‑93 (MQ=255)

                                   ggAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATc                                                      <  2:3144963/93‑1 (MQ=255)

                                    gAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCg                                                     >  1:3823089/1‑93 (MQ=255)

                                     aaGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGc                                                    <  2:4265635/93‑1 (MQ=255)

                                       gTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAg                                                  <  1:100285/93‑1 (MQ=255)

                                        tGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGt                                                 <  2:3006666/93‑1 (MQ=255)

                                          ggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTc                                               <  2:3350508/93‑1 (MQ=255)

                                          ggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTc                                               <  2:830388/93‑1 (MQ=255)

                                          ggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTc                                               >  2:2527007/1‑93 (MQ=255)

                                          ggTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTc                                               >  1:3666107/1‑93 (MQ=255)

                                            gaCAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTa                                             <  2:1963498/92‑1 (MQ=255)

                                             aCAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTaa                                            <  1:428008/93‑1 (MQ=255)

                                              cAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAg                                           >  1:987/1‑93 (MQ=255)

                                                ggTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcc                                         >  2:4052539/1‑92 (MQ=255)

                                                ggTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGca                                         >  2:4556259/1‑93 (MQ=255)

                                                ggTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGACTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGca                                         >  1:4238474/1‑93 (MQ=255)

                                                 gTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaa                                        >  2:3095561/1‑93 (MQ=255)

                                                 gTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaa                                        >  2:1419772/1‑93 (MQ=255)

                                                  ttACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaac                                       >  1:481120/1‑93 (MQ=255)

                                                  ttACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaac                                       <  1:946582/93‑1 (MQ=255)

                                                    aCTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaa                                     >  2:3196880/1‑93 (MQ=255)

                                                    aCTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaa                                     >  2:2343923/1‑93 (MQ=255)

                                                    aCTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaa                                     <  2:1864069/93‑1 (MQ=255)

                                                    aCTTCGATAAAGGGACGACTATTATTCAGGGGGGCACGCCTTCAGTCACCGCGCTGACCGGCTGCACTGAAGAATCGCAGTTCTAAGcgacca                                     >  2:3653446/1‑91 (MQ=255)

                                                     cTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaac                                    <  2:3563598/93‑1 (MQ=255)

                                                     cTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaac                                    >  2:2082831/1‑93 (MQ=255)

                                                     cTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaac                                    >  1:1560730/1‑93 (MQ=255)

                                                     cTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaaa                                    >  2:1455016/1‑92 (MQ=255)

                                                      ttCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaaca                                   <  1:2238351/93‑1 (MQ=255)

                                                      ttCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaaca                                   >  2:1060697/1‑93 (MQ=255)

                                                       tCGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGcaacaacaa                                  >  2:2298995/1‑93 (MQ=255)

                                                        cGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCaacaacaac                                 <  1:4620339/93‑1 (MQ=255)

                                                        cGATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCaacaacaac                                 <  1:2334132/93‑1 (MQ=255)

                                                         gATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAAcc                                >  1:342617/1‑93 (MQ=255)

                                                         gATAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAAcc                                <  2:3967012/93‑1 (MQ=255)

                                                          aTAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCg                               >  2:3533366/1‑93 (MQ=255)

                                                          aTAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCg                               >  1:764654/1‑93 (MQ=255)

                                                           tAAAGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGt                              >  2:2763143/1‑93 (MQ=255)

                                                             aaGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTg                            >  1:1956276/1‑93 (MQ=255)

                                                             aaGTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAAAAACAGTTg                            >  1:3336032/1‑93 (MQ=255)

                                                               gTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGCCGCATCGCAGTTCTAAGCAACAACAACCGTAGCt                          >  1:527302/1‑93 (MQ=255)

                                                               gTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCt                          >  1:3615708/1‑93 (MQ=255)

                                                               gTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCt                          <  1:731237/93‑1 (MQ=255)

                                                               gTGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCt                          >  1:3960548/1‑93 (MQ=255)

                                                                tGACGACTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTg                         <  1:3868044/93‑1 (MQ=255)

                                                                  acgacTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGAc                       >  1:2596567/1‑93 (MQ=255)

                                                                  acgacTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGAc                       <  2:670272/93‑1 (MQ=255)

                                                                   cgacTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACt                      >  1:2882565/1‑93 (MQ=255)

                                                                    gacTATTATTCAGGACGGCACGCCTTCAGTCACCGCGCTGTCCTGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTg                     <  2:4114755/93‑1 (MQ=255)

                                                                      cTATTATTCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTa                   >  2:128770/1‑93 (MQ=255)

                                                                        attattCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAgg                 >  2:487070/1‑93 (MQ=255)

                                                                        attattCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAgg                 >  2:63501/1‑93 (MQ=255)

                                                                        attattCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAgg                 >  1:1730383/1‑93 (MQ=255)

                                                                        attattCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAgg                 >  1:4486581/1‑93 (MQ=255)

                                                                          tattCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGcc               >  2:2540662/1‑93 (MQ=255)

                                                                           attCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCg              >  2:2611931/1‑93 (MQ=255)

                                                                           attCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCg              <  1:915393/93‑1 (MQ=255)

                                                                           attCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCg              >  1:640416/1‑93 (MQ=255)

                                                                            ttCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCgg             >  1:1436327/1‑93 (MQ=255)

                                                                            ttCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCgg             >  1:4574918/1‑93 (MQ=255)

                                                                             tCAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGa            >  1:3602458/1‑93 (MQ=255)

                                                                              cAGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGAt           <  1:159890/93‑1 (MQ=255)

                                                                               aGGGCGGCACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATa          >  2:732664/1‑93 (MQ=255)

                                                                                  gcggcACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAgg       <  1:657053/93‑1 (MQ=255)

                                                                                   gggcACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGc      <  1:440787/92‑1 (MQ=255)

                                                                                      cACGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGtt   <  1:3192640/93‑1 (MQ=255)

                                                                                      cACGCCTTCAGTCACCGCGCTGACCGGCTACACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGtt   <  2:4604487/93‑1 (MQ=255)

                                                                                       aCGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTc  <  2:2806617/93‑1 (MQ=255)

                                                                                       aCGCCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTc  >  1:3784314/1‑93 (MQ=255)

GCGAAAGATGGCTATACCTTCGTATCTCACCAGCAGGAAGTGGGTACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTC  >  NC_000913/4215855‑4216034

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