BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,961,588	C→G	R86T (AGG→ACG)	gppA ←	guanosine pentaphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,961,588	0	C	G	100.0%	2494.5 / NA	724	R86T (AGG→ACG)	gppA	guanosine pentaphosphatase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (391/333); total (391/333)

Only 100 of 738 total aligned reads displayed.

TCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATT  >  NC_000913/3961498‑3961677

 ccTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCg                                                                                        >  1:740041/1‑93 (MQ=255)

 ccTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCg                                                                                        >  2:2476252/1‑93 (MQ=255)

  ctctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGc                                                                                       >  1:4658548/1‑93 (MQ=255)

   tctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGcc                                                                                      >  2:7199811/1‑93 (MQ=255)

   tctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGcc                                                                                      >  2:1936909/1‑93 (MQ=255)

   tctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGcc                                                                                      >  2:3116798/1‑93 (MQ=255)

    ctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCg                                                                                     <  1:5711967/93‑1 (MQ=255)

    ctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCg                                                                                     >  1:4437967/1‑93 (MQ=255)

      tCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTa                                                                                   >  1:6799316/1‑93 (MQ=255)

        cccGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGc                                                                                 <  1:4829253/92‑1 (MQ=255)

        aCCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGc                                                                                 <  1:4215660/93‑1 (MQ=255)

           gCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCACTAAAATCACCCGCATTGACGGCAAGGCGTAACATCGCCGTAGCGaa                                                                              >  1:6617692/1‑92 (MQ=255)

           gCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGAc                                                                              >  2:8331932/1‑93 (MQ=255)

           gCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGAc                                                                              <  2:7804124/93‑1 (MQ=255)

            cTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACa                                                                             >  1:2057026/1‑93 (MQ=255)

             tGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                                            >  1:628937/1‑93 (MQ=255)

              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           >  2:4231597/1‑93 (MQ=255)

              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           <  1:6959320/93‑1 (MQ=255)

              gATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                                           <  1:167350/93‑1 (MQ=255)

                tCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                                         >  2:1042127/1‑93 (MQ=255)

                   ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGaa                                                                      >  2:5273921/1‑93 (MQ=255)

                   ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGaa                                                                      >  1:7803426/1‑93 (MQ=255)

                   ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGaa                                                                      <  2:2081959/93‑1 (MQ=255)

                   ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGaa                                                                      >  2:2825449/1‑93 (MQ=255)

                   ccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGaa                                                                      <  2:3936042/93‑1 (MQ=255)

                    cTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAt                                                                     <  2:5182918/93‑1 (MQ=255)

                     tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAtt                                                                    >  2:7368855/1‑93 (MQ=255)

                       tACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTg                                                                  >  1:5352901/1‑93 (MQ=255)

                        aCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGc                                                                 <  1:5001629/93‑1 (MQ=255)

                           ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgag                                                              >  2:4629506/1‑93 (MQ=255)

                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggagg                                                         >  1:2506994/1‑93 (MQ=255)

                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggagg                                                         >  2:5616096/1‑93 (MQ=255)

                                 ccGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGggagga                                                        >  1:1254251/1‑93 (MQ=255)

                                  cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGaa                                                       <  1:557617/93‑1 (MQ=255)

                                   gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAat                                                      <  2:4911317/93‑1 (MQ=255)

                                   gagaATTGCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAat                                                      <  1:796863/93‑1 (MQ=255)

                                     gaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAatat                                                    >  2:812188/1‑93 (MQ=255)

                                      aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATc                                                   <  1:3948273/93‑1 (MQ=255)

                                      aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATc                                                   <  2:7006276/93‑1 (MQ=255)

                                       aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCt                                                  <  2:2694978/93‑1 (MQ=255)

                                        tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                                                 >  2:1175774/1‑93 (MQ=255)

                                        tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                                                 >  1:2672828/1‑93 (MQ=255)

                                          tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGc                                               <  2:3344045/93‑1 (MQ=255)

                                          tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGc                                               <  2:7982984/93‑1 (MQ=255)

                                           ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAACGCGTAACGTCGCCGCAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCa                                              >  2:4440970/1‑93 (MQ=255)

                                            cTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAg                                             >  1:5981490/1‑93 (MQ=255)

                                            cTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAg                                             <  1:3843516/93‑1 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCGAGGCGTAACGTCGCACAAGAGACAACGCGAATTTGCGAGGGAGGAATATGTTGccgc                                            >  1:8374046/1‑90 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCTTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGa                                            <  1:3130281/93‑1 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGa                                            <  2:8587327/93‑1 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTCGCGAGGGAGGAATATCTTGCAGa                                            >  2:4763681/1‑93 (MQ=255)

                                              gcgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGAc                                           >  1:7507886/1‑93 (MQ=255)

                                              gcgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGAc                                           <  1:3427407/93‑1 (MQ=255)

                                                gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                         >  1:2941454/1‑93 (MQ=255)

                                                gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                         <  2:8193202/93‑1 (MQ=255)

                                                gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                         <  2:84616/93‑1 (MQ=255)

                                                 cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGtt                                        <  2:3788138/93‑1 (MQ=255)

                                                  tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                       <  2:2172529/93‑1 (MQ=255)

                                                   ttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCa                                      <  2:3408404/93‑1 (MQ=255)

                                                   ttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCa                                      <  2:3047295/93‑1 (MQ=255)

                                                   ttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCa                                      <  2:8806162/93‑1 (MQ=255)

                                                   ttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCa                                      <  2:1598338/93‑1 (MQ=255)

                                                       cAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaa                                  >  2:5486420/1‑93 (MQ=255)

                                                        aaTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaaa                                 >  1:2409394/1‑93 (MQ=255)

                                                         aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAAc                                <  2:5854185/93‑1 (MQ=255)

                                                          tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa                               >  2:3224783/1‑93 (MQ=255)

                                                           aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                              <  2:6016291/93‑1 (MQ=255)

                                                            aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                             <  1:1569423/93‑1 (MQ=255)

                                                            aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                             <  1:8307714/93‑1 (MQ=255)

                                                             aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCTACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            <  1:5139006/93‑1 (MQ=255)

                                                             aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            <  2:7287086/93‑1 (MQ=255)

                                                              aTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcg                           >  2:2364017/1‑93 (MQ=255)

                                                               tCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcgc                          <  2:8067364/93‑1 (MQ=255)

                                                                  cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaac                       >  2:6286432/1‑92 (MQ=255)

                                                                  cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                       <  1:1858086/93‑1 (MQ=255)

                                                                  cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                       >  1:4746401/1‑93 (MQ=255)

                                                                   ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCACACGTTCAGAAAACAGGCGCaaaa                      >  2:6553269/1‑92 (MQ=255)

                                                                     gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                    >  1:6140596/1‑93 (MQ=255)

                                                                      cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                   <  1:8514535/93‑1 (MQ=255)

                                                                      cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                   <  1:2288910/93‑1 (MQ=255)

                                                                       aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTg                  <  2:1380347/93‑1 (MQ=255)

                                                                         tGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGcc                <  2:3508480/93‑1 (MQ=255)

                                                                           gcGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCaa              <  2:7617574/92‑1 (MQ=255)

                                                                             ggCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAcc            <  2:1667753/93‑1 (MQ=255)

                                                                              gCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAccc           >  1:7189156/1‑92 (MQ=255)

                                                                              gCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcg           >  1:100083/1‑93 (MQ=255)

                                                                              gCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcg           >  1:5530567/1‑93 (MQ=255)

                                                                              gCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcg           >  1:3716382/1‑93 (MQ=255)

                                                                              gCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcg           >  1:6525611/1‑93 (MQ=255)

                                                                              gCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcg           <  1:1484749/93‑1 (MQ=255)

                                                                                aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcgcg         >  1:2783623/1‑93 (MQ=255)

                                                                                aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcgcg         >  2:6788306/1‑93 (MQ=255)

                                                                                 aGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc        >  2:6128820/1‑93 (MQ=255)

                                                                                 aGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc        >  1:8210270/1‑93 (MQ=255)

                                                                                 aGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc        <  1:8716525/93‑1 (MQ=255)

                                                                                  ggCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCt       <  1:766944/93‑1 (MQ=255)

                                                                                  ggCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCt       <  1:5259472/93‑1 (MQ=255)

                                                                                  ggCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCt       >  1:4068972/1‑93 (MQ=255)

                                                                                    cGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTcc     <  1:3311015/93‑1 (MQ=255)

                                                                                     gTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCa    >  2:3633407/1‑93 (MQ=255)

TCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATT  >  NC_000913/3961498‑3961677

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