BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,917,420	(C)_5→6	coding (846/864 nt)	atpG ←	ATP synthase F1 complex subunit gamma

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,917,415	1	.	C	99.6%	904.2 / ‑3.3	228	coding (851/864 nt)	atpG	ATP synthase F1 complex subunit gamma
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (140/87); total (141/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

Only 100 of 265 total aligned reads displayed.

CCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGG‑CCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTGCCGCCATTG  >  NC_000913/3917315‑3917510

ccTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGcggcgg                                                                                                  >  1:5157528/1‑101 (MQ=255)

    aGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGG‑cccc                                                                                             >  2:5157532/1‑101 (MQ=255)

      ggAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGG‑cccccc                                                                                           >  1:5157534/1‑100 (MQ=255)

        aaTTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGG‑ccccccga                                                                                         >  1:5157538/1‑98 (MQ=255)

         aTTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGG‑ccccccgag                                                                                        >  2:5157540/1‑97 (MQ=255)

           tCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGcggcgg                                                                                                  >  2:5157544/1‑90 (MQ=255)

           tCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGAc                                                                                       >  1:5157543/1‑101 (MQ=255)

            cGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACg                                                                                      <  2:5157002/101‑1 (MQ=255)

             gACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGa                                                                                     >  1:5157548/1‑101 (MQ=255)

              aCGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGAt                                                                                    >  2:5157549/1‑101 (MQ=255)

                  cAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCg                                                                                <  2:5157168/101‑1 (MQ=255)

                  cAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCg                                                                                >  2:5157554/1‑101 (MQ=255)

                  cAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCg                                                                                <  1:5157252/101‑1 (MQ=255)

                   aaCTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACg                                                                                      >  1:5157556/1‑94 (MQ=255)

                     cTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTg                                                                             >  2:5157557/1‑101 (MQ=255)

                      tACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGa                                                                            <  1:5157393/101‑1 (MQ=255)

                          gcgcCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTc                                                                        >  2:5157561/1‑101 (MQ=255)

                            gcCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCt                                                                      >  2:5157565/1‑101 (MQ=255)

                             ccGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTg                                                                     >  1:5157566/1‑101 (MQ=255)

                              cGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGa                                                                    >  2:5157568/1‑101 (MQ=255)

                                 ttACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTa                                                                 >  1:5157571/1‑101 (MQ=255)

                                  tACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTg                                                                     >  1:5157575/1‑96 (MQ=255)

                                   aCCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAAt                                                               >  2:5157576/1‑101 (MQ=255)

                                    ccTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATg                                                              <  1:5157580/101‑1 (MQ=255)

                                    ccTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATg                                                              <  1:5157165/101‑1 (MQ=255)

                                     cTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTaa                                                                <  1:5157216/98‑1 (MQ=255)

                                      tGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCt                                                            >  2:5157583/1‑101 (MQ=255)

                                       ggACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTg                                                           <  1:5157406/101‑1 (MQ=255)

                                       ggACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTg                                                           >  1:5157584/1‑101 (MQ=255)

                                       ggACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTg                                                           <  1:5156946/101‑1 (MQ=255)

                                          cAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAg                                                                   >  2:5157588/1‑90 (MQ=255)

                                          cAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGcc                                                        <  2:5157334/101‑1 (MQ=255)

                                           aaTCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGAc                                                                                       <  2:5157128/69‑1 (MQ=255)

                                           aaTCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCt                                                       >  2:5157593/1‑101 (MQ=255)

                                           aaTCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCt                                                       >  1:5157591/1‑101 (MQ=255)

                                           aaTCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGG‑cc                                                                                               >  1:5157596/1‑60 (MQ=255)

                                            aTCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTg                                                      <  1:5157102/101‑1 (MQ=255)

                                            aTCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTg                                                      <  2:5157281/101‑1 (MQ=255)

                                                ttCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGa                                                  >  1:5157600/1‑101 (MQ=255)

                                                 tCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAg                                                 >  2:5157602/1‑101 (MQ=255)

                                                   cAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCt                                               >  1:5157605/1‑101 (MQ=255)

                                                    aGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCtt                                              >  2:5157607/1‑101 (MQ=255)

                                                     gTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACTAGCttt                                             <  2:5157275/101‑1 (MQ=255)

                                                     gTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCttt                                             <  2:5157189/101‑1 (MQ=255)

                                                     gTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCctt                                             <  1:5157260/101‑1 (MQ=255)

                                                      tAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCt                                               <  2:5157360/98‑1 (MQ=255)

                                                      tAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTttg                                            >  1:5157611/1‑101 (MQ=255)

                                                        gCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTttgtt                                          >  1:5157614/1‑101 (MQ=255)

                                                         ccATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTttgttg                                         >  2:5157616/1‑101 (MQ=255)

                                                         ccATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTttgttg                                         >  1:5157619/1‑101 (MQ=255)

                                                          cATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTtgttgt                                        >  2:5157621/1‑101 (MQ=255)

                                                          cATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTtgttgt                                        <  2:5157152/101‑1 (MQ=255)

                                                           aTCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGta                                       >  1:5157628/1‑101 (MQ=255)

                                                           aTCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGta                                       >  2:5157627/1‑101 (MQ=255)

                                                           aTCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGta                                       >  1:5157625/1‑101 (MQ=255)

                                                             cTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTg                                                                             <  1:5157068/61‑1 (MQ=255)

                                                             cTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGtata                                     >  2:5157632/1‑101 (MQ=255)

                                                               tAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATAcc                                   >  1:5157634/1‑101 (MQ=255)

                                                                aaaTCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCa                                  >  2:5157636/1‑101 (MQ=255)

                                                                aaaTCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCa                                  <  1:5157272/101‑1 (MQ=255)

                                                                 aaTCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCaa                                 <  2:5157077/101‑1 (MQ=255)

                                                                  aTCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAAc                                >  1:5157639/1‑101 (MQ=255)

                                                                   tCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACt                               <  2:5157173/101‑1 (MQ=255)

                                                                   tCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACt                               <  1:5157118/101‑1 (MQ=255)

                                                                    ccTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTg                              >  2:5157645/1‑101 (MQ=255)

                                                                    ccTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTg                              >  1:5157643/1‑101 (MQ=255)

                                                                      tctACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCgg                                                                               >  1:5157649/1‑50 (MQ=255)

                                                                      tctACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCa                            <  1:5157347/101‑1 (MQ=255)

                                                                       ctACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAg                           >  2:5157651/1‑101 (MQ=255)

                                                                       ctACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAg                           <  1:5157234/101‑1 (MQ=255)

                                                                        tACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGc                          >  1:5157652/1‑101 (MQ=255)

                                                                        tACGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGc                          >  2:5157654/1‑101 (MQ=255)

                                                                         aCGAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGct                         <  2:5157253/101‑1 (MQ=255)

                                                                           gAAATAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGctct                       >  2:5157658/1‑101 (MQ=255)

                                                                             aaTAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCttt                     <  1:5157340/101‑1 (MQ=255)

                                                                              aTAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTa                    >  1:5157663/1‑101 (MQ=255)

                                                                              aTAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTa                    <  2:5157099/101‑1 (MQ=255)

                                                                              aTAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTa                    <  1:5157365/101‑1 (MQ=255)

                                                                               tAACCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGa                                                                            >  2:5157666/1‑44 (MQ=255)

                                                                                aaCCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAAt                  >  1:5157667/1‑101 (MQ=255)

                                                                                aaCCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAAt                  <  1:5157204/101‑1 (MQ=255)

                                                                                 aCCTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGc                                                >  1:5157670/1‑70 (MQ=255)

                                                                                  ccTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCa                >  2:5157671/1‑101 (MQ=255)

                                                                                  ccTGTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCa                <  2:5157214/101‑1 (MQ=255)

                                                                                     gTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGc             <  2:5157225/101‑1 (MQ=255)

                                                                                     gTTTAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGc             >  1:5157676/1‑101 (MQ=255)

                                                                                      tttAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTg                              >  2:5157677/1‑83 (MQ=255)

                                                                                      tttAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCt            <  2:5157268/101‑1 (MQ=255)

                                                                                        tAAACCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTgc          >  2:5157681/1‑101 (MQ=255)

                                                                                         aaaCCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTgcc         >  1:5157684/1‑101 (MQ=255)

                                                                                          aagcgcgGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTgccg        >  1:5157686/4‑101 (MQ=255)

                                                                                          aaCCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTgccg        <  2:5157302/101‑1 (MQ=255)

                                                                                          aaCCGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTgccg        >  2:5157687/1‑101 (MQ=255)

                                                                                            ccGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTgccgcc      <  2:5157263/101‑1 (MQ=255)

                                                                                            ccGCGGCGGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTgccgcc      <  1:5157299/101‑1 (MQ=255)

                                                                                              gcggcgGCCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTGCCGCCAt    >  2:5157691/1‑101 (MQ=255)

                                                                                               cggcggcCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTGCCGCCAtt   <  1:5157324/101‑1 (MQ=255)

                                                                                               cggcggcCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTGCCGCCAtt   <  1:5157296/101‑1 (MQ=255)

                                                                                                ggcggcCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTGCCGCCATTg  >  1:5157698/1‑101 (MQ=255)

                                                                                                ggcggcCCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTGCCGCCATTg  >  1:5157694/1‑101 (MQ=255)

CCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGG‑CCCCCGAGACGATCTCGGTGAGTTCCTGAGTAATGCTGGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTGCCGCCATTG  >  NC_000913/3917315‑3917510

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