BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,597,072	T→C	D73G (GAC→GGC)	b3460 ←	LivJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,597,072	0	T	C	99.1%	1404.4 / ‑3.7	451	D73G (GAC→GGC)	b3460	LivJ
Reads supporting (aligned to +/- strand): ref base T (2/2); new base C (241/205); total (244/207)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 457 total aligned reads displayed.

GAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGG  >  NC_000913/3596982‑3597161

 aagaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCat                                                                                        <  2:5989996/93‑1 (MQ=255)

  agaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCata                                                                                       <  1:4801201/93‑1 (MQ=255)

   gtaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCatat                                                                                      <  1:6145226/91‑1 (MQ=255)

   gaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCatat                                                                                      >  2:6642345/1‑93 (MQ=255)

    aaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAtt                                                                                     >  1:3259637/1‑93 (MQ=255)

    aaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAtt                                                                                     >  1:2543658/1‑93 (MQ=255)

     aCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAttt                                                                                    >  1:245029/1‑93 (MQ=255)

      cAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAtttt                                                                                   <  2:1531644/93‑1 (MQ=255)

       agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  >  1:585038/1‑93 (MQ=255)

        gagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 >  2:2411517/1‑93 (MQ=255)

        gagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 >  2:2997330/1‑93 (MQ=255)

        gagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 >  2:6446804/1‑93 (MQ=255)

         agGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTGCGCCACCGCACCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACg                                                                                <  1:186743/93‑1 (MQ=255)

         agGTGACCAATCACATATTTAATGCCGTCGCTAACGACTTTGTTCGCCACCGCAACCGCCCGCTTCGGGTCACAGGCATCGCCATATTTTACg                                                                                >  2:4866903/1‑93 (MQ=25)

          ggggACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACGGCAACCGCCAGTTTCGGGTCACAGGGATCGCATTATTTTACgt                                                                               >  2:1760568/4‑92 (MQ=255)

          ggTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACgt                                                                               >  1:5091401/1‑92 (MQ=255)

           gTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGAt                                                                              <  1:637245/93‑1 (MQ=255)

              aCCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTg                                                                           <  1:430641/93‑1 (MQ=255)

                cAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTgag                                                                         >  1:2670755/1‑91 (MQ=255)

                cAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCa                                                                         <  2:3391455/93‑1 (MQ=255)

                cAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCa                                                                         >  1:3872745/1‑93 (MQ=255)

                cAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCa                                                                         >  2:2636123/1‑93 (MQ=255)

                cAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCa                                                                         >  1:2237995/1‑93 (MQ=255)

                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  2:3581403/1‑93 (MQ=255)

                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGAATCGCCATATTTTACGATTTGCAg                                                                        >  1:4545771/1‑93 (MQ=255)

                     acaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttg                                                                    >  2:4766348/1‑93 (MQ=255)

                      caTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgt                                                                   >  1:120772/1‑93 (MQ=255)

                       atatTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgtt                                                                  <  2:5252093/93‑1 (MQ=255)

                        tatTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCGGTttgttg                                                                 >  2:7624039/1‑93 (MQ=255)

                        tatTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgttg                                                                 >  2:896030/1‑93 (MQ=255)

                        tatTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgttg                                                                 >  2:5745236/1‑93 (MQ=255)

                        tatTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgttg                                                                 >  2:6126547/1‑93 (MQ=255)

                         atatAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGc                                                                <  2:4565662/90‑1 (MQ=255)

                         atTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGc                                                                >  1:5993813/1‑93 (MQ=255)

                         atTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGc                                                                <  2:6834457/93‑1 (MQ=255)

                          tttATTGCCGTCGTTCACGTCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGCTTTGCAGTTTGTTGcc                                                               <  2:7027449/93‑1 (MQ=255)

                          tttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGct                                                               >  1:4873280/1‑92 (MQ=255)

                          tttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGcc                                                               <  1:5355598/93‑1 (MQ=255)

                           ttACTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCt                                                              >  2:5542545/1‑93 (MQ=255)

                           ttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCt                                                              <  1:1662237/93‑1 (MQ=255)

                           ttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCt                                                              <  1:7096484/93‑1 (MQ=255)

                            tAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCtt                                                             <  1:2301775/93‑1 (MQ=255)

                              aTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTa                                                           >  2:1661223/1‑93 (MQ=255)

                              aTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTa                                                           >  1:771241/1‑93 (MQ=255)

                                  cgtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgc                                                       >  2:71800/1‑93 (MQ=255)

                                  cgtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgc                                                       <  1:854244/93‑1 (MQ=255)

                                  cgtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgc                                                       >  2:2807205/1‑93 (MQ=255)

                                  cgtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgc                                                       <  1:5310665/93‑1 (MQ=255)

                                  cgccgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGTCATCGCCATATTTTACGATCTGCAGTTTGTTGCCTTTAATgc                                                       <  1:1330000/90‑1 (MQ=255)

                                   gtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgcc                                                      >  2:4054001/1‑93 (MQ=255)

                                    tcgtTCGCGCCTTTGTTCGCCCCCGCCGCCGCCTGTTTCGGGTCCCGGGCCTCGCCATTTTTTCCGCTTTGCCGTTTGTTGCCTTTCCTgccg                                                     <  2:6037226/93‑1 (MQ=255)

                                    tcgtTAACGACTTTGTTCGCCACCGCACCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCCGTTTGTTGCCTTTAATgacg                                                     >  1:3789461/1‑93 (MQ=255)

                                      gtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgccgcc                                                   >  2:4464706/1‑93 (MQ=255)

                                       ttAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCt                                                  >  2:3878988/1‑93 (MQ=255)

                                       ttAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCt                                                  >  2:1386688/1‑93 (MQ=255)

                                        tAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCtt                                                 <  2:1142801/93‑1 (MQ=255)

                                        tAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGATTGTTGACTTTAATGCCGCCtt                                                 >  2:5161063/1‑93 (MQ=255)

                                         aaCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCttt                                                <  2:7756469/93‑1 (MQ=255)

                                          aCGCCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTa                                               <  2:2237378/93‑1 (MQ=255)

                                          aCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTa                                               <  1:2107957/93‑1 (MQ=255)

                                            gACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGc                                             >  1:2899815/1‑93 (MQ=255)

                                             aCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCg                                            >  1:5246657/1‑93 (MQ=255)

                                             aCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCg                                            >  2:7406969/1‑93 (MQ=255)

                                             aCTTTGTTCGACACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGct                                            >  2:2070287/1‑92 (MQ=255)

                                              cTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGt                                           >  1:3672897/1‑93 (MQ=255)

                                               tttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGtt                                          >  2:1311274/1‑93 (MQ=255)

                                               tttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGtt                                          <  1:310364/93‑1 (MQ=255)

                                                ttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTg                                         <  2:4280356/93‑1 (MQ=255)

                                                   ttCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGata                                      >  1:3398126/1‑93 (MQ=255)

                                                     cGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATc                                    <  2:7298689/93‑1 (MQ=255)

                                                      gCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATcc                                   >  1:6638949/1‑93 (MQ=255)

                                                        cACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGc                                 >  1:5713250/1‑93 (MQ=255)

                                                        cACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGc                                 >  2:3973810/1‑93 (MQ=255)

                                                          ccGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCaa                               <  1:5542546/93‑1 (MQ=255)

                                                             cAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCg                            >  2:1924098/1‑93 (MQ=255)

                                                              aaCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGc                           >  1:6226653/1‑93 (MQ=255)

                                                              aaCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCGTTAGCGTTGATATCCGCAACCGc                           >  2:4344492/1‑93 (MQ=255)

                                                               aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGcc                          >  1:3537669/1‑93 (MQ=255)

                                                               aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGcc                          >  1:1969869/1‑93 (MQ=255)

                                                                ccGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCt                         <  2:3745431/93‑1 (MQ=255)

                                                                 cGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTg                        >  1:6873650/1‑93 (MQ=255)

                                                                  gCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGc                       >  1:7347258/1‑93 (MQ=255)

                                                                    cTGTTTCGGGTGACAGGCATCGACATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCTGCCTTTAGCCTTGATATCCGCGACCGCCTGctc                     >  2:5405001/1‑93 (MQ=255)

                                                                    cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     >  1:6544963/1‑93 (MQ=255)

                                                                    cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     <  2:4671079/93‑1 (MQ=255)

                                                                     tGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctct                    >  1:4677680/1‑93 (MQ=255)

                                                                     tGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctct                    <  2:7838754/93‑1 (MQ=255)

                                                                     tGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctct                    <  1:7879588/93‑1 (MQ=255)

                                                                          cGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCATCCGCCTGCTCTGCGcc               <  1:17693/93‑1 (MQ=255)

                                                                               cacaGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTaa          <  1:4951379/93‑1 (MQ=255)

                                                                                    gCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTcc     <  2:1012084/93‑1 (MQ=255)

                                                                                    gCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTcc     <  2:4978036/93‑1 (MQ=255)

                                                                                     cATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCt    >  1:3985694/1‑93 (MQ=255)

                                                                                     cATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCt    >  1:5757813/1‑93 (MQ=255)

                                                                                     cATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCt    <  2:3108065/93‑1 (MQ=255)

                                                                                       tCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  >  1:7428740/1‑93 (MQ=255)

                                                                                       tCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  <  2:2734165/93‑1 (MQ=255)

                                                                                       tCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  >  2:3363923/1‑93 (MQ=255)

                                                                                       tCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  >  1:1000657/1‑93 (MQ=255)

                                                                                       tCGCCATATTTTACGATTTGCAGCTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  >  2:6315642/1‑93 (MQ=255)

GAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGG  >  NC_000913/3596982‑3597161

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