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%	958.5 / NA	296	R86T (AGG→ACG)	gppA	guanosine pentaphosphatase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (146/150); total (146/150)

Only 100 of 299 total aligned reads displayed.

TCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCATT  >  NC_000913/3961498‑3961677

tCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTc                                                                                         <  2:297798/93‑1 (MQ=255)

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

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

   tctTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGcc                                                                                      <  2:1528628/93‑1 (MQ=255)

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

     ttCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGt                                                                                    <  1:506746/93‑1 (MQ=255)

     ttCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGt                                                                                    >  1:1551842/1‑93 (MQ=255)

      tCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTa                                                                                   <  2:3351277/93‑1 (MQ=255)

      tCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTa                                                                                   <  1:5610834/93‑1 (MQ=255)

        aCCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTatc                                                                                 >  2:5772719/1‑91 (MQ=255)

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

        aCCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGc                                                                                 >  2:4988642/1‑93 (MQ=255)

         ccGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCg                                                                                >  2:4183773/1‑93 (MQ=255)

         ccGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCg                                                                                >  1:3942134/1‑93 (MQ=255)

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

             tGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                                            <  1:2306127/93‑1 (MQ=255)

             tGATCACCTGTACCGGACAACCGAGAATTTCATGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                                            <  2:4409494/93‑1 (MQ=255)

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

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

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

                  cccTGTACCGGACAACCGAGAATGTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGa                                                                       <  2:4897481/92‑1 (MQ=255)

                    cTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAt                                                                     >  1:3476593/1‑93 (MQ=255)

                     tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCGCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAtt                                                                    <  2:2282738/93‑1 (MQ=255)

                     cgTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAtt                                                                    <  2:5833201/92‑1 (MQ=255)

                          cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCga                                                               <  2:2673954/93‑1 (MQ=255)

                           ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgag                                                              >  2:2168607/1‑93 (MQ=255)

                            gACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagg                                                             <  1:4899159/93‑1 (MQ=255)

                            gACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagg                                                             <  2:521566/93‑1 (MQ=255)

                             aCAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgaggg                                                            <  2:5652507/93‑1 (MQ=255)

                              cAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgaggga                                                           <  1:6013923/93‑1 (MQ=255)

                               aaCCGAGAATTACCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCATGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggag                                                          <  1:3383399/93‑1 (MQ=255)

                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggagg                                                         >  1:4740176/1‑93 (MQ=255)

                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggagg                                                         >  1:1062889/1‑93 (MQ=255)

                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggagg                                                         >  2:4496364/1‑93 (MQ=255)

                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggagg                                                         >  2:4763561/1‑93 (MQ=255)

                                aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCgagggagg                                                         <  1:1312966/93‑1 (MQ=255)

                                 ccGAGAATTTCCTGCGCTTTGGCGATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGggagga                                                        <  1:5793298/93‑1 (MQ=255)

                                 ccGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGggagga                                                        >  2:3384833/1‑93 (MQ=255)

                                  cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGaa                                                       >  2:3737498/1‑93 (MQ=255)

                                   gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAat                                                      >  1:2984421/1‑93 (MQ=255)

                                    agaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAata                                                     <  1:2215403/93‑1 (MQ=255)

                                    agaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAata                                                     <  2:1040489/93‑1 (MQ=255)

                                    agaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAata                                                     >  1:1261812/1‑93 (MQ=255)

                                     gaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAatat                                                    >  1:4354510/1‑93 (MQ=255)

                                      aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATc                                                   <  1:1955922/93‑1 (MQ=255)

                                      aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATc                                                   <  1:5921221/93‑1 (MQ=255)

                                       aattCCCGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                  <  2:6285886/91‑1 (MQ=255)

                                       aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                  <  1:2747541/93‑1 (MQ=255)

                                       aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCt                                                  >  2:198652/1‑93 (MQ=255)

                                        tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCtt                                                 <  1:2919646/93‑1 (MQ=255)

                                        tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCtt                                                 <  1:2995844/93‑1 (MQ=255)

                                        tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCtt                                                 <  2:2899249/93‑1 (MQ=255)

                                         ttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTg                                                <  2:1199237/93‑1 (MQ=255)

                                           ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCa                                              <  1:647437/93‑1 (MQ=255)

                                            cTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAg                                             <  1:5224338/93‑1 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGa                                            >  2:5582653/1‑93 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGa                                            <  1:1127509/93‑1 (MQ=255)

                                             tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGa                                            <  2:6078206/93‑1 (MQ=255)

                                                gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                         <  2:234541/93‑1 (MQ=255)

                                                gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                         >  1:4333461/1‑93 (MQ=255)

                                                gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGt                                         <  2:3154759/93‑1 (MQ=255)

                                                 cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGtt                                        <  1:4033346/93‑1 (MQ=255)

                                                  tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                       <  2:1652106/93‑1 (MQ=255)

                                                  tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                                       <  2:206926/93‑1 (MQ=255)

                                                    tGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAg                                     <  1:1657579/93‑1 (MQ=255)

                                                     ggCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGc                                    >  2:1304104/1‑93 (MQ=255)

                                                      gCAATAAAATCACCCGCATTAACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                                   >  2:2770607/1‑93 (MQ=255)

                                                      gCAATAAAATCACCCGCATGGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                                   <  1:3536133/93‑1 (MQ=255)

                                                       cAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaa                                  >  2:3813722/1‑93 (MQ=255)

                                                          tAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGATCAGCAAACa                               >  2:4610204/1‑93 (MQ=255)

                                                           aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                              <  1:1240886/93‑1 (MQ=255)

                                                           aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg                              <  2:4707840/93‑1 (MQ=255)

                                                            aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg                             >  1:3227633/1‑93 (MQ=255)

                                                             aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            >  1:4134281/1‑93 (MQ=255)

                                                             aaTCACCCACATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgc                            >  2:6152803/1‑93 (MQ=255)

                                                              aTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcg                           >  2:1434422/1‑93 (MQ=255)

                                                               tCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcgc                          <  2:624048/93‑1 (MQ=255)

                                                                  cccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                       <  1:3681522/93‑1 (MQ=255)

                                                                  cccGCATTGACGGCAAGGCGTAACGTCGCCGCAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa                       >  1:6174711/1‑93 (MQ=255)

                                                                   ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAAc                      <  2:788504/93‑1 (MQ=255)

                                                                     gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                    >  1:2808152/1‑93 (MQ=255)

                                                                     gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt                    <  2:902030/93‑1 (MQ=255)

                                                                      cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAtt                   >  1:3726684/1‑93 (MQ=255)

                                                                       aTTGACGGCAAGGTGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTg                  <  1:3602025/93‑1 (MQ=255)

                                                                       aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTg                  <  1:2115433/93‑1 (MQ=255)

                                                                        ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGc                 >  2:5271400/1‑93 (MQ=255)

                                                                        ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGc                 >  1:2472819/1‑93 (MQ=255)

                                                                          gACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCa               >  2:4118144/1‑93 (MQ=255)

                                                                            cGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAc             <  1:5529177/93‑1 (MQ=255)

                                                                              gCAAGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAAccc           >  1:4526028/1‑92 (MQ=255)

                                                                                aaGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACcgcg         >  1:2611813/1‑93 (MQ=255)

                                                                                 aGGCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCgcgc        >  2:4353848/1‑93 (MQ=255)

                                                                                  ggCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCt       <  1:5653751/93‑1 (MQ=255)

                                                                                   gCGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTc      <  1:5016678/93‑1 (MQ=255)

                                                                                    cGTAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTcc     >  1:5881667/1‑93 (MQ=255)

                                                                                      tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGATGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTACAt   >  2:540936/1‑93 (MQ=255)

                                                                                      tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt   <  1:3767608/93‑1 (MQ=255)

                                                                                      tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt   <  2:254453/93‑1 (MQ=255)

                                                                                      tAACGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAt   >  1:5348477/1‑93 (MQ=255)

                                                                                       aaCGTCGCCGTAGCGACAACACGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATTGCCAACCGCGCTCCAtt  <  1:3855345/93‑1 (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: `‑`