BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,562,600	T→C	T274A (ACC→GCC)	yjiM ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,562,600	0	T	C	99.4%	948.2 / ‑4.1	309	T274A (ACC→GCC)	yjiM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (176/131); total (177/132)
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.96e-01

Only 100 of 330 total aligned reads displayed.

AGGAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGTGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGA  >  NC_000913/4562510‑4562689

 ggAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGgcgcc                                                                                         >  1:5694049/1‑92 (MQ=255)

 ggAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCa                                                                                         <  2:5829294/93‑1 (MQ=255)

  gAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAg                                                                                        >  2:4702654/1‑93 (MQ=255)

      agcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGtttt                                                                                    <  2:733302/93‑1 (MQ=255)

      agcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGtttt                                                                                    >  2:6707793/1‑93 (MQ=255)

       gcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTc                                                                                   >  2:7855188/1‑93 (MQ=255)

           aTCGCAAAATATTTATCCGCCAGCGCGTCGTAAGCATCGCCCGTTTCTGCCGCGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATaa                                                                               <  2:1486170/93‑1 (MQ=255)

            tCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAAc                                                                              <  2:6307393/93‑1 (MQ=255)

              gCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAAcag                                                                            >  1:1275727/1‑91 (MQ=255)

              atcAGATATTGATCCGCCAGCGCGACGTAGACATCGCCCGTTTCTGCCACGCATGGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCg                                                                            <  2:2974063/91‑1 (MQ=255)

               ccAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGa                                                                           >  2:7930252/1‑93 (MQ=255)

               ccAGATATTTATACGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGa                                                                           <  1:1003205/93‑1 (MQ=255)

                cAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGAc                                                                          >  2:4427112/1‑93 (MQ=255)

                cAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGAc                                                                          >  2:3551748/1‑93 (MQ=255)

                 aGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACa                                                                         >  1:3719493/1‑93 (MQ=255)

                 aGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGCTTTCATAACCGACa                                                                         >  1:7345948/1‑93 (MQ=255)

                  gATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACaa                                                                        >  1:4349297/1‑93 (MQ=255)

                   atatTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAAc                                                                       <  1:3580832/93‑1 (MQ=255)

                   atatTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAAc                                                                       <  1:1489695/93‑1 (MQ=255)

                    tatTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAAcc                                                                      >  1:5433987/1‑93 (MQ=255)

                     atTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAAccc                                                                     >  2:2436186/1‑93 (MQ=255)

                      tttATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATCACCGACAACCCa                                                                    <  2:3714215/93‑1 (MQ=255)

                      tttATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCa                                                                    <  1:1178949/93‑1 (MQ=255)

                      tttATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCa                                                                    >  1:3136492/1‑93 (MQ=255)

                            cGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑gccgcc                                                             <  2:5406833/93‑1 (MQ=255)

                               cAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCAtt                                                          <  1:7611431/93‑1 (MQ=255)

                                aggacgtggtAGGCATCCACCGTTTCCGCCCCGCATTGCTCGGTGCCTTTCGCCCGGGCGCAGTTTTCATAACCGAAAGCCCA‑GCCGCCATtc                                                         <  2:5978841/89‑1 (MQ=255)

                                aGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATtc                                                         >  2:1637335/1‑93 (MQ=255)

                                    cgtcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCtt                                                     <  2:5047142/93‑1 (MQ=255)

                                     gtcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTc                                                    <  1:5189008/93‑1 (MQ=255)

                                     gtcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTc                                                    <  1:1644045/93‑1 (MQ=255)

                                     gtcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTc                                                    <  1:1812826/93‑1 (MQ=255)

                                     gtcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTc                                                    <  1:6933893/93‑1 (MQ=255)

                                     gtcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTc                                                    >  1:5475428/1‑93 (MQ=255)

                                      tcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCa                                                   >  1:68107/1‑93 (MQ=255)

                                      tcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCa                                                   <  1:2033329/93‑1 (MQ=255)

                                      tcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCa                                                   >  2:4805563/1‑93 (MQ=255)

                                      tcgtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCa                                                   >  2:6386418/1‑93 (MQ=255)

                                        gtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAAt                                                 <  1:5832269/93‑1 (MQ=255)

                                        gtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAAt                                                 >  2:2719742/1‑93 (MQ=255)

                                         tAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATc                                                >  2:3986451/1‑93 (MQ=255)

                                         tAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATc                                                >  2:5637433/1‑93 (MQ=255)

                                         tAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATc                                                >  2:2882114/1‑93 (MQ=255)

                                         tAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATc                                                >  1:3229234/1‑93 (MQ=255)

                                         tAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATc                                                >  1:4653561/1‑93 (MQ=255)

                                          aGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATcg                                               >  2:4141308/1‑93 (MQ=255)

                                           gACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATcgc                                              <  1:2522352/93‑1 (MQ=255)

                                           gACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATcgc                                              >  2:2599918/1‑93 (MQ=255)

                                           gACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATcgc                                              >  2:6857312/1‑93 (MQ=255)

                                           gACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCGTCAATcgc                                              >  1:62039/1‑93 (MQ=255)

                                             cATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATACCCGACAACCCA‑GCCGCCATTCTCTTCAATcgcgc                                            <  1:4279104/93‑1 (MQ=255)

                                             cATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATcgcgc                                            <  2:2040132/93‑1 (MQ=255)

                                                cGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGca                                         >  1:2673325/1‑93 (MQ=255)

                                                cGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGca                                         >  1:3050884/1‑93 (MQ=255)

                                                  cccGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGcacc                                       >  2:5760499/1‑93 (MQ=255)

                                                  cccGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGcacc                                       >  1:3811366/1‑93 (MQ=255)

                                                  cccGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGcacc                                       >  1:665526/1‑93 (MQ=255)

                                                   ccGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGcacca                                      >  2:7757682/1‑93 (MQ=255)

                                                   ccGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGcacca                                      >  1:4960079/1‑93 (MQ=255)

                                                     gTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACt                                    >  1:2093586/1‑93 (MQ=255)

                                                      tttCTTCCACGCATTGCTCGGTCGCTTTCGCCCCGGCTCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACtt                                   <  1:4074691/93‑1 (MQ=255)

                                                      tttCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACtt                                   >  2:710782/1‑93 (MQ=255)

                                                       ttCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACttt                                  <  1:3906224/93‑1 (MQ=255)

                                                        tCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACtttt                                 <  2:5127413/93‑1 (MQ=255)

                                                         ttGCCCCGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACttttt                                <  2:4608479/92‑1 (MQ=255)

                                                         cTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACttttt                                >  1:7296264/1‑93 (MQ=255)

                                                          tGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTc                               >  1:3034100/1‑93 (MQ=255)

                                                           gCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTct                              <  2:1290562/93‑1 (MQ=255)

                                                           gCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTct                              <  1:2746609/93‑1 (MQ=255)

                                                           gCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTct                              >  1:226196/1‑93 (MQ=255)

                                                             cACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGGTTTTATAAACCACAAACCACCCCCCCCTCCTCTTCATCTCCCCCCACCACtttttttc                             >  2:6642887/1‑90 (MQ=255)

                                                             cACGCATTGCTCGGTCGCTTTCGACCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgc                            <  2:4353403/93‑1 (MQ=255)

                                                              aCGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgct                           >  2:6089940/1‑93 (MQ=255)

                                                              aCGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgct                           >  2:5903476/1‑93 (MQ=255)

                                                                 cATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTATCATACCCGCCAACCCA‑GCCGCCGTTCTCTTCGATCGCGCGCACCACTTTTTCTGCTgcg                        <  2:7222646/93‑1 (MQ=255)

                                                                   ttGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGcgcc                      >  2:52381/1‑93 (MQ=255)

                                                                   ttGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGcgcc                      >  1:7026130/1‑93 (MQ=255)

                                                                      cTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCgccgcc                   >  1:6736166/1‑93 (MQ=255)

                                                                       tCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCa                  >  2:1874178/1‑93 (MQ=255)

                                                                       tCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCa                  >  1:781126/1‑93 (MQ=255)

                                                                          gTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATc               >  2:3133125/1‑93 (MQ=255)

                                                                          gTCGCTTTCGCCCCGGCGAAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATc               <  1:4206504/93‑1 (MQ=255)

                                                                           tCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCg              >  1:4561289/1‑93 (MQ=255)

                                                                           tCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCg              >  2:6906487/1‑93 (MQ=255)

                                                                           tCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCg              <  2:1287222/93‑1 (MQ=255)

                                                                             gCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCggg            >  2:227674/1‑93 (MQ=255)

                                                                             gCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCggg            <  1:1469706/93‑1 (MQ=255)

                                                                              cTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCCGGc           >  1:4438290/1‑93 (MQ=255)

                                                                              cTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCC‑CCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCa          <  2:3367384/93‑1 (MQ=255)

                                                                               tttCGCCCCGGCGCAGTTTTCATACCCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCa          <  1:121063/93‑1 (MQ=255)

                                                                               tttCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCa          <  1:7651395/93‑1 (MQ=255)

                                                                                ttCGCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAg         >  2:4623421/1‑93 (MQ=255)

                                                                                   gCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCg      >  1:5350068/1‑93 (MQ=255)

                                                                                   gCCCCGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCg      >  1:6369893/1‑93 (MQ=255)

                                                                                    ccccGGTGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCgg     >  1:6516849/1‑93 (MQ=255)

                                                                                    ccccGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCgg     >  1:869942/1‑93 (MQ=255)

                                                                                      ccGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTg   >  1:6333122/1‑93 (MQ=255)

                                                                                      ccGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTg   >  2:4077391/1‑93 (MQ=255)

                                                                                      ccGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTg   >  2:1003399/1‑93 (MQ=255)

                                                                                       cGGCGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGa  >  2:551136/1‑93 (MQ=255)

AGGAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGTGCAGTTTTCATAACCGACAACCCA‑GCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGA  >  NC_000913/4562510‑4562689

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: `‑`