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%	1062.3 / ‑5.0	341	T274A (ACC→GCC)	yjiM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (183/155); minor base G (2/0); total (186/155)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 359 total aligned reads displayed.

AGGAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGTGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGA  >  NC_000913/4562510‑4562689

 ggAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTTCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCa                                                                                        >  1:4788225/1‑93 (MQ=255)

 ggAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCa                                                                                        >  2:4174586/1‑93 (MQ=255)

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

       gcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCGGTTTTc                                                                                  >  1:6368349/1‑93 (MQ=255)

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

       gcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTc                                                                                  <  2:7437107/93‑1 (MQ=255)

       gcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTc                                                                                  <  2:605571/93‑1 (MQ=255)

       gcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTc                                                                                  <  2:2693639/93‑1 (MQ=255)

       gcCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTc                                                                                  <  2:1584037/93‑1 (MQ=255)

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

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

        ccAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCa                                                                                 <  1:6384538/93‑1 (MQ=255)

         cAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCAt                                                                                >  2:3932988/1‑93 (MQ=255)

          aaTCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCGCGCAGTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATa                                                                               <  2:6585645/93‑1 (MQ=255)

          aaTCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATa                                                                               <  2:5385029/93‑1 (MQ=255)

          aaTCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATa                                                                               <  1:1313227/93‑1 (MQ=255)

          aaTCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATa                                                                               >  1:6993406/1‑93 (MQ=255)

           aTCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATaa                                                                              <  2:1120564/93‑1 (MQ=255)

           aTCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATaa                                                                              >  2:2143129/1‑93 (MQ=255)

           aTCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGACCAGTTTTCATaa                                                                              >  1:4664896/1‑93 (MQ=255)

            tCGCAAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAAc                                                                             <  2:7395446/93‑1 (MQ=255)

                cAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGAc                                                                         <  2:514888/93‑1 (MQ=255)

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

                 agagATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACa                                                                        <  2:7628877/93‑1 (MQ=255)

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

                  gATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGTTCCGGCGCAGTTTTCATAACCGACaa                                                                       >  1:9941243/1‑93 (MQ=255)

                  aatatTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACaa                                                                       <  2:3227112/92‑1 (MQ=255)

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

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

                        tATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgc                                                                 >  1:5197549/1‑93 (MQ=255)

                        tATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgc                                                                 >  2:4688964/1‑93 (MQ=255)

                         aTCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgcc                                                                <  2:5700549/93‑1 (MQ=255)

                          tCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgccg                                                               >  2:7819978/1‑93 (MQ=255)

                           ccGCCAGCGCGTCGTAGATATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAgccgc                                                              <  1:5737565/93‑1 (MQ=255)

                              ccAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCAt                                                           <  2:5792876/93‑1 (MQ=255)

                              ccAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCAt                                                           <  2:8455974/93‑1 (MQ=255)

                              ccAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCAt                                                           >  1:8303031/1‑93 (MQ=255)

                              ccAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCAt                                                           >  2:642571/1‑93 (MQ=255)

                                aGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATtc                                                         >  1:7341945/1‑93 (MQ=255)

                                aGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATtc                                                         <  1:10044484/93‑1 (MQ=255)

                                 gcgcGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATtct                                                        >  1:7939938/1‑93 (MQ=255)

                                 gcgcGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATtct                                                        <  1:3768519/93‑1 (MQ=255)

                                   gcgTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTctct                                                      <  1:2909620/93‑1 (MQ=255)

                                   gcgTCGTAGACATCGCCCGTTTCTCCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTctct                                                      >  1:5529583/1‑93 (MQ=255)

                                      ccctAGACATCGCCCGTTTCTCCCACGCATTGCCCGGTCGCTTTCGCCCCGGCGCAGTTCCCATAACCGACAACCCAGCCGCCATTCTCTTCa                                                   <  1:2245738/90‑1 (MQ=255)

                                        gtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAAt                                                 <  1:7755597/93‑1 (MQ=255)

                                        gtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAAt                                                 <  2:4990327/93‑1 (MQ=255)

                                        gtAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAAt                                                 <  2:8304517/93‑1 (MQ=255)

                                         tAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATc                                                >  2:7798182/1‑93 (MQ=255)

                                         tAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATc                                                <  2:1766560/93‑1 (MQ=255)

                                         tAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATc                                                >  1:276003/1‑93 (MQ=255)

                                         tAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATc                                                >  1:6859903/1‑93 (MQ=255)

                                          aGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATcg                                               <  1:3113069/93‑1 (MQ=255)

                                          aGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATcg                                               >  2:4495183/1‑93 (MQ=255)

                                            aCATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATcgcg                                             >  2:9620410/1‑93 (MQ=255)

                                            aCATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATcgcg                                             <  1:6780080/93‑1 (MQ=255)

                                            aCATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATcgcg                                             <  1:2177612/93‑1 (MQ=255)

                                               tCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCgcgcgc                                          >  2:2460589/1‑93 (MQ=255)

                                               tCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCgcgcgc                                          >  1:4026279/1‑93 (MQ=255)

                                                  cccGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGGCGCCATTCTCTTCAATCGCGCGcacc                                       >  2:8745264/1‑93 (MQ=255)

                                                    cGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCACCGCGCAGTTTTCATAACCGACAACACAGCCGCCATTCTCTTCAATCTCGCGcaccac                                     >  1:9315021/1‑93 (MQ=255)

                                                     gTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACt                                    <  1:9613157/93‑1 (MQ=255)

                                                      tttCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACtt                                   <  2:9270225/93‑1 (MQ=255)

                                                       ttCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACttt                                  <  1:7891142/93‑1 (MQ=255)

                                                        tCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACACCCCAGCCGCCATTCTCTTCAATCGCGCGCACCACtttt                                 <  1:6212569/93‑1 (MQ=255)

                                                         cTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACttttt                                >  1:8710544/1‑93 (MQ=255)

                                                          tGCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTc                               >  2:3443947/1‑93 (MQ=255)

                                                           gCCACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCACGCGCACCACTTTTTct                              >  2:3064661/1‑93 (MQ=255)

                                                             cACGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgc                            >  2:8079514/1‑93 (MQ=255)

                                                              aCGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgct                           >  2:8218632/1‑93 (MQ=255)

                                                               cGCATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTctgctg                          <  1:945866/93‑1 (MQ=255)

                                                                 cATTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTgcg                        >  2:6863004/1‑93 (MQ=255)

                                                                  aTTGCTCGGTCGCTTTCGCCCCGGGGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGACCCCCTTTTTCttctgccc                       >  1:8404741/1‑91 (MQ=255)

                                                                  aTTGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTgcgc                       >  2:8053147/1‑93 (MQ=255)

                                                                   ttGCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGcgcc                      >  1:6570569/1‑93 (MQ=255)

                                                                     gCTCGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGcgccgc                    <  2:2922143/93‑1 (MQ=255)

                                                                      ctcgctcgCTATCGCCCCGTCGTAGTTTTCAAACCCGACACCCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCtccgcc                   <  2:49710/93‑1 (MQ=255)

                                                                        cGGTCTCTTTCGCCCCGGCGCAGTTTTCATAACCGACACCCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCaa                 <  1:2601428/93‑1 (MQ=255)

                                                                        cGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCaa                 <  1:5082880/93‑1 (MQ=255)

                                                                        cGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGAGCACCACTTTTTCTGCTGCGCCGCCaa                 >  1:3694250/1‑93 (MQ=255)

                                                                        cGGTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCCCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCaa                 <  2:7165575/93‑1 (MQ=255)

                                                                         ggTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAAt                <  1:3465762/93‑1 (MQ=255)

                                                                          gTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATc               <  2:973043/93‑1 (MQ=255)

                                                                          gTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATc               >  2:9047122/1‑93 (MQ=255)

                                                                          gTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATc               >  1:9658909/1‑93 (MQ=255)

                                                                          gTCGCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATc               >  1:1977133/1‑93 (MQ=255)

                                                                             gCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCggg            <  2:2012450/93‑1 (MQ=255)

                                                                             gCTTTCGCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCggg            <  2:8634965/93‑1 (MQ=255)

                                                                                   gCCCCGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCg      <  2:4860434/93‑1 (MQ=255)

                                                                                    ccccGGCGCAGTTTTCATAACCGACACCCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCgg     <  2:1481388/93‑1 (MQ=255)

                                                                                    ccccGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCgg     >  1:6077454/1‑93 (MQ=255)

                                                                                    ccccGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCAAATCGGGAAGCCgg     >  1:3037087/1‑93 (MQ=255)

                                                                                     cccGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGt    >  1:4401133/1‑93 (MQ=255)

                                                                                     cccGGCGCAGTTTTAATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGt    <  2:920851/93‑1 (MQ=255)

                                                                                     cccGGCGCAGTTTGCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGt    <  1:62703/93‑1 (MQ=255)

                                                                                      ccGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTg   >  2:1038628/1‑93 (MQ=255)

                                                                                      ccGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTg   >  1:9153454/1‑93 (MQ=255)

                                                                                       cGGCGCAGTTTTCATACCCAACACCCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTATTCTGCTGCGCCGCCAATCGGACAGCCGGTGa  <  1:8926914/93‑1 (MQ=255)

                                                                                       cGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGa  >  1:2879037/1‑93 (MQ=255)

                                                                                       cGGCGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGa  >  1:1143363/1‑93 (MQ=255)

AGGAGCAGCCAATCGCCAGATATTTATCCGCCAGCGCGTCGTAGACATCGCCCGTTTCTGCCACGCATTGCTCGGTCGCTTTCGCCCCGGTGCAGTTTTCATAACCGACAACCCAGCCGCCATTCTCTTCAATCGCGCGCACCACTTTTTCTGCTGCGCCGCCAATCGGGCAGCCGGTGA  >  NC_000913/4562510‑4562689

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 17 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`